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ABSTRACT 

In  Part  I  a  semiempirical  approach  is  used  to  predict  fire 
spread  beyond  the  area  directly  ignited  by  thermal  radiation  from 
nuclear  weapons.  Mathematical  models,  both  stochastic  and  deter¬ 
ministic,  describe  the  progress  of  fires  in  two-dimensional  or  one¬ 
dimensional  space.  Application  of  each  model  to  appropriate  cases 
is  discussed.  Empirical  data  needed  for  evaluation  of  parameters 
are  specified  and  methods  for  acquiring  these  data  suggested.  Ob¬ 
served  data,  accumulated  over  many  years  in  records  of  past  fires, 
have  proved  valuable  in  determining  some  of  the  parameters.  The 
remaining  parameters  require  further  observed  data.  At  present 
approximate  prediction  can  be  made  by  use  of  a  specially  designed 
version. 

Part  II  presents  the  results  of  a  statistical  study  on  ob¬ 
served  rate  of  spread  data  and  discusses  a  number  of  specific 
problems  that  must  be  worked  out  before  the  method  can  be  used 
for  assessing  the  fire  damage  from  nuclear  attacks. 
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Section  1 
INTRODUCTION 


PREVIOUS  STUDY 

Mass  fires  can  produce  staggering  damage  both  ir.  wartime  and 
in  peacetime.  The  advancement  of  nuclear  weapons  has  introduced 
a  powerful  means  for  destroying  large  areas  by  direct  ignition  and 
subsequent  fire  spread.  The  great  consequences  of  mass  fires  fol¬ 
lowing  nuclear  attacks  have  stimulated  considerable  efforts  in  in¬ 
vestigating  various  aspects  of  the  behavior  of  mass  fires,  despite 
their  highly  complex  nature. 

One  of  these  efforts  was  a  small  study  performed  at  URS 
Corporation  (formerly  United  Research  Services)  in  1960  under  the 
sponsorship  of  the  U.S.  Forest  Service  (Contract  No.  12-11-005“ 
21911).  This  study  had  two  major  objectives,  namely,  to  develop 
a  simple  calculation  method  for  predicting  the  total  burnout  area 
following  a  nuclear  attack  on  the  basis  of  the  U.S.  Forest  Service 
1957  Damage  Assessment  Study  and  to  explore  the  possibilities  of 
developing  improved  models  for  predicting  the  total  area  of  spread. 

NEEDS  FOR  IMPROVED  MODEL 

This  study  indicated  that  significant  improvements  could  still 
be  made  in  fire  spread  prediction  models.  It  was  also  concluded 


*  Jewell,  W.  S.  and  A.  B.  Willoughby,  A  Study  to  Analyze  and  Improve 
Procedures  for  Fire  Damage  Assessment  Following  Nuclear  Attack. 

BRC  167-1,  Parts  I  and  II,  Final  Report,  October  1960  (AD  286  943). 
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that,  in  addition  to  predicting  the  ultimate  burnout  area,  future 
studies  should  attempt  to  describe  the  dynamics  of  fire  spread. 

For  gross  damage  assessment  in  simulated  or  actual  nuclear  attacks, 
there  seems  to  be  a  real  need  for  information  about  the  time- 
dependent  behavior  of  fires,  particularly  in  estimating  personnel 
casualties,  since  the  population  distribution  will  presumably  not 
be  static  in  the  early  postattack  period.  Such  information  is 
also  useful  at  the  local  level  for  planning  evacuation  routes  and 
for  deciding  whether  fire  fighting  is  feasible  in  protecting  crit¬ 
ical  installations  and  shelters.  When  the  use  of  a  model  is  ex¬ 
tended  to  peacetime  forest  or  urban  fires,  it  is  obvious  that  the 
knowledge  of  fire  dynamics  is  even  more  critical  since  it  permits 
rational  use  of  manpower  and  equipment  in  fire  fighting. 
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Section  2 

OBJECTIVES  AND  SCOPE 

DEFINITION  OF  MAIN  OBJECTIVE 

The  objective  of  the  present  study  is  to  develop  an  improved 
prediction  method  describing  the  extent  of  fire  spread  beyond  the 
initial  area  ignited  by  nuclear  attack  and  to  specify  the  data  re 
quired  for  use  with  the  method  and  ways  of  acquiring  these  data. 
Thus,  no  effort  is  devoted  to  prediction  of  the  initially  ignited 
area  as  a  function  of  weapon  characteristics,  meteorological  con¬ 
ditions,  and  fuel  characteristics,  a  problem  which  is  essential 
in  the  assessment  of  fire  damage  from  nuclear  attacks  but  which 
has  been  appropriately  considered  as  a  separate  task. 

SPECIFICATION  OF  REQUIRED  DATA 

Since  the  models  are  intended  for  both  present  and  future 
uses,  their  implementation  must  be  immediately  practical  with  the 
existing  input  data  and  be  capable  of  subsequent  improvement  as 
more  and  better  data  become  available.  Therefore,  after  models 
have  been  selected,  the  required  input  data  must  be  indicated  and 
procedures  for  gathering  existing  data  for  immediate  implementa¬ 
tion  and  for  future  acquisition  of  more  extensive  data  must  be 
suggested. 
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Section  3 

OUTPUT  REQUIREMENTS 


SPECIFIC  APPLICATIONS 


The  type  of  fire  spread  model  and  the  method  of  implementa¬ 
tion  must  be  determined  by  the  desired  output  information  and, 
consequently,  by  the  particular  uses  intended  for  the  model.  The 
general  purpose  of  fire  spread  models  in  civilian  defense  is  to 
provide  estimates  of  damage  to  the  nation’s  resources  from  the 
fire  effects  of  a  massive  nuclear  attack.  These  estimates  are 
needed  primarily  for  preattack  planning,  although  there  also  ex¬ 
ists  a  requirement  for  early  postattack  damage  assessment.  Pre¬ 
attack  and  postattack  applications  can  be  performed  at  three  geo¬ 
graphical  levels  as  follows : 


Preattack  Planning 

•  National.  To  help  assess  the  nation’s  ability  to  survive 
massive  nuclear  attack. 

•  Regional.  To  provide  a  basis  for  various  preattack  counter 
measures  such  as  stockpiling,  protection  of  critical  re¬ 
sources,  and  development  of  postattack  survival  and 
recovery  plans. 

•  Local.  To  provide  a  basis  for  local  preattack  thermal  coun 
termeasures  and  development  of  specific  survival  plans. 


Postattack  Indirect  Damage  Assessment 

•  National.  To  help  provide  early  assessment  of  the  nation's 
overall  survival  posture. 

•  Regional.  To  help  select  one  of  several  postattack  sur¬ 
vival  plans. 
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•  Local,  To  provide  a  basis  for  specific  actions  by  the 
population  and  local  civil  defense  units, 

REQUIRED  LEVELS  OF  DETAIL  IN  THE  OUTPUT 

It  is  clear  that  each  type  of  application  calls  for  a  differ¬ 
ent  level  of  detail  in  the  output  information.  National  uses  seem 
to  require  only  gross  information  —  for  example,  the  fraction  of 
each  resource  still  available  after  an  attack  of  hundreds  of  nu¬ 
clear  weapons.  Since  primary  concern  is  with  the  average  result 
for  a  large  number  of  weapons,  a  great  deal  of  uncertainty  in  pre¬ 
dicting  damage  resulting  from  each  individual  weapon  can  be  toler¬ 
ated,  provided  systematic  errors  are  not  serious. 

For  regional  uses,  a  somewhat  greater  level  of  detail  is 
needed.  For  example,  to  decide  where  and  how  much  to  stockpile 
certain  critical  supplies  demands  reasonably  accurate  knowledge 
of  the  expected  damage  averaged  over  areas  of  a  few  hundred  miles 
on  a  side.  Local  uses  would  require  by  far  the  greatest  level  of 
detail.  The  fire  spread  for  each  individual  nuclear  detonation 
must  be  predicted  with  fairly  good  accuracy. 

In  the  original  planning  for  the  present  study,  it  was  assumed 
that  models  for  national,  and  regional  application  were  most  ur¬ 
gently  needed.  Hence,  most  effort  was  concentrated  in  these  types 
of  application, and  only  preliminary  considerations  were  given  to 
models  for  local  use. 

CONVERSION  OF  FIRE  SPREAD  INFORMATION  TO  DAMAGE  INFORMATION 

For  damage  assessment  purposes,  the  needed  information  is 
the  damage  by  fire  or  the  fraction  of  resources  destroyed  rather 
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than  the  fire  boundaries  per  se.  As  with  predictions  of  damage 
from  other  nuclear  weapons  effects,  fire  damage  can  be  derived 
from  fire  spread  information  by  assigning  vulnerability  factors 
to  various  resources .  These  factors  give  a  measure  of  the  frac¬ 
tion  of  resources  destroyed  or  the  probability  of  their  destruc¬ 
tion  if  the  area  is  passed  over  by  fire. 
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Section  4 

DESCRIPTION  OF  THE  PROBLEM 

DIRECTLY  IGNITED  AREA 

Fire  behavior  in  the  area  directly  ignited  by  nuclear  weapons 
would  constitute  a  useful  and  interesting  study.  By  combining 
data  on  intensity  of  thermal  radiation  with  a  fire  spread  model 
especially  designed  for  handling  the  spread  of  multiple  ignition 
centers,  it  would  be  possible  tc  estimate  the  size  and  contour  of 
the  directly  ignited  area,  to  predict  how  fast  the  individual 
fires  grow,  merge,  and  burn  out  and  how  weather  and  fuel  affect 
such  a  burning  process.  Although  such  a  study  does  not  constitute 
the  subject  matter  of  the  present  program,  it  seems  appropriate  to 
have  a  qualitative  picture  of  what  happens  in  this  area  shortly 
after  ignition. 

The  radius  of  the  directly  ignited  area,  i.e.,  the  distance 
beyond  which  the  thermal  energy  is  insufficient  to  ignite  kindling 
fuel,  depends  on  three  groups  of  variables:  the  weapon  character¬ 
istics  (including  yield  and  height  of  burst),  the  transmission 
characteristics  of  the  atmosphere,  and  the  ignition  characteris¬ 
tics  of  fuels. 

Although  contours  of  equal  thermal  energies  are,  in  general, 
concentric  circles  with  center  at  ground  zero  (Fig.  1),  the  shape 
of  the  initially  ignited  area  will  not  be  circular  unless  the 
kindling  fuels  are  uniformly  distributed,  which  is  not  necessarily 
the  case.  Furthermore,  the  ignition  radii  give  no  information  on 
the  density  of  ignition  centers  (number  of  ignition  centers  per 
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unit  area)  and  its  variation  with  distance  from  ground  zero.  In 
general,  however,  we  can  expect  the  density -of  ignition  centers 
to  decrease  with  increasing  distance  from  ground  zero.  This  is 
due  to  the  fact  that  the  thermal  intensity  decreases  with  increas¬ 
ing  distance.  Consequently,  a  decreasing  number  of  kindling  fuels 
actually  occurring  in  the  area  can  be.  ignited. 


The  evolution  of  the  initially  ignited  area  can  be  qualita¬ 
tively  described  by  three  successive  stages: 

(1)  Provided  sufficient  fuel  is  available,  some  q.f  the 
ignition  centers  spread  to  adjacent  fuels  and  build 
up  to  small  and  sustained  fires.  The  rest  go  out  due 

.o  their  tenuity  or  unfavorable  surroundings. 

(2)  In  the  next  stage,  these  small  fires  merge  together  to 
form  larger  fires ,  which  —  with  further  merging  —  may 
take  mass  fire  proportions.  The  time  for  merging 
depends  on  many  factors,  of  which  the  density  of  small 
fires  is  perhaps  the  most  important.  Since  this  den¬ 
sity  is  higher  near  the  center,  the  fires  in  the  central 
portion  generally  merge  and  burn  out  earlier  than  those 
near  the  periphery.  The  possibility  of  a  fire  storm 
comes  into  the  picture  at  this  stage.  Unfortunately, 

we  know  little  about  the  factors  which  produce  fire 
storms  and  how  fire  storms  modify  the  evolution  of  mass 
fires.  We  know,  however,  that  there  are  mass  fires 
without  genuine  fire  storms. 

(3)  In  the  absence  of  a  fire  storm,  a  burnout  zone  first 
appears  at  the  central  portion  and  expands  toward  the 
periphery.  The  third  stage  begins  when  the  burning  area 
is  reduced  to  a  more  or  less  continuous  ring  which  keeps 
moving  outward.  With  a  fire  storm  the  behavior  is  less 
clear.  However,  if  one  assumes  that  the  fire  storm  does 
not  extinguish  the  fire,  it  seems  reasonable  to  again 
assume  a  ring  fire  after  the  fire  storm  has  died  out. 

It  is  the  spread  of  this  ring  that  we  undertook  to  study 
in  the  present  program. 
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RELATIONS  BETWEEN  FIRE  SPREAD  PHENOMENA 

Because  of  the  intimate  relations  betV7een  various  character¬ 
istics  of  fire  spread,  the  study  of  burnout  area  (required  for 
estimating  fire  damage  to  resources)  cannot  be  separated  from  that 
of  other  quantities  equally  fundamental  in  understanding  fire 
behavior.  Such  quantities  are  rate  of  spread,  evolution  of  heat 
and  flame,  size  of  burning  area  at  a  given  time,  and  percentage 
of  burned  fuel  in  a  general  burnout  area. 

Rate  of  spread  is  naturally  derivable  from  any  model  that 
represents  fire  spread  both  in  time  and  in  space.  Heat  and  flame 
output,  per  unit  burning  fuel,  are  determined  primarily  by 
the  type  of  fuel  and  possibly  by  its  moisture  content.  Since 
higher  heat  output  raises  the  temperature  of  exposed  fuel  more 
rapidly,  the  rate  of  spread  is  expected  to  depend  on  heat  output 
and,  hence,  on  fuel  type. 

The  instantaneous  area  of  the  burning  fire  may  play  an  impor¬ 
tant  role.  The  ignition  rate  and  the  spread  rate  depend  on  the 
total  heat  flux  (plus  embers  and  firebrands),  which  for  a  given 
set  of  fuel  and  weather  conditions  is  determined  by  the  area  of 
the  burning  fire.  This  total  flux  will  become  independent  of  the 
burning  area  only  when  the  fire  has  reached  a  dimension  such  that 
the  heat  from  the  rear  portion  cannot  reach  the  unignited  fuel  in 
front  of  the  fire.  Furthermore,  a  very  large  burning  fire  may 
conceivably  modify  its  environment  (create  its  own  weather)  and 
consequently  -behave  in  a  manner  not  predictable  on  the  basis  of 
conditions  initially  specified.  Therefore  it  is  important  to 
know  how  large  a  burning  fire  is  at  a  given  moment  and  under  what 
conditions  a  big  fire's  behavior  is  influenced  by  new  weather 
factors  generated  by  itself. 
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The  percentage  of  fuel  ultimately  burned  out  is  determined 
by  the  degree  of  homogeneity  of  the  fuel.  Since  fuel  distribu¬ 
tion  is  an  important  factor  in  setting  up  and  applying  a  fire 
spread  model,  it  is  probably  worthwhile  to  elaborate  upon  this 
topic  somewhat  further. 

An  area  containing  a  single  type  of  fuel  with  uniform  dis¬ 
tribution  is  either  completely  burned  out  or  completely  unburned. 
Uniform  distribution  is  frequently  encountered  in  wildland  fuel. 

A  large  area  is  normally  made  up  of  sections  of  various  sizes, 
each  with  a  rather  uniform  and  continuous  fuel  distribution.  For 
urban  areas  the  degree  of  homogeneity  is  rather  low.  With  the 
possible  exception  of  residential  tracts,  structural  materials, 
sizes,  and  spacings  may  vary  considerably.  Therefore  some  struc¬ 
tures  may  survive  a  fire  which  destroys  others. 

As  output  information  from  the  model,  the  burnout. percentage 
depends  on  the  level  of  detail  of  the  input  data.  To  illustrate 
this  point,  suppose  we  consider  predicting  fire  spread  in  a  small 
area,  say  20  to  30  square  miles,  for  local  damage  assessment  pur¬ 
pose.  This  special  application  requires  that  fuel  characteristics 
be  determined  with  sufficient  geographic  details.  For  example,  a 
set  of  fuel  characteristics  might  be  specified  for  each  unit  of 
one-hundredth  of  a  square  mile.  This  unit  would  then  constitute 
what  might  be  called  the  resolution  of  the  model  output.  Any 
prediction  regarding  the  burnout  event  will  refer  to  the  unit  as 
a  whole, and  the  burnout  percentage  within  each  unit  cannot  be 
given. 

When  very  large  areas,  such  as  major  sections  of  the  country 
or  the  whole  country  itself,  are  considered,  two  factors  will 
limit  the  level  of  detail  in  the  output  information,  namely,  the 
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size  of  the  unit  area  (or  cell)  and  the  size  of  areas  for  which 
fuel  characteristics  can  be  specified.  The  maximum  number  of 
cells  and,  hence,  the  minimum  size  of  the  cell  is  determined  by 
the  amount  of  computation  that  can  be  afforded.  A  similar  eco¬ 
nomical  reason  will  restrict  the  specification  of  input  data  to 
a  single  set  of  fuel  conditions  for  each  section  (herein  called 
"big  square")  containing  a  large  number  of  cells.  A  big  square 
might  be  20  or  25  miles  on  each  side.  Thus  the  model  will  tell 
whether  a  cell  has  burned  out  or  give  the  probability  of  the 
burnout  event  but  cannot  tell  what  fraction  of  the  cell  has  burned 
out.  It  may  give  a  different  answer  for  each  cell  in  a  given  big 
square,  but  the  difference  is  due  to  locations  of  the  cells  rela¬ 
tive  to  the  origin  of  the  fire  rather  than  to  variations  in  fuel 
characteristics , 

FIRE  SPREAD  VARIABLES 

Many  factors  of  different  nature  combine  their  action  to 
determine  whether  a  fire  will  spread,  how  fast  and  how  far  it  will 
spread,  and  how  large  the  burning  area  will  be  at  a  given  moment. 
These  factors  are  called  fire  spread  variables  and  may  be  grouped 
into  uhree  general  classes:  fu^l  variables,  weather  variables, 
and  topography  variables.  We  will  now  attempt  to  list  in  each 
class  all  variables  that,  according  to  our  experience  and  judge¬ 
ment,  may  influence  fire  spread  and  to  qualitatively  indicate  how 
this  influence  is  brought  about. 

Fuel  Variables 


Fuel  Type 

Some  general  fuel  types  are:  grass,  brush,  and  timber  for 
wildland  areas,  industrial,  commercial,  and  residential  s  true  tur  e  s 
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for  urban  areas.  Each  general  type  may  include  a  number  of 
varieties.  It  is  rather  obvious  that  fuel  types  may  not  be 
equally  igni table  and  when  burning  may  have  different  igniting 
capability. 

Fuel  Buildup 

This  is  the  thickness  of  the  fuel  layer  as  measured  by  the 
number  of  stories  in  urban  fuel  or  the  fuel  height  in  wildland 
fuel.  Fuels  with  high  buildup  probably  have  longer  burning  time 
and  put  out  more  heat  and  more  flames  per  unit  burning  area. 
However,  the  effects  of  buildup  on  ignition,  spread  rate,  and 
long-range  spotting  are  not  so  obvious. 

Fuel  Density 

Our  terminology  assigns  this  term  to  the  fraction  of  total 
area  occupied  by  the  fuel.  Low  density  has  an  opposing  effect  on 
fire  spread  because  less  heat  is  released  per  unit  burning  area. 

If  the  fuel  distribution  is  specified  stochastically,  there  exists 
for  any  fuel  density  a  finite  chance  that  a  fire  will  spread 
indefinitely  through  the  area.  This  chance,  however,  becomes 
increasingly  small  as  the  density  decreases.  Therefore,  the 
mean  final  spread  distance  decreases  with  fuel  density. 

Spatial  Distribution  of  Fuel 

The  fuel  in  two  areas  of  the  same  fuel  density  may  not  be 
distributed  in  the  same  way.  The  fraction  of  fuel  area  may 
spread  out  and  mix  uniformly  with  the  nonfuel  fraction  or,  at 
the  other  extreme,  may  form  a  single  continuous  section.  Obvi¬ 
ously  the  latter  mode  of  distribution  greatly  enhances  the  proba¬ 
bility  of  stopping.  Fuel  distribution  may  have  other  types  of 
effects  that  are  less  apparent. 
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Fuel  Age 

Urban  fuels  probably  become  more  ignitable  after  reaching  a 
certain  age,  particularly  if  they  are  not  maintained  properly. 
Wildland  fuels  seem  to  decrease  in  overall  moisture  content  as 
they  grow  older  and,  consequently,  may  become  more  ignitable. 
Under  fuel  age  we  might  include  green  state  and  dead  state  of 
fuel.  The  age  effect  may  also  operate  through  fuel  fineness. 

Fuel  Fineness 


This  factor  contributes  to  ignitability  and  spread  rate.  It 
may  promote  intense  burning  and  reduce  the  burning  time. 

Weather  Variables 


Moisture  Content 

Although  it; certainly  has  an  effect  on  ignitability,  rate  of 
spread,  and  probability  of  stopping,  this  variable  is  complicated 
and  difficult  to  define.  Rural  fuel  usually  consists  of  a  mix¬ 
ture  of  light,  medium,  and  heavy  components,  and  each  component 
can  be  either  green  or  dead.  At  any  time,  the  moisture  content 
may  he  different  in  each  case,  because  the  time  for  fuel  moisture 
to  come  into  equilibrium  with  the  surroundings  depends  on  each 
component  and  its  state.  The  finer  the  fuel,  the  faster  it  reaches 
moisture  equilibrium.  Thus  for  very  fine  fuels,  the  current  rela¬ 
tive  humidity  might  be  a  good  measure  of  the  moisture  content. 

Heavy  dead  fuels  dry  out  slowly,  with  the  result  that  their  mois¬ 
ture  content  at  any  time  depends  on  the  precipitation  and  atmos¬ 
pheric  moisture  for  a  number  of  days  preceding  the  time  of  interest. 
Also,  in  larger  fuels  there  is  probably  a  moisture  gradient  from 
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the  surface  to  the  center.  Thus  factors  such  as  days  since  last 
rain,  relative  humidity,  and  temperature  history  may  be  important. 

Wind 


Wind  carries  the  heat  into  unignited  fuel,  supplies  oxygen 
to  the  burning  fuel,  and  promotes  direct  flame  contact  and  spotting. 
This  variable  appears  to  be  simple,  easily  defined  by  speed  and 
azimuth.  Complications  arise,  however,  because  of  the  possibility 
of  vertical  wind  gradients  and  modifications  of  wind  speed  and 
direction  by  local  topography  and  fuel. 

Temperature 

The  higher  the  atmospheric  temperature,  the  less  heat  is 
required  to  raise  the  fuel  to  the  ignition  point.  However,  vari¬ 
ations  in  ordinary  atmospheric  temperatures  do  not  significantly 
affect  heat  requirements  for  ignition  in  mass  fires.  Thus,  except 
for  extreme  cases  perhaps,  temperature  by  itself  may  play  only  a 
minor  role  in  fire  spread. 

Topography  Variables 

Slope 

Since  hot  gases  tend  to  go  up,  the  slope  of  the  fuel  area 
favors  heat  convection  either  in  the  direction  of  unignited  fuel 
or  in  the  opposite  direction,  depending  on  whether  the  fire  is 
spreading  up  or  down  slope.  Thus  slope  effect  is  similar  to  wind 
effect  in  both  mechanism  and  directionality.  However,  for  a 
given  area,  wind  varies  with  time,  whereas  slope  effect  varies 
with  space  only.  The  slope  effect  may  be  complicated  or  counter- 
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balanced  by  concomitant  phenomena,  such  as  rolling  of  burning  fuel 
or  changes  in  spotting  conditions. 

Altitude 

As  such,  altitude  should  play  no  role  in  fire  spread  unless 
it  is  associated  with  other  variables,  such  as  vertical  wind  gradi¬ 
ents,  fuel  types,  moisture  content,  etc. 

Relief 

By  this  term,  we  mean  the  combined  effect  of  the  inequalities 
of  land  surface,  their  sizes,  and  their  mode  of  distribution. 

With  the  same  total  area,  the  hills  can  be  large  and  few  or  small 
|  and  many.  They  can  also  be  scattered  throughout  a  region  or  closely 

drawn  together  in  one  section.  It  is  conceivable  that  these  factors 
|  have  a  bearing  tn  fire  spread,  though  it  is  difficult  to  speculate 

i 

*  on  how  important  they  are  and  in  what  direction  they  operate. 

Terms  such  as  "probably,” "perhaps ,"  and  "may”  have  been  used 
abundantly  in  the  preceding  paragraphs  to  emphasize  the  fact  that 
without  experimental  data,  our  knowledge  on  fire  spread  is  highly 
uncertain  and  purely  qualitative  at  best.  The  effect  of  a  vari- 
j  able  may  be  real  but  indiscernible  among  other  more  important  ones, 

|  or  it  may  remain  sensibly  constant  over  the  usual  range  of  fire 

spread  conditions.  A  possible  example  of  constant  variable  might 
be  the  fuel  build-up  iifi  wildland  areas.  For  a  given  type  of  fuel, 
the  variation  in  buildup  is  probably  not  significant.  Some  vari¬ 
ables  operate  in  opposite  direction  and  partially  or  completely 
cancel  each  other's  effect.  Others  are  not  strictly  independent 
but  are  frequently  associated  in  their  occurrence.  For  example, 
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fuel  type  may  be  associated  with  slope;  fuel  buildup  and  moisture 
content  with  fuel  fineness,  etc. 
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Section  5 

METHODS  OF  APPROACH 

THREE  METHODS  OF  APPROACH 

Possible  models  for  prediction  of  fire  spread  can  be  classi¬ 
fied  into  three  general  types:  purely  theoretical,  purely  em¬ 
pirical,  and  semiempirical* 

Purely  theoretical  models  would  be  characterized  by  an  attempt 
to  predict  fire  spread  on  the  basis  of  physical  laws,  without  re¬ 
sorting  to  observed  data  on  actual  fires.  This  approach  is  prob¬ 
ably  one  that  we  would  like  to  use  eventually;  however,  physical 
laws  needed  to  describe  various  phenomena  in  macroscopic  fires 
are  still  incomplete  or  nonexistent.  Most  of  the  existing  laws 
have  been  obtained  by  calculation  and  experimentation  under  highly 
simplified  and  rigorously  controlled  conditions,  which  are  rarely 
encountered  in  actual  cases.  When  applied  to  large-scale  phenom¬ 
ena,  these  laws  are  likely  to  yield  questionable  results.  Thus 
much  basic  research  on  macroscopic  fires  is  still  to  be  done  be¬ 
fore  a  purely  theoretical  approach  is  practical. 

A  purely  empirical  approach  would  be  to  analyze  data  from 
actual  fires  and  obtain  a  series  of  tables  or  graphs  showing  how 
a  given  quantity  depends  on  variation  of  one  or  more  controlling 
factors.  The  given  quantity  might  be  the  final  burnout  area, 
the  rate  of  spread,  or  some  fire  behavior  that  might  be  of  in¬ 
terest  in  the  future.  Prediction  of  potential  fires  could  be 
made  by  comparing  the  fuel  characteristics,  weather,  and  topog¬ 
raphy  conditions  with  these  tables  or  graphs  and  probably  by 
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doing  some  interpolation  or  extrapolation.  Because  of  the  complex¬ 
ity  of  the  phenomena  and  the  large  number  of  variables  involved,  it 
would  be  difficult  with  this  type  of  approach  to  organize  empirical 
data  logically  without  some  theoretical  backup.  Furthermore  the 
procedures  for  setting  up  such  a  model  and  for  using  it  are  too 
mechanical  to  yield  much  understanding  of  mass  fires. 

Selection  of  the  Semiempirical  Approach 

In  view  of  the  present  state  of  the  art,  a  semiempirical 
approach  seems  to  be  the  most  appropriate  and  has  been  adopted 
in  the  present  study.  In  this  approach,  the  fire  spread  process 
is  mathematically  represented  in  terms  of  a  mechanism  such  as 
suggested  for  similar  phenomena  in  epidemiology.  The  mechanism 
involves  some  basic  properties  which  may  depend  on  many  variables 
and  which  appear  in  the  mathematical  description  as  parameters. 
Parameters  must  be  provided  by  a  study  of  records  of  actual  fires 
or  by  direct  experimentation.  Thus  the  basis  for  semiempirical 
prediction  of  fire  spread  lies  both  in  the  body  of  experimental 
data  and  in  the  selected  mechanism. 

Treatment  of  Fire  Spread  Variables 

Detailed  procedures  in  mathematical  representation  are 
largely  determined  by  the  stochastic  or  deterministic  character 
of  the  variables.  At  certain  scales  of  time  and  space,  all  varia¬ 
bles  are  stochastic,  i.e.,  all  undergo  random  variations  of  such 
amplitudes  as  to  significantly  affect  the  process  of  fire  spread. 
Since  a  large  number  of  stochastic  variables  may  render  the  model 
excessively  complicated,  it  is  important  to  decide  which  variables 
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can  be  treated  as  deterministic  and  to  use  only  those  stochastic 
variables  required  in  each  particular  situation.  To  help  visual¬ 
ize  the  problem  let  us  qualitatively  represent  the  tendency  to 
fire  spread  as  a  function  of  time  T  and  space  X  (one  dimensional) . 
Tendency  to  fire  spread  is  an  undefined  but  unambiguous  quantity 
that  we  shall  denote  by  Z.  (See  Fig.  2.) 

At  a  given  T,  time- dependent  variables  (weather)  have  the 
same  values  throughout  a  rather  large  area.  Therefore  the  quan¬ 
tity  Z  will  vary  with  space  only,  i.e.,  with  space- dependent 
variables.  Its  variations  may  be  represented  by  a  curve  such 
as  A.  There  are  three  types  of  fluctuations:  small  ripples 
within  short  distances  in  the  range  of  20  —50  ft,  due  to  lack 
of  perfect  uniformity  in  a  given  fuel  type;  wider  fluctuations, 
spanning  distances  of  a  quarter  of  a  mile  or  more,  corresponding 
to  changes  in  general  fuel  type  or  slopes  (in  going  from  brush 
to  timber,  or  from  residential  area  to  industrial  area)  ;  and 
occasional  sharp  dips  to  zero  in  empty  spaces. 

At  a  given  distance  X,  the  quantity  Z  varies  with  time 
according  to  some  curve  such  as  B.  Fluctuations  are  of  two 
general  types:  high-frequency  fluctuations  due  to  unsteady 
wind  and  low-frequency  fluctuations  due  to  slow  changes  in 
moisture  content,  temperature,  and  long-term  value  of  wind 
velocity.  The  instantaneous  wind  velocity  is  very  unstable, 
even  for  periods  of  10  or  20  min.  The  effect  of  this  in¬ 
stability  on  fire  spread  is  depicted  on  curve  B  (Fig.  2)  by 
fluctuations  of  small  amplitude.  The  latter,  however,  should 
not  be  taken  too  seriously,  because  most  useful  quantities  in 
fire  spread  are  not  truly  instantaneous  values  but  rather 
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averages  over  short  periods  of  time.  Most  other  variables  in  the 
weather  category,  including  long-term  value  of  wind  velocity, 
change  with  a  much  slower  tempo.  The  Z- value  resulting  from  the 
combined  effect  of  these  variables  is  likely  to  hold  rather 
steady  for  intervals  of  6  to  12  hr. 

It  thus  appears  that  in  many  cases,  weather  variables  can  be 
treated  as  deterministic  quantities.  In  these  cases,  their  aver¬ 
age  values  would  be  specified  for  each  time  interval  (6,  12,  or 
24  hr) . 

The  X  variables,  on  the  contrary,  are  considerably  more  random, 
especially  in  urban  fuel  areas.  Deterministic  treatment  of  these 
variables  would  require  specification  of  any  variation  in  fuel 
type,  density,  buildup,  or  spacing  and  in  empty  areas.  Whether 
or  not  this  is  practical  depends  on  the  level  of  effort  devoted 
to  fuel  survey  and  numerical  c.omputation.  The  level  of  effort 
does  not  seem  to  be  excessive  if  fire  spread  prediction  is  re¬ 
stricted  to  an  area  of  25  or  30  square  miles.  For  much  larger 
areas,  the  effort  may  not  be  justified  by  the  detail  require¬ 
ment  in  the  output  information.  In  this  case,  the  probability 
distributions  of  space  variables  must  be  given  and  a  stochastic 
model  must  be  used. 
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Section  6 

MATHEMATICAL  MODELS 

CLASSIFICATION  OF  MODELS 

Figure  3  shows  two  families  of  models:  fire-front  and  fuel- 
state.  In  each,  a  number  of  specific  versions  have  been  derived 
from  the  parent  model.  Specific  versions  serve  a  variety  of  pur¬ 
poses:  illustrating  particular  features  of  the  parent  model, 
deriving  some  general  properties  which  could  not  be  conveniently 
studied  with  the  parent  model,  providing  more  convenience  in 
specific  applications. 

In  the  following  paragraphs,  each  model  or  specific  version 
will  be  formulated  in  some  detail,  and  the  more  important  results 
will  be  explicitly  indicated.  All  mathematical  derivations  will 
be  presented  in  Appendixes. 

FIRE  FRONT  MODELS 

The  fire  front  model  and  all  its  variations  treat  the  fire 
front  (or  the  entire  fire  itself)  as  a  special  random  walker  mov¬ 
ing  along  a  row  of  cells  or  small  square  areas.  In  each  short  time 
interval  dt,  the  random  walker  may: 

•  Die  or  stop  permanently  with  probability  Hfdt 

9  Move  one  cell  forward  with  probability  \dt 

•  Pause  or  stay  where  it  is  with  probability  1  -  (ji^+\)dt 
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FIRE  FRONT  MODELS 


FUEL  STATE  MODELS 


Fig.  3.  Two  Groups  of  Models 
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The  parameters  and  X  are  probability  rates  or  probabili¬ 
ties  per  unit  time  for  the  extinction  and  spread  processes.  This 
model  and  the  following  ones  describe  the  fire  spread  process  by 
using  only  two  or  three  parameters.  Consequently  the  latter  are  by 
necessity  complicated  functions  of  all  fire  spread  variables. 

Their  role  and  significance  will  be  discussed  later  in  this  report. 
For  the  present,  only  their  general  definitions  as  given  above 
are  required. 

The  mechanism  does  not  allow  the  fire  to  move  more  than  one 
cell  during  the  interval  dt.  The  underlying  assumption  can  be 
made  plausible  by  choosing  appropriate  cell  dimension  and  con¬ 
sidering  only  cases  in  which  long-range  spotting  does  not  occur. 

A  version  to  be  discussed  later  is  sufficiently  general  to  include 
the  spotting  process. 

Denoting  by  pn(t)  and  QR(t)  the  probability  that  at  time  t 
the  fire  front  is  at  cell  n  or  has  died  at  cell  n,  respectively, 
the  basic  equations  can  be  written  as  follows: 

BP  (t) 

“at  -  -<*  +  <» 

3P  (t) 

—  ^n-1^)  **"  ^f)^n(t)  n  =  2,  3,  .  .  .  (2) 

BQ  (t) 

«  HfPn(t)  n  =  0,  1,  2,  .  .  .  (3) 


(4) 
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The  solutions  are 
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where 


0  =  +  X 


Variations  of  Pn(t)  and  Qn(t)  are  illustrated  in  Fig.  4. 

As  t  approaches  infinity,  the  burnout  distances  follow  a 
geometric  distribution 


V">  -  T  III" 


The  mean  ultimate  burnout  distance  is 


n(°°) 


L  nQn(”) 


n=o 


cell  widths 


(7) 


(8) 


At  time  t,  the  mean  position  of  the  fire  front,  given  that  it 
is  still  alive,  is  given  by 

TnP  (t) 

np(t)  =  ^YJt)  =  Xt  cel1  widths  (9) 

Thus,  the  mean  velocity  of  the  burning  front  is  X  cell  widths 
per  unit  time.  Another  important  result  is  that  the  mean  lifetime 
of  the  fire  is  l/p^. 
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APPLICATION  TO  INHOMOGENEOUS  CELLS  (MODEL  IA) 


If  fuel  characteristics  are  specified  for  individual  cells, 

the  above  model  can  be  modified  by  taking  a  value  p  for  each 

cell  n  and  a  value  X  for  each  pair  of  cells  n,n+l.  The 

n,n+l 

equations  become 


aPo^> 

St 


~^o  +  X0pPo^ 


(10) 


SP  (t) 

St  =  Xn-1  ,nPn-l  ’  ^Xn,n+1  +  Mn^Pn^ 


(ID 


3<L(t) 


n 


at  '  =  ^nPn(t) 


(12) 


A  closed  expres 


O  JlKJLL  L.11CI  L-. 


p  n-1  X. 

Q  (ro)  =  — r - — -  n  - - 

xnv  '  p  +  X  ,i  .  |i.  +  X.  . 

n  n,n+l  i=o  i  i,i+1 


(13) 


Solutions  for  other  quantities  are  more  conveniently  per¬ 
formed  by  computer,  although  some  complicated  analytical  expres¬ 
sions  can  be  obtained. 

If  the  fuel  area  contains  several  uniform  regions  having 
different  fuel  characteristics,  the  basic  model  can  be  applied 
to  each  region,  and  the  resulting  distributions  can  be  adjusted 
to  match  each  other  at  the  discontinuities. 
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APPLICATION  TO  SPOTTING  (MODEL  IB) 


From  any  cell  n  the  fire  can  jump  to  cells  n+1,  n+2,  n+3,  .  .  . 
with  decreasing  probabilities  Xn  n+^dt,  Xn  n+2<^t>  Xn  n+3<^t’  •  •  •• 

If  all  cells  are  stochastically  identical,  a  single  value  applies 
and 

X01  =  X12  =  X23  =  •  '  *  =  X1 

X02  =  x13  =  x24  =  *  *  *  =  X2  (14) 

Xo,i  =  X  1,1+1  88  X2,i+2  =  *  *  *  Xi 

Calculation  of  Qn(°°)  can  be  performed  recursively  by  using 
the  relation 


x-,Q_  i  (“)  +  x9Q  9(“)  +  •  •  •  +  x  Q (°°) 

Q  M  =  B:1 - - - - s_ 2 -  (15) 

n 


Mf  + 


Si 


i=l 


Details  on  mathematical  derivations  are  presented  in 
Appendix  A. 

FUEL  STATE  MODELS 


The  burning  process  is  evidently  a  gradual  and  continuous 
one.  Unless  the  element  of  area  or  cell  under  consideration  is 
very  small  (a  few  square  feet),  instantaneous  and  uniform  igni¬ 
tion  of  the  entire  cell  is  of  rare  occurrence.  In  general  the 
cell  is  ignited  at  one  or  more  points  and  the  fire  builds  up  by 
intra-cell  spread.  The  chance  of  extra-cell  spread  (spread  out¬ 
side  the  cell  of  origin)  becomes  increasingly  greater,  then  begins 
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to  decline  after  the  fire  has  passed  its  peak  intensity.  This 
chance  reaches  the  zero  point  at  some  uncertain  moment  which 
coincides  with  or  precedes  total  extinction  of  the  cell.  Whether 
extra-cell  spread  actually  takes  place  depends  obviously  on  the 
condition  of  the  surrounding  cells  as  well  as  on  the  fire  intensity 
in  the  cell  under  consideration. 

This  continuous  train  of  events  is  somewhat  arbitrarily  divided 
into  three  distinct  stages  or  states,  the  unignited  state,  the 
flaming  state,  and  the  burnout  state.  The  flaming  state  is  defined 
as  one  capable  of  extra-cell  spread  regardless  of  whether  this  type 
of  spread  actually  takes  place  or  not,  and  it  is  left  undetermined 
as  to  what  phase  the  combustion  must  reach  or  what  fraction  of  the 
cell  must  burn  in  order  to  pass  to  F-state.  The  states  preced¬ 
ing  and  following  the  flaming  state  are  defined  as  unignited  and 
burnout  states,  respectively.  Thus  a  cell  which  has  been  ignited 
but  still  cannot  spread  fire  to  its  neighbors  is  by  definition  in 
the  unignited  state. 

The  concept  of  fuel  states  originated  in  the  striking  simi¬ 
larity  between  epidemic  breakout  and  fire  spread.  Many  phenomena 
in  epidemics,  such  as  susceptibility,  incubation,  infection, 
removal  of  victims  by  isolation  or  death,  dependence  of  infection 
on  intensity  of  epidemic,  age,  state  of  health,  population  den¬ 
sity,  etc.,  have  their  obvious  counterparts  in  fire  spread. 

The  fuel  state  models  attempt  to  predict  the  state  of  each 
cell  at  various  times  (deterministic  versions)  or  the  time -dependent 
probabilities  of  the  states  of  each  cell  (stochastic  versions). 
Models  II,  IIA,  IIB,  and  IIC  are  indicated  in  Fig.  3  and  will  be 
discussed  in  that  order. 
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Model  II 


Two -dimensionality ,  discrete  space  and  probabilistic  approach 
are  the  main  features  of  this  model.  The  probabilities  that  at 
time  t  cell  ij  is  in  the  unignited,  flaming,  and  burnout  states 
are  denoted  in  that  order  by  U^.(t),  F^(t),  and  (t)  satisfying 
the  relation: 

U. .  (t)  +  F .  .  (t)  +  B..(t)  =  1  (16) 

1 3  ij  ij 

Given  that  cell  ij  is  in  F-state,  the  probability  for  transition 

from  F-state  to  B-state  is  |i  dt  for  each  time  increment  dt.  The 

s 

decay  parameter  n  depends  on  some  or  all  fire  spread  variables. 

s 

However,  for  a  given  set  of  variable  values,  it  is  assumed  to  be 
a  constant,  independent  of  previous  happenings.  More  extensive 
discussions  will  be  devoted  to  this  and  other  parameters  in  sub¬ 
sequent  sections  of  this  report. 

If  cell  ij  is  in  U-state  and  its  neighbors  in  F-state,  the 
probability  that  ij  is  ignited  in  time  dt  is  denoted  by  A- .dt. 

(See  Fig.  5.)  Suppose,  for  the  moment,  that  the  cell  dimension 
has  been  chosen  so  that  a  burning  cell  can  ignite  the  immediate 
neighbor  cells  but  not  those  lying  further  away.  Then  an  unig¬ 
nited  cell  ij  can  be  ignited  by  one  or  more  of  its  8  immediate 
neighbors.  If  the  chance  of  ignition  by  each  neighbor  separately 
denoted  by  A^dt,  A^dt,  .  .  .,  Agdt,  then  the  chance  of  ignition  by 
one  or  more  neighbors,  i.e.,  A^dt,  is  given  by 

A..dt  =  1  -  (1  -  Andt) (1  -  A„dt)  .  .  .  (1- Audt) 

XJ  1  z  o 

=  (A1+A2  +  .  .  .  +  Ag)dt  -  (A1A2  +  A1A3  +  .  .  . )dt2  +  .  .  .  (17) 


I 
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2  3 

Neglecting  terms  in  dt  ,  dt  ,  .  .  . ,  we  have 


A^jdt  *  (A^  +  ^  +  .  .  .  +  Ag)dt 


(18) 


For  each  pair  of  igniting  cell  and  exposed  cell,  the  spread 
parameter  A  is  the  rate  of  the  ignition  probability.  Its  value 
is  determined  by  the  fuel  characteristics  of  the  two  cells,  the 
weather,  and  topography  variables. 

The  equations  of  the  model  can  be  written  as 


SF. . (t) 
St 


1 

U..(t)  )  F  (t) a 
lj  N  uv  u 

U=1 


SU..(t) 

.  .  ..1.1 ... 

St 


|i  F.  .  (t) 
s  IJ 


8 


u  *  1 ,  2 ,  .  .  . ,  8 


(19) 


■Ve>  Z  fu(t)Au 

u— 1 


(20) 


SB..(t) 

-  ^sFij  <*> 


(21) 


where  the  u’s  indicate  the  eight  cells  surrounding  i j . 

For  initial  conditions,  one  or  more  cells  may  be  specified  as 
being  in  F-state  and  all  others  in  U-state  at  t  *  0.  For  example 


F.  .(0)  =  1  for  i  =  0  j  =  0 
ij  J 

ILj  (0)  =  1  for  all  other  ij  combinations 


(22) 


Solutions  for  U..(t),  F..(t),  B..(t)  and  other  useful  quan- 
tities  must  be  obtained  by  numerical  or  machine  computation. 
However  approximate  expressions  can  be  obtained  if  we  consider 
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only  the  Initial  period  during  which  IL  ^ (t)  has  not  dropped  far 
below  1.  These  expressions  permit  some  qualitative  study  of  the 
system  before  extensive  computation  is  undertaken.  An  exact  solu¬ 
tion  for  B^,.  (»)  is  presented  in  Appendix  B. 

Detailed  Spread  Pattern 


In  the  above,  we  assumed  that  an  unignited  cell  can  be 
ignited  only  by  its  immediate  neighbors  with  spread  parameters 
^i »  A2>  .  .  .  Ag,  which  may  or  may  not  be  equal  or  symmetrical, 
since  they  depend  on  some  directional  factors  such  as  wind  and 
ground  slope.  If  the  cell  dimension  is  small  or  long-range  spot¬ 
ting  is  important,  a  cell  may  be  directly  ignited  by  more  remote 
neighbors.  In  this  case  the  spread  pattern  should  include  more 
than  one  layer  of  cells. 


The  equations  are  still  the  same,  except  the  summation  is 
extended  to  cover  all  possible  ignition  sources,  for  example 


aF..(t) 
ij  v 

at 


8*16+24 


=  U,,(f)  )  F  (t) A  -  u  F.. 

1J  us  U  S  lj 

U=1 


(t) 


(23) 


Linear  Version 


The  equations  become  somewhat  more  tractable  when  the  model 
is  applied  to  a  row  of  cells  in  which  lateral  and  back  spread  can 
be  ignored.  If  fire  can  spread  only  from  one  cell  to  the  next, 
we  can  write 


3F  (t) 
nv  ' 

at 


(24) 
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The  following  exact  solutions  can  be  obtained: 


"aBn-l(") 

B  (»)  -  1  -  e  a  - 

nN 

A/ns 

(25) 

■|Jst 

Fo(t)  -  e 

(26) 

~  -LUt  ‘  p?" 

F1(t)  -  ae  e  [ugt  +  a(l  -  e  )  +  2.2, 

-2nst  . 

(1  -  e  )  +  .  .  .| 

(27) 

B1(t)  =  ng  /  F1(T)dT 

(28) 

o 


Calculation  for  F^Ct)  is  most  conveniently  carried  out  with 
the  aid  of  a  computer. 

The  curves  for  F  (t)  and  B  (t)  are  somewhat  similar  to  those 

nx  nN 

for  Pn(t)  and  Q^(t)  in  the  fire  front  model.  (See  Fig.  4.) 

It  will  be  shown  analytically  in  Appendix  B  that  the  mean 
jump  time  is  1/a  and  the  mean  rate  of  spread  is  A  cell  widths  per 
unit  ti.ine.  This  result  has  been  confirmed  by  computer  calculation. 
For  this  purpose,  the  mean  location  n(t)  of  the  fire  is  defined  as: 

£  nF„(t) 

n ^  =  £F  (t) 


and  the  mean  velocity  is  given  by 


r  =  n(t  +  At)  -  n(t)  = 


At 


A  cell  widths 


/  O  A\ 


I 
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For  any  set  of  conditions  (i.e.,  any  value  of  a),  the  actual 
spread  distance  for  t  =  »  can  be  zero,  finite,  or  infinite.  There 
is  however  a  probability  for  each  distance,  and  the  mean  spread 
distance  D  given  by 


D  =  VbC)  (31) 

o 

is  a  characteristic  of  the  fuel  area  and  the  weather  during  the 
fire.  In  Appendix  B  it  will  be  shown  that  D  is  infinite  for 
a  £  1.  When  a  <  1,  the  mean  burnout  distance  is  to  a  good  approx¬ 
imation  given  by 


D  = 


£»„<•> 


(32) 


As  in  the  two-dimensional  model,  a  detailed  spread  pattern 
can  be  described  by  including  more  than  one  spread  parameter  A  and 
writing  the  equation  as 


f 

t 

K 

i  . 

I 

4  » 

I 

T 

I 

t 

■ 

I 

a 

■ 


a 


“aT"  -  V^k-l^l  +  Fn-2(t)A2  +  •••]-  MsFn(t)  (33> 

A 

Deterministic  Version 

I 

When  the  fuel  characteristics  and  weather  conditions  can  be 
specified  deterministically,  as  might  be  in  the  case  in  local  j| 

applications,  a  deterministic  model  is  desirable.  To  construct 
a  w.odel  of  this  type  we  start  with  the  following  considerations:  || 

•  The  ignition  rate  and,  consequently,  the  spread  rate  dX/'dt 

of  the  fire  front  depends  on  the  total  heat  flux  (including  J| 

I 
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embers  and  firebrands)  from  the  burning  zone.  The  total 
heat  flux,  in  turn,  depends  on  the  heat  released  per  unit 
burning  fuel  and  the  width  of  the  fire  front. 

•  When  the  width  of  the  fire  front  exceeds  a  certain  value 
W0,  the  heat  flux  from  the  rear  cannot  reach  the  unignited 
fuel  at  the  front.  WQ  depends  on  factors  such  as  wind  veloc¬ 
ity,  fuel  density,  etc. 


The  general  equation  is  then 


dX 

dt 


cp[W(t),k] 


(34) 


where  X  is  the  position  of  the  leading  edge  of  the  fire  front, 
W(t)  the  width  at  time  t,  k  a  spread  constant  and  a  measure  of 
the  heat  released  per  unit  burning  fuel,  and  cp  some  function  of 
W(t)  and  k.  The  exact  form  of  this  should  be  determined  empiri¬ 
cally,  perhaps  with  some  theoretical  backup.  However,  extensive 
work  in  this  direction  is  not  justified  at  this  stage,  and  a 
linear  relation  seems  to  be  a  fair  description  of  the  process. 
Thus 


f  -  kW(t) 
=  kW 

o 


from  which  we  derive 


W(t)  ^  wo 
W(t)  >  Wq 


W(t) 


dX(9) 

d9 


d9  = 


w(e)de 


t-T 


or 

-  kfW(t)  -  W(t-T)] 


(35) 


(36) 


(37) 


where  t  is  the  burning  time  of  the  fuel. 
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The  final  burnout  distance  X  is  given  by 

D  =  W  kT  (38) 

o  1  -  kT  ' 

as  derived  in  Appendix  E. 

Instantaneous  values  of  W(t)  and  dX/dt  have  been  obtained 
by  numerical  methods  and  plotted  in  Fig.  6. 

It  is  interesting  to  note  that  the  condition  for  infinite 
spread  in  this  model,  i.e.,  kT  ^  1,  is  equivalent  to  the  condi¬ 
tion  A/u  £  1  in  the  linear  stochastic  model.  Also  for  kT  <  1 
s 

the  final  spread  distance  X  is  a  most  exactly  equal  to  the  final 

mean  distance  in  that  model.  Also  k  and  t  are  equivalent  to  the 

stochastic  A  and  1/p  (the  mean  burning  time). 

s 

For  kT  <  1,  fire  front  width  decreases  steadily  to  zero 
regardless  of  whether  its  initial  value  W(0)  is  greater  or 
smaller  than  Wq.  The  velocity  however  remains  constant  when¬ 
ever  the  front  width  is  greater  than  Wq. 

For  kT  =  1,  the  front  width  decreases  or  increases  from 

its  initial  value  to  a  steady  value,  which  is  Wq.  The  velocity 

also  tends  to  the  constant  value  kW  . 

o 

For  kT  >  1,  the  front  width  increases  until  it  reaches  a 
constant  value  given  by  WQkT  while  the  velocity  increases  to 
kWQ  and  remains  unchanged  thereafter. 

If  information  on  W  and  the  function  cp  exists,  it  should 

o 

be  possible  to  determine  the  minimum  fuel  gap  required  to  stop 
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the  fire  and  to  investigate  the  effect  of  small  empty  spots 
(fuel  density)  on  the  velocity,  the  front  width,  and  the  final 
burnout  distance. 

Model  IIC 


The  problem  is  divided  into  two  parts:  first,  prediction 
of  whether  fire  can  spread  through  the  area  under  consideration; 
and  second,  if  it  does,  following  its  progress  in  time  and  space 
by  means  of  observed  data  on  rate  of  spread. 

Spread  Capability 

One  important  result  of  the  fuel  state  model  is  that  ex¬ 
cept  for  a  narrow  range  of  conditions,  i.e.,  of  A/p.  values, 

3 

a  fire  will  spread  a  negligible  distance  or  an  infinite  dis¬ 
tance  depending  on  the  particular  set  of  conditions.  Thus  to 
a  good  approximation,  spread  conditions  can  be  divided  into  two 
groups,  a  spread  group  and  a  no-spread  group. 

In  rural  fuel  weather  plays  a  dominant  role  in  making  a 

fire  go  or  stop  because  the  fuel  is  normally  continuous  and  the 

influence  of  commonly  occurring  topography  is  minor. 

« 

The  spread  capability  in  urban  fuel  depends  on  fuel  dis¬ 
tribution  more  strongly  than  on  weather  conditions.  Following 
are  two  alternative  methods  for  determining  spread  conditions 
in  urban  areas.  It  will  be  seen  that  for  each  type  of  urban 
fuel,  fire  spreads  indefinitely  if  the  fuel  density  (fraction 
of  area  occupied  by  fuel)  exceeds  a  certain  value. 

The  following  data  will  be  required  for  each  type  of  urban 
area:  the  spacing  distribution  a(x),  i.e.,  the  probability  that 
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an  adjacent  structure  is  at  a  distance  between  x  and  x  +  dx,  the 
average  structure  dimension  s  and  the  probability  PgR(x)  that  dur¬ 
ing  its  burning  time  a  burning  structure  S  ignites  a  structure  R 
at  distance  x. 

First  let  us  apply  the  linear  stochastic  model  to  an  urban 
area  divided  into  parallel  strips,  which  serve  as  cells.  It  is 
specified  that  Fq(0)  =  1 *  The  fraction  of  structures  burning 
in  cell  0  at  t  =  0  is  determined  by  the  spacing  distribution, 
size  distribution,  and  by  PCD(x)  but  does  not  have  to  be  known 
explicitly.  (See  p.  6-8.)  The  same  applies  to  the  fraction 
of  structures  ignited  at  any  time  in  other  cells. 

Let  P„  be  the  probability  that  cell  i,  if  ignited,  will 
ignite  cell  j  during  its  burning  time.  Since  each  cell  acts  as 
a  single  unit,  P. .  is  also  the  probability  that  cell  i,  if  ig- 
nited,  will  ignite  each  element  of  cell  j .  As  shown  in  con¬ 
nection  with  the  linear  stochastic  model,  P„  is  related  to  the 
parameter  a  (for  the  present  system)  by 

P..  =  1  -  e"a  (39) 

ij 

The  mean  final  spread  distance  is  infinite  if  a  £  1  or 

P. .  >  1  -  e  \  Thus  the  value  of  P. .  will  indicate  whether 
ij  i  J 

the  fire  will  spread  through  an  urban  area. 

P. .  can  be  calculated  as  follows.  Each  structure  in  j  is 
exposed  to  a  number  of  structures  in  i  (the  burning  cell),  with 
spacings  distributed  according  to  o>(x) .  For  an  average  pair, 
the  ignition  probability  must  be  weighted  according  to  this 
distribution:  ~ 

p  sr  =  J  PSR^X^  ^(x)dx  (40) 

o 
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To  estimate  the  average  number  of  structures  in  i  that  can 
be  seen  from  a  given  structure  in  j ,  we  line  up  structures  of 
average  size  s  on  a  half  circle  of  radius  s/2  +  x  where  x  is  the 
average  sparing.  Let  us  call  this  number  v.  Then  Pc.  is  given 
by 


Pij  i  -  (i  -  pSR) 


(41) 


The  requirement 


-1 


P. .  s  1  -  e 

for  infinite  spread  distance  can  now  be  expressed  as 

■  Vu  * e_1 

or 


PSR  *  1  -  6 


-<l/v) 


(42) 


(43) 


The  same  criterion  (a  £  1  or  PCE>  £  1  -  e  ^'v^)  Can  also 

oK 

be  obtained  by  an  alternative  derivation.  Fire  spread  is  pic¬ 
tured  as  proceeding  by  generations .  The  original  generation  is 
a  fire  burning  at  t  =  0,  which  may  be  a  straight  fire  front  or 
the  marginal  zone  of  an  area  fire  of  considerable  dimension. 
This  original  generation  ignites  a  group  of  structures,  called 
first  generation,  which  in  turn  ignites  a  second  generation  and 
so  on.  It  is  assumed  that  structures  in  generation  n  do  not 
participate  in  the  ignition  of  generation  n+2. 


The  structure  density  of  the  original  generation  (number  of 
structures  per  unit  area)  is,  for  convenience,  assumed  to  be  the 
same  as  the  structure  density  DQ  of  the  fuel  area.  The  original 
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generation  is  separated  from  the  rest  of  the  fuel  by  a  highly 
irregular  boundary.  As  an  average,  each  burning  structure 
faces  v  exposed  structures  and  each  exposed  structure  faces  the 
same  number  of  burning  structures.  We  have  derived  v  in  terms  of 
average  spacing  x  and  average  structure  dimension  s. 

The  mean  probability  that  a  structure  of  the  original  gener¬ 
ation  will  ignite  an  exposed  structure  is  Pg^  as  defined  above. 
The  probability  pQ  that  an  exposed  structure  is  ignited  by  one 
or  more  structures  of  the  original  generation  is 

Po  =  1  ■  (1  -  PSR)V  (44) 

which  has  been  also  denoted  by'P... 

The  structures  of  the  first  generation  are,  of  course, 
randomly  distributed,  and  it  is  rather  obvious  that  their 
density  is  given  by 

D,  =  D  p  (45) 

1  o  o 

After  the  first  generation  is  ignited,  each  unignited 
structure  faces  v  D^/Dq  burning  structures  instead  of  v  and 
has  the  ignition  probability 

Pl  =  1  '  (1  -  PSR)Vdl  <46) 

where  d^  is  D-^/D^  so  t^ie  structure  density  of  the  second 

generation  is 


D  p, 
orl 


(47) 
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Putting  (1  -  PrD)v  =  P,  the 

OK 

generations  can  be  written  as 


relative  densities  of  successive 


(48) 


(  •  I 


d 

n 


It  should  be  noted  that  this  model  by  itself  does  not  say 
anything  about  the  rate  of  spread  or  the  spread  distance.  Pos¬ 
sibly  the  density  of  burned  structures  after  the  end  of  the  fire 
decreases  gradually  at  increasing  distance  from  the  initial  fire, 
with  the  result  that  the  burnout  area  has  no  sharp  boundary. 
However,  if  we  define  a  mean  spread  area  as  one  filled  with  the 
total  number  of  burned  structures,  then  the  mean  spread  distance 

oc 

is  given  by  2  d  units  of  distance, 
o  n 

To  investigate  the  convergence  of  the  series  1  +  PQ  +  p^  + 
p^  4*  ...»  we  use  the  relations 

P  =  1  -  Pq  *  1  -  (1  -  e'a)  =  e"a  (49) 


and  express  d  as 
n 


-ad 


d  =  1  -  e 
n 


n-1 


which  is  of  the  same  form  as  the  relation 

"aBn-l^ 

B  (»)  =  1  -  e  1 
n 


(50) 


(51) 
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in  the  stochastic  linear  model.  It  has  been  shown  that  the  series 
S  B  diverges  for  a  s  1.  The  same  conclusion  applies  to  the  sum- 


(52) 

(53) 

as  condition  for  infinite  spread. 

One  feature  of  the  second  model  is  that  it  applies  equally 
well  to  two-dimensional  and  linear  fuel  areas.  Special  cases, 
such  as  that  of  a  single  structure  constituting  the  original  fire, 
can  be  readily  handled. 

Spread  Computation 

As  in  Model  II,  the  fuel  area  is  divided  into  equal-size 
cells.  At  any  time  each  cell  is  in  one  of  five  states:  suscept¬ 
ible  -(ignitable  but  unignited),  immune  (unigni table) ,  ignited 
(burning  but  still  incapable  of  extra-cell  spread),  burning 
(capable  of  spreading  to  adjacent  cells),  and  burned  out.  Usually 
a  susceptible  cell  is  ignited  on  one  edge  and  some  length  of  time 
is  required  for  the  fire  to  span  the  cell  dimension  and  reach 
adjacent  cells.  This  is  the  difference  between  the  ignited  state 
and  burning  state.  The  travel  time  across  the  cell,  also  called 
incubation  time,  is  determined  by  cell  size  and  rate  of  spread. 


n 


mat ion  2  d^.  For  a  z  1,  Eq.  (49)  gives 


P  £  e 


-1 


and  by  definition  of  P,  we  again  have 


-  Wv  * e'1 

PSRa  1  '  e'U/V) 


The  following  information  is  required  to  implement  this 
model : 

•  The  type  of  fuel  in  each  cell 

•  A  spread  vs  no-spread  table  giving  the  range  of  weather 
conditions  for  eac^i  type  of  rural  fuel  and  fuel  con¬ 
ditions  for  each  type  of  urban  fuel  under  which  the  fuel 
is  either  susceptible  (fire  will  spread  through  the  fuel 
indefinitely) ,  or  immune  (fire  will  not  spread  through 
the  fuel) 

•  A  rate  of  spread  table,  giving  the  rate  of  fire  spread 
for  each  fuel  type  as  a  function  of  the  pertinent  weather 
and  topography  conditions 

•  A  burning- time  table,  giving  the  burning  time  for  each 
fuel  type  as  a  function  of  pertinent  weather  variables 

•  For  rural  fuel,  weather  conditions  for  seven  days  prior 
to  the  assumed  starting  time  and  a  forecast  of  future 
weather  conditions  (historical  data  may  be  used  if 
adequately  forecasted  information  is  not  available) 

•  The  initial  burning  conditions  (the  cells  that  are 
initially  ignited,  burning,  or  burned  out) 


The  first  step  in  applying  the  model  is  to  establish  the 
initial  state  of  each  cell,  which  can  be  accomplished  from  the 
fuel  and  weather  data,  the  spread  vs  no-spread  table,  and  the 
specified  initial  burning  conditions. 

The  state  of  each  cell  is  then  redetermined  at  the  end  of 
each  successive  small  ti^e  increment  At  until  a  weather  change 
occurs  using  the  following  rules: 

•  A  cell  in  the  immune  state  will  remain  in  the  immune 
state 

•  A  cell  in  the  susceptible  state  will  instantly  change 

to  the  ignited  state  if  the  cell  is  immediately  adjacent 
to  a  cell  in  the  burning  state 
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•  A  cell  in  the  ignited  state  will  change  to  the  burning 
stace  when  the  time  after  ignition  is  equal  to  or  greater 
than  the  incubation  time 

•  A  cell  in  the  burning  state  will  change  to  the  burned- out 
state  when  the  time  after  start  of  the  burning  state  is 
equal  to  or  greater  than  the  burning  time  of  the  fuel 

•  A  cell  in  the  burned-out  state  will  remain  in  the  burned- 
out  state 

Figure  7  illustrates  the  changes  in  a  fire  that  can  occur  in 
time  interval  At. 

When  weather  conditions  change  significantly,  all  cells  of 
rural  fuel  must  be  reexamined  for  possible  change  of  state;  also, 
for  both  rural  and  urban  fuel,  rates  of  spread  and  incubation 
times  must  be  adjusted. 
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A.  Fire  front  and  state  of  each  cell  at  time  t 


B  =  burning  S  «  susceptible 

0  —  burned  out  M  **  immune 

I  *  ignited 


B.  Fire  front  and  state  of  each  cell  at  time  (t  +  At) 


r 

L 


Fig.  7.  Two-Dimensional  Deterministic  Model 
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Section  7 
PARAMETERS 

IMPLICIT  ROLE  OF  BASIC  VARIABLES  IN  MATHEMATICAL  MODELS 

All  characteristics  of  a  fire  are  determined  by  three  groups 
of  basic  variables,  fuel,  weather,  and  topography.  However,  these 
variables  do  not  appear  explicitly  in  the  mathematical  expressions. 
Instead  they  are  regarded  as  joint  determinants  of  two  or  three 
fundamental  properties  which  enter  the  equations  as  parameters 
and  are  sufficient  to  provide  complete  representation  of  fire 
spread  phenomena.  The  derivation  of  these  parameters  from  the 
basic  variables  constitutes  3  major  step  in  building  semiempirical 
models.  In  the  following  paragraphs,  we  shall  discuss  the  sig¬ 
nificance  of  these  parameters  and  the  method  of  determining 
their  values  from  empirical  data. 

COMPARISON  OF  PARAMETERS 

We  recall  that  in  the  fire  front  model  the  entire  fire  is 
regarded  as  a  random  walker  which  moves  along  a  row  of  cells  and 
which,  during  the  infinitesimal  time  dt,  has  the  chance  ji^dt  of 
dying,  the  chance  Xdt  of  jumping  forward,  and  the  chance 
1  ”  (m.£  +  X)dt  of  staying  where  it  is.  The  probability  rate  X 
represents  the  dynamic  aspect  of  the  fire,  its  tendency  to  move 
forward;  hence  it  is  called  the  spread  parameter.  The  probability 
rate  on  the  other  hand  is  a  measure  of  the  tendency  of  the  fire 
to  go  out  and  is  termed  decay  parameter.  It  is  intuitively  clear 
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that  each  of  these  parameters  is  determined  to  various  extents 
by  the  type,  density,  and  moisture  content  of  the  fuel,  the  veloc¬ 
ity  of  the  wind,  and  the  topographic  features  of  the  area. 

Superficially,  the  parameters  A  and  p  in  the  Models  II  and 

s 

IIA  are  similar  to  X  and  p^  in  Models  I,  IA,  and  IB.  Fundamental 

differences  exist  however  between  these  two  pairs  of  parameters. 

These  differences  arise  from  the  fact  that  II  and  IIA  describe 

the  fire  in  each  cell  separately,  whereas  in  I,  IA,  and  IB  the 

whole,  fire  front  is  an  entity  which  can  be  in  one  cell  at  a  time. 

Although  A  and  X  are  numerically  equivalent,  A  refers  to  the 

ignition  of  cell  j  by  cell  i  as  independent  from  what  happens  to 

the  fire  in  cell  i,  in  contrast  with  X,  which  is  the  chance  that 

the  fire  front  quits  cell  i  to  go  to  cell  j .  Similarly  p  refers 

s 

to  the  extinction  probability  of  the  fire  in  one  cell  and  p^.  to 
the  probability  of  permanent  extinction  of  the  whole  fire  front. 

It  can  be  expected  therefore  that,  in  general,  the  two  parameters 
have  different  numerical  values . 

The  deterministic  model  IIB  uses  the  three  parameters  k,  t, 

and  W  .  The  first  one  is  a  measure  of  the  heat,  flames,  and 
o 

firebrands  evolved  per  unit  area  of  burning  fuel.  Consequently 

it  is  the  deterministic  counterpart  of  A  and  X.  The  parameter  t, 

the  burning  time  of  the  fuel,  corresponds  to  l/.p  .  Again  it  is 

s 

intuitively  clear  that  both  k  and  t  are  functions  of  the  basic 

fire  spread  variables.  W  is  the  result  of  a  mathematical  sim- 

o 

plification.  The  igniting  flux  per  unit  burning  fuel  naturally 
decreases  as  the  distance  from  the  leading  edge  increases.  The 
assumption  is  that  the  igniting  flux  can  be  made  uniform  within 
a  distance  from  0  to  W  . 
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The  model  IIC  is  a  practical  method  for  following  the  prog¬ 
ress  of  the  fire  by  direct  application  of  data  on  basi.c  variables. 
The  mathematical  framework  is  that  of  models  IX  and  IIA,  includ¬ 
ing  all  parameters  used  in  these  models. 

DETERMINATION  OF  PARAMETERS 

Each  parameter  combines  the  effects  of  several  basic  varia¬ 
bles.  One  characteristic  of  semiempirical  models  is  that  param¬ 
eter  values  are  derived  from  empirical  sources,  such  as  observa¬ 
tions  or  measurements  of  actual  fires,  and  experimentation  on 
reduced-scale  fires.  In  the  following  paragraphs  we  shall  dis¬ 
cuss  how  existing  records  of  past  fires  and  certain  subjective 
estimates  by  Forest  Service  personnel  are  used  to  derive  some  of 
the  parameters.  It  will  appear  in  this  discussion  that  many 
important  quantities  in  fire  spread  have  not  been  made  the 
object  of  quantitative  observation  in  the  past.  For  parameters 
which  cannot  be  evaluated  immediately  because  of  the  lack  of 
pertinent  data,  we  shall  briefly  indicate  procedures  for  making 
useful  observations  on  future  fires ,  thus  accelerating  the 
build-up  of  our  empirical  knowledge  on  mass  fire  behavior. 

Spread  Parameters  X  and  A 


It  would  be  very  difficult  to  design  experiments  or  to  find 
observed  data  which  can  directly  yield  the  probability  rates  A 
and  X  as  such.  We  have  shown  analytically,  however,  that  X  and 
A  are  numerically  equal  to  the  mean  velocity  of  the  fire,  given 
that  it  is  still  alive.  If  we  observe  a  fire  front  during  a 
sufficiently  large  number  n  of  time  increments  At  before  it  dies. 
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we  would  see  that,  on  thfe  average,  it  makes  a  jump  during  each  of 
nXdt  time  increment,  that  during  the  other  n(l  -  Xdt)  time  incre¬ 
ments  it  stays  where  it  is,  and  that  the  jumps  and  stays  are 
randomly  distributed.  Thus,  the  rate  of  motion  from  cell  i  to 
cell  j  is  a  random  variable,  but  if  the  number  of  cells  traveled 
after  a  large  number  of  time  increments  is  recorded,  then  the 
mean  velocity  and  X  can  be  computed. 

Decay  Parameters  tig  and  ps 

It  is  shown  in  Appendix  that  according  to  the  mechanism 
used  in  the  fire  front  model,  l/p^  is  the  mean  lifetime  of  the 
fire.  Thus  if  the  mean  lifetime  of  the  fire  under  each  set  of 
basic  variables  is  known,  p  can  be  easily  derived.  Unfortunately 
the  only  direct  method  for  determining  the  mean  lifetime,  except 
for  a  very  special  case  mentioned  in  Appendix  A,  is  to  measure 
the  actual  lifetimes  of  a  large  number  of  fires  which  bum  under 
known  weather,  fuel,  and  topography  conditions,  and  unhampered 
by  fire  control  measures.  There  are  not  many  real  fires  satis¬ 
fying  these  conditions,  and  direct  experimentation  is  quite 
impractical.  An  indirect  method  involves  calculating  the  mean 
lifetime  from  some  other  fire  spread  model.  The  reciprocal  of 
the  mean  lifetime  thus  determined  is  taken  as  p^.  In  some  par¬ 
ticular  cases  the  mean  lifetime  can  be  calculated  from  fuel  dis¬ 
tribution,  as  discussed  at  the  end  of  Appendix  A. 

Since  1/jjl  is  the  mean  burning  time  of  cells  in  a  given 
s 

area,  the  determination  of  p  requires  the  burning  time  of  each 

s 

type  of  cell  and  their  distribution  in  a  given  area.  If  the 
weather  is  reasonably  constant  and  all  cells  are  identical  with 
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respect  to  fuel  and  topography,  the  lifetime  is  the  same  for  all 
cells,  and  p.  is  a  deterministic  parameter.  The  fuel,  in  any 
reasonably  large  area,  however,  is  likely  to  have  a  nonuniform 
distribution.  Some  cells  are  empty,  the  rest  contain  fuel  at 
various  densities  (fraction  of  total  area  actually  occupied  by 
fuel)  and  various  degrees  of  buildup  (thickness  of  fuel  layer) . 
Furthermore  any  area  is  likely  to  contain  more  than  one  type  of 
fuel.  Weather  and  topography  are  probably  not  critical  factors 
for  burning  times.  Their  small  effects,  however,  fluctuate  both 
in  time  and  in  space.  These  considerations  lead  to  the  view¬ 
point  that,  for  large  areas,  burning  time  is  predominantly 
stochastic.  In  the  fuel  state  models,  the  treatment  of  burning 
time  as  a  stochastic  variable  is  based  on  the  assumption  that 
the  cumulative  distribution  of  actual  burning  times  (or  of  cells 
according  to  their  burning  times)  is  exponential,  i.e.,  the 

probability  of  a  cell  or  the  percentage  of  cells  having  a  bum- 

-p.t 

ing  time  equal  to  or  greater  than  t  is  e  s  .  It  should  be 
pointed  out  that  this  assumption  is  not  essential  to  the  model. 
While  the  exponential  distribution  is  mathematically  convenient, 
other  types  of  distribution,  either  postulated  or  experimentally 
determined,  can  be  used  if  the  computation  is  carried  out  by 
machine.  It  has  been  shown  that  the  type  of  distribution  used 
has  no  strong  effect  on  the  final  results  as  long  as  the  mean 
burning  time  remains  essentially  the  same. 

Determination  of  the  mean  burning  time  must  be  carried  out 
in  two  steps:  first,  measuring  the  deterministic  burning  time 
of  fuel;  second,  evaluating  the  fuel  distribution  among  the  cells 
of  an  extended  region. 
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In  the  first  step,  the  fuel  may  be  divided  into  several  types 
such  as  brush,  grass,  conifer,  timber,  hard-wood  timber,  commer¬ 
cial  structures,  industrial  structures,  and  residential  struc¬ 
tures.  Each  fuel  type  may  be  broken  down  into  several  sub-types 
corresponding  to  structural  characteristics  and  degree  of  buildup 
(as  defined  earlier) .  The  correct  classification  can  be  deter¬ 
mined  only  after  study  of  the  relative  importance  of  fuel  type, 
structural  characteristics,  and  degree  of  buildup.  The  burning 
time  of  each  fuel  subtype  is  measured,  and  the  effects  of  weather 
and  topography  upon  this  time  are  investigated.  This  measurement 
should  be  carried  out  in  a  manner  consistent  with  the  mechanism 
employed  in  the  model.  Thus  the  burning  time  should  be  the  in¬ 
terval  during  which  the  cell  is  capable  of  igniting  an  adjacent 
cell.  The  mode  of  cell  ignition  should  approximately  duplicate 
cell  ignition  that  takes  place  during  actual  fires.  The  cell 
width  has  been  defined  as  the  distance  from  the  ignition  source 
beyond  which  the  probability  of  normal  ignition  (by  heat  and 
short-range  embers,  excluding  long-range  spotting)  is  negligibly 
small.  For  urban  fuel  a  cell  is  a  structure  (square  cell)  or  a 
row  of.  structures  (strip  cell),  and  the  exact  cell  width  is  de¬ 
termined  by  the  fuel  density  and  the  average  dimension  of  struc¬ 
ture.  Cell  dimension  in  the  essentially  continuous  rural  fuel 
is  selected  in  accordance  with  maximum  range  of  normal  ignition 
which  is  a  characteristic  of  fuel  type. 

In  the  second  step,  it  is  seen  that  whereas  burning  times 
can  be  measured  once  for  all,  the  spatial  distribution  of  fuel  is 
a  characteristic  of  the  area  under  consideration  and,  strictly 
speaking,  must  be  determined  in  each  case  by  counting  the  numbers 
nl,  n2’  n3J  ***  ceH-s  with  burning  times  T^,  ,  T~  and 
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calculating  the  mean  burning  time  t  from  the  formula 


t 


Sn.T. 
1  i. 

2n. 


(54) 


For  very  large  areas,  cell  counting  would  require  excessive  work. 
A  more  practical  method  is  as  follows.  A  number  of  typical  areas 
are  selected  from  maps  of  cities,  suburbs,  and  wildland.  The 
fuel  distribution  is  established  in  detail,  and  the  mean  burning 
time  is  calculated  for  each  area  type.  Mean  burning  times  are 
tabulated  for  future  use  in  each  model  application.  To  obtain  a 
grand  mean  burning  time,  it  is  required  only  to  estimate  the  dis¬ 
tribution  of  area  types  in  the  total  area  under  consideration. 


Deterministic  Parameters  t,  k,  and  Wn 

In  model  IIB  the  parameter  t  is  the  deterministic  burning 
time  of  the  uniform  fuel  or  uniform  mixture  of  fuel,  and  it  can 
be  measured  as  described  earlier  in  connection  with  the  burning 
time  of  various  fuel  types.  The  other  two  parameters,  k  and  WQ, 
can  be  obtained  if  certain  quantities  are  measured  on  actual 
fires . 

Rate  of  spread  data,  already  accumulated  in  substantial 
quantity,  correspond  to  an  essentially  stationary  state,  i.e., 
a  state  in  which  the  width  of  the  fire  front  and  the  velocity 
remain  constant.  According  to  the  simplified  formulation  of 
model  IIB,  these  quantities  are  given  by 

W(t)  =  kW  r  and  ■—  * 
x  7  o  dt  o 


(55) 
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Thus  for  each  area  type  and  each  combination  of  variables  the 
product  kWQ  is  equal  to  the  constant  rate  of  spread.  If  the 
width  kWQT  is  also  observed,  the  two  equations  can  be  solved 
for  t.  For  cases  in  which  t  has  been  measured  independently, 
this  will  provide  a  good  check. 

To  obtain  Wq  and  k  separately,  the  rate  of  spread  and  front 
width  at  the  initial  stage  must  also  be  observed.  The  data  will 
be  used  to  calculate  k  from  the  equation 

ff  =  kW(t)  (56) 

Observation  of  the  early  stage  of  fire  spread  may  be  ex¬ 
tended  to  fires  which  go  out  soon  after  start.  In  this  case 
the  initial  velocity  should  increase  with  the  initial  width  W(0) 
and  become  constant  when  W(0)  is  equal  to  Wq.  The  parameter  k 
is  obtained  from  the  instantaneous  value  of  velocity  and  front 
width  according  to  the  above  equation. 

It  is  possible  that  for  a  given  set  of  conditions,  k  is 
not  a  true  constant.  This  would  be  due  to  inaccuracy  of  the 
data  and  especially  to  the  assumption  that  velocity  is  propor¬ 
tional  to  the  front  width.  Unless  the  relation  between  these 
quantities  is  established  experimentally  or  theoretically,  the 
value  of  k  must  be  averaged  over  a  sufficiently  large  number 


of  observations. 
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Section  8 

APPLICATION  OF  MODELS 

Provided  all  parameters  can  be  determined,  we  now  have  a 
group  of  semiampirical  models  to  represent  spread  character¬ 
istics  of  mass  fires  and  obtain  information  required  for  damage 
prediction.  Some  of  these  models  are  stochastic,  some  determin¬ 
istic,  and  each  type  includes  both  linear  and  two-dimensional 
varieties.  The  question  now  arises  as  to  which  of  these  models 
is  most  suitable  for  a  given  type  of  fire. 

DETERMINISTIC  MODELS  VERSUS  STOCHASTIC  MODELS 

We  said  earlier  that  weather  variables  can  be  treated  de¬ 
terministically  because  of  their  relatively  slow  variations  and 
because  of  the  possibility  of  specifying  average  weather  con¬ 
ditions  which  remain  applicable  for  periods  of  6  or  12  hours. 
Furthermore,  investigation  of  fire  spread  data  in  rural  fuel 
has  shown  that  the  opposing  effects  of  upslopes  and  downslopes 
tend  to  cancel  each  other  to  such  an  extent  that  the  net  result 
over  large  distances  cannot  be  detected  in  the  best  data  now 
available.  Thus  the  choice  between  stochastic  models  and  de¬ 
terministic  models  is  essentially  determined  by  fuel  variables, 
i.e.,  type,  density,  degree  of  buildup,  and  distribution.  It 
seems,  worthwhile  to  examine  each  of  these  variables  for  rural 
fuel  and  urban  fuel  separately  at  various  dimension  scales. 

The  main  types  of  rural  fuel  are  brush,  conifer  timber,  hard 
wood  timber,  and  grass.  There  is  evidence  that  rate  of  spread  in 
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each  of  these  types  is  approximately  the  same.  Their  burning 
times,  however,  vary  considerably.  The  fuel  density  or  fraction 
of  total  area  occupied  by  fuel  is  normally  very  high.  Varia¬ 
bility  in  buildup  or  vertical  dimension  of  fuel  depends  on  the 
number  of  different  fuel  types  in  the  area  under  consideration 
and,  therefore,  on  the  size  of  the  area.  Fuel  distribution,  a 
function  of  density  and  variability  in  buildup,  is  generally 
uniform  within  relatively  large  expanses.  A  rural  area  tends 
to  be  continuously  occupied  by  the  fuel  except  for  some  patches 
of  empty  or  almost  empty  area,  e.g.,  roads,  creeks,  or  rivers. 

Urban  fuel  is  essentially  discrete.  It  is  concentrated  in 
separate  units  whose  size,  structural  characteristics,  and  spac- 
ings  are  highly  variable.  The  density  in  general  increases  from 
residential  areas  to  commercial  areas  and,  except  for  certain 
residential  tracts,  is  highly  variable  locally.  Buildup  is 
naturally  uniform,  at  least  within  a  given  zone  although  it  may 
vary  over  a  wide  range  in  commercial  areas. 

It  thus  appears  that  rural  fuel  lends  itself  easily  to 
deterministic  treatment,  at  least  within  areas  of  a  few  square 
miles.  Moderately  large  areas  can  also  be  treated  in  this 
manner,  provided  input  data  can  be  collected  for  each  cection 
with  different  fuel  characteristics.  For  very  large  areas,  a 
stochastic  approach  might  spare  some  of  the  efforts  required 
for  data  collection  and  model  computation. 


For  urban  fuel,  the  variability  of  sizes,  structural  charac¬ 
teristics,  and  spacings  —  all  of  which  are  important  factors  in 
fire  spread  —  calls  for  stochastic  models  except  when  the  fire 


spread  investigation  is  restricted  to  a  small  area. 


t 
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LINEAR  MODELS  VERSUS  TWO-DIMENSIONAL  MODELS 


Obviously  a  two-dimensional  model  can  be  applied  to  any  fire 
spread  problem  regardless  of  the  spread  pattern,  and  situations 
frequently  arise  in  which  such  a  model  is  required.  In  some  cases, 
however,  a  linear  model  is  sufficient  and  due  to  the  simple  compu¬ 
tation  it  requires,  can  produce  substantial  savings  in  time  and 
effort . 


A  linear  model  describes  the  progress  of  a  fire  along  the 
main  spread  direction  (on  the  assumption  that  the  spread  be¬ 
havior  is  approximately  the  same  along  straight  lines  perpen¬ 
dicular  to  the  direction  of  spread) .  If  space  is  divided  into 
equal  cells,  the  forward  dimension  of  each  cell  is  set  by  the 
requirement  that  fire  cannot  jump  over  one  or  more  cells  by  the 
normal  mode  of  ignition.  The  lateral  dimension  can  take  any 
value,  provided  fuel  characteristics  (deterministically  or 
stochastically  specified)  do  not  vary  significantly  along 
this  dimension.  The  assumption  of  uniform  spread  throughout 
the  lateral  dimension  of  the  cell  is  applicable  in  two  cases: 

(1)  The  wind  velocity  is  such  that  the  lateral  component 
of  spread  velocity  is  small  compared  with  the  forward 
component.  One  such  case  is  pictured  in  Fig.  8A.  The 
size  of  the  initial  fire  is  immaterial  in  this  case. 

(2)  The  wind  velocity  is  low,  the  initial  fire  front  is 
straight  or  only  slightly  curved  and  has  a  consider¬ 
able  length  (Fig.  8B) .  Under  these  conditions,  a 
certain  portion  of  the  fire  front  is  approximately 
straight  and  remains  parallel  to  itself  in  its  for¬ 
ward  motion. 

A  case  requiring  special  attention  in  using  a  linear  model 
is  pictured  in  Fig.  8C.  The  model  is  applied  to  a  strip  of  fuel 
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Fig.  8.  Application  of  Linear  Model 


URS  641-6 


8-5 


area  extending  from  the  initial  fire  to  a  fuel  break  (open  space) . 
If  sufficiently  wide,  the  break  should  stop  the  fire  (going  along 
this  strip)  with  certainty.  However,  the  possibility  for  the 
fire  to  go  around  the  edge  of  the  fire  break  and  start  again  on 
the  other  side  must  always  be  examined. 

If  the  position  of  the  fuel  break  is  stochastic,  the  entire 
area  is  stochastically  uniform  and  can  be  treated  as  a  general 
case . 


When  the  situation  permits  the  use  of  a  linear  model,  the 
general  procedure  is  as  follows: 

(1)  From  input  data,  determine  whether  the  fire  will  stop 
at  a  finite  distance  (A/us  or  kT  smaller  than  1)  or 
spread  indefinitely  (A/us  or  kr  equal  to  or  greater 
than  1) . 

(2)  In  the  first  case  calculate  the  final  spread  distance 
using  appropriate  formulas  and  from  the  final  spread 
distance  and  average  rate  of  spread  calculate  the 
lifetime  of  the  fire. 

(3)  In  the  second  case  the  fire  will  spread  until  it  en¬ 
counters  a  spatial  barrier  or  a  weather  barrier.  A 
spatial  barrier  is  a  completely  empty  or  low-density 
area  whose  size  and  location  are  specified  in  the 
input  data.  A  weather  barrier  is  a  set  of  wind  and 
humidity  conditions  which  inhibits  fire  spread  in 
the  fuel  area  under  consideration.  In  view  of  pre¬ 
vious  discussions  on  fuel  characteristics,  urban 
fires  are  stopped  by  spatial  barriers  more  frequently 
than  by  weather  barriers,  and  the  reverse  is  true  with 
rural  fires. 

(4)  If  fire  spread  is  expected  to  be  multidirectional 
(low  wind  velocity) ,  determine  the  final  burnout 
contour  from  spread  distances  in  various  directions. 
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(5)  Position  of  the  fire  at  a  given  moment  can  be  found 
from  the  rate  of  spread  and  time  after  start. 

Thus  in  many  cases  linear  models  only  require  derivation  of 

parameter  values  from  input  data  and  numerical  computation  of 

certain  simple  quantities  such  as  A/p  or  kr.  Both  of  these 

s 

tasks  can  be  accomplished  with  or  without  the  help  of  a.  computer 
depending  on  the  particular  situation. 

Table  1  summarizes  the  characteristics  and  applicability 
of  various  models . 

MODEL  APPLICATION  AS  PERMITTED  BY  EXISTING  DATA 

As  appears  in  Table  1,  only  some  parameters  can  be  empiri¬ 
cally  established  from  existing  data.  Consequently,  some  models 
are  not  immediately  ready  for  actual  application.  Thus  for 
models  1,  IA,  and  IB,  the  parameter  would  have  to  be  derived 
from  the  mean  lifetime  of  the  fire,  a  quantity  which  cannot  be 
obtained  from  any  observed  data.  In  the  case  of  models  II  and 
IIA,  the  determination  of  p„  (reciprocal  of  mean  burning  time) 
would  require  information  on  burning  time  of  various  types  of 
fuels,  which  would  be  combined  with  spatial  distribution  of 
fuel  to  give  the  mean  burning  time.  Some  estimates  of  burn¬ 
ing  time  do  exist,  but  their  reliability  is  not  considered 
sufficient  for  this  purpose.  Burning  time  is  also  a  parameter 
in  model  IIB. 

Until  further  observed  data  become  available  for  evaluation 
of  these  missing  parameters,  the  problem  is  how  to  perform  fire 
spread  predictions  on  the  basis  of  existing  data. 
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Table  1 

SYNOPSIS  OF  MODELS 


Model' 


Characteristics 


>ropriate  Fuel  Area 


Fire  front  mechanism.  Stochas¬ 
tic.  One  dimensional.  Param¬ 
eters  are  X  and  the  latter 
not  available ;  required  data 
difficult  to  obtain. 


Preferably  urban 


As  above  with  more  detail  input  Preferably  small  or 
information  on  fuel  character-  moderate  urban  areas 


lSClCS . 


As  above  with  provision  for  long  As  above 
range  spotting  possibility. 


Fuel  state  mechanism.  Stochas¬ 
tic.  Two  dimensional.  Param¬ 
eters  are  A  and  p.s .  Data  for 
evaluating  ( iig  not  available. 

As  above  but  one  dimensional. 


General 


General,  preferably 
urban 


Deterministic,  one  dimensional. 
Existing  data  for  evaluating 
parameters  t ,  k,  W0  are  in¬ 
complete  . 


General,  preferably 


rural 


Deterministic.  Two  dimensional. 
Application  uses  existing  data 
supplemented  by  some  subjective 
estimates . 


Rural  and  urban 


See  Fig.  3. 


8-8 


URS  641-6 


The  solution  is  provided  by  model  IIC.  It  is  assumed  that 
for  a  given  area  and  a  given  type  of  weather,  fire  either  spreads 
indefinitely  or  does  not  spread  at  all,  neglecting  the  rather  few 
cases  in  which  fire  spreads  to  a  significant  distance,  then  stop. 
For  rural  fuel,  personnel  of  the  Forest  Service  have  estimated 
weather  conditions  for  no-spread.  Limited  data  from  past  urban 
fires  have  provided  approximate  information  on  the  ignition  prob¬ 
ability  as  a  function  of  spacing,  from  which  fuel  conditions  for 
no-spread  have  been  derived.  If  fire  is  predicted  to  spread 
through  an  area,  its  position  at  any  time  can  be  determined  on 
the  basis  of  rate  of  spread  data. 


Part  II 


Work  described  in  Part  I  concerns  the  development  of  semi- 
empirical  models  fc'*  predicting  the  gross  spread  behavior  of 
mass  fires,  the  investigation  of  model  parameters,  and  methods 
for  obtaining  their  numerical  values  from  observed  data. 

Part  II  will  present  the  results  of  studies  on  the  required 
input  data  (rate  of  spread,  weather  conditions  for  spread  and  for 
extinction  in  rural  areas,  fuel  conditions  for  spread  in  urban 
areas,  and  specification  of  weather  characteristics)  and  on  the 
vulnerability  of  resources.  Finally,  some  preliminary  considera¬ 
tions  are  made  regarding  application  of  fire  spread  models  on 
the  local  scale. 
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Section  9 

ANALYSIS  OF  RATE -OF- SPREAD  DATA  IN  RURAL  FUEL 


RAW  DATA 

The  analysis  consists  in  examining  how  rate  of  spread  depends 
on  fuel  type  and  density,  wind  velocity  and  direction,  relative 
humidity,  slope  and  slope  distribution.  The  use  of  relative 
humidity  instead  of  moisture  content  is  justified  by  the  com¬ 
monly  observed  fact  that  in  natural  rural  fuel,  fire  spreads 
through  the  finer  components,  whose  moisture  content  follows 
the  fluctuations  of  relative  humidity  rather  closely. 

A' sample  of  raw  data,  as  supplied  by  the  U.S.  Forest  Service, 
is  shown  in  Table  2.  Columns  1  and  2  of  the  table  identify  the 
fires  and  the  section  of  the  fires  where  the  measurement  was 
made.  Column  3  indicates  the  time  of  start  and  Column  4  the 
time  during  which  fire  spread  was  measured.  Weather  variables 
are  stated  in  Columns  5  to  9,  the  last  three  of  which  refer  to 
relative  humidity,  stick  moisture  content,  and  burning  index, 
respectively.  Fuel  type  is  shown  in  Column  10.  Topography 
variables  are  given  in  Columns  11  to  16.  Rate  of  spread  and 
direction  in  which  spread  was  measured  are  indicated  in  Columns 
17  and  18.  Column  19  indicates  the  manner  of  spreading  in  the 
area  where  the  rate  of  spread  was  measured.  H  is  a  head  fire, 


*  In  addition  to  profile  sketches  in  Column  16,  slop,-  angles  for 
each  section  of  a  profile  were  provided  separately  but  not  in¬ 
cluded  in  this  report. 
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Table  2 

SAMPLE  FIRE  SPREAD  DATA 
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R  a  rear  or  backing  fire;)  F  a  flank  fire,  and  0  a  circular  fire. 
Detailed  explanations  of  this  table  has  been  given  in  Ref.  1. 

ANALYSIS 

The  first  problem  examined  was  the  effect  of  measured  time 
(of  spread)  on  rate  of  spread.  Preliminary  analysis  had  indicated 
that  short  measured  times  are  frequently  related  to  high  rates  of 
spread.  In  fact,  plots  of  rate  of  spread  against  measured  time 
for  various  ranges  of  wind  velocity,  wind  angle,  and  relative 
humidity  showed  that  this  is  true  for  times  shorter  than  8  or 
9  hr.  In  this  range,  the  apparent  dependence  of  rate  of  spread 
on  measured  time  is  possibly  due  to  a  tendency  to  take  data  on 
fire  location  more  frequently  for  fast  fires  than  for  slow  fires. 
To  minimize  possible  errors  due  to  this  effect,  data  with  meas¬ 
ured  time  less  than  9  hr  in  duration  were  rejected. 

The  topography  problem  was  taken  up  next.  In  general,  up- 
slopes  are  expected  to  increase  the  rate  of  spread  because  they 
favor  flame  contact  and  convection  heating.  In  contrast,  down- 
slopes  are  expected  to  slow  down  the  spread.  An  exception  is 
when  the  downs lope  is  steep  enough  so  that  burning  fragments 
tend  to  roll  down  and  spread  fire  ahead  of  the  fire  front. 

However,  in  some  cases,  the  fire  pockets  created  by  rolling 
fragments  develop  slowly  and  are  caught  up  by  the  main  fire 
front,  so  that  the  ultimate  effect  of  these  fragments  might 
not  be  noticeable. 

To  test  the  effect  of  slopes,  pairs  of  rate  measurements  on 
single  slopes  were  selected  for  which  all  basic  variables  except 
slope  were  the  same  for  both  members  of  the  pair.  Thus 
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any  difference  in  rate  of  spread  could  be  attributed  to  differ¬ 
ence  in  slope  only  and  not  to  any  other  measured  variable.  Each 
pair  was  marked  mR"  if  the  rate  of  spread  was  greater  for  the 
pair  member  with  either  greater  positive  slope  or  less  negative 
slope.  If  the  variation  in  rate  of  spread  was  opposite  to  that 
of  slope,  the  pair  was  marked  "W."  Of  a  total  of  52  pairs, 
selected  from  widely  different  ranges  of  variables,  24  were 
marked  R  and  28  were  marked  W.  This  result  is  taken  to  indi¬ 
cate  that,  within  the  accuracy  of  these  rate  of  spread  measure¬ 
ments,  the  effect  of  topography  is  insignificant.  Therefore, 
slope  and  slope  distribution  were  ignored  in  the  subsequent 
treatment  of  rate-of -spread  data. 

After  eliminating  topography,  four  variables  remain: 
fuel  type,  wind  velocity  (w) ,  wind  direction  (a)  (angle  between 
spread  direction  and  wind  direction) , and  relative  humidity  (r) . 


*  Another  reason  for  neglecting  topography  is  that  in  most 
practical  cases,  the  path  of  fire  spread  in  non-flat  areas 
during  time  periods  of  interest  will  not  be  on  a  single 
slope,  but  rather  over  a  number  of  up  and  down  slopes. 

Thus,  any  effect  of  one  type  of  slope  tends  to  be  compen¬ 
sated  by  that  of  another  type.  There  are  two  possible 
exceptions  to  this  compensating  tendency:  (1)  in  the  case 
of  a  series  of  small  steep  ridges,  the  fire  might  jump  from 
the  top  of  one  rising  slope  to  the  next  rising  slope  with¬ 
out  burning  the  falling  slope;  (2)  the  fire  might  remain 
stationary  for  some  time  at  the  top  of  a  ridge  with  the 
result  that  the  overall  rate  of  spread  for  a  combination 
of  upslopes  and  downslopes  is  lower  than  the  rate  of  spread 
for  a  flat  area.  However,  such  special  cases  are  beyond 
the  scope  of  our  considerations  at  present. 
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The  use  of  relative  humidity  as  a  variable  rather  than  fuel 
moisture  content  is  justified  by  the  fact  that  fire  tends  to 
spread  through  the  fine  component  of  the  fuel,  whose  moisture 
content  rather  closely  follows  the  variations  in  relative 
humidity. 

The  method  used  to  determine  the  effect  of  these  variables 
consisted  in  obtaining,  for  each  fuel  type,  a  plot  of  the  rate 
of  spread  R  versus  one  variable  for  various  ranges  of  the  other 
two  variables.  More  specifically, 

(1)  R  vs  a  for  all  possible  combinations  of  four  r-ranges 
(0  —  15,  15~30,  30  —  45,  45—  60%)  and  six  w-ranges 
(0-5,  5-  10,  10-  15,  15-20,  20~  25,  25-30  mph) 

(2)  R  vs  r  for  all  possible  combinations  of  six  w-deg 
ranges  and  three  c.-ranges  (0  —45,  45  —135,  135—  180  deg) 

(3)  R  vs  w  for  all  possible  combinations  of  four  r-ranges 
and  three  ? -ranges 

Brush  was  considered  first  because  data  for  this  fuel  type 
are  more  abundant  than  for  any  other.  However,  there  are  sev¬ 
eral  subtypes  and  several  mixtures  of  brush:  brush,  heavy  brush, 
medium  brush,  light  brush,  brush-chamise-chaparal ,  chamise, 
chamise-chaparal .  To  discover  possible  differences  between  these 
subtypes,  the  R  vs  n  plots  were  made  with  each  subtype  repre¬ 
sented  by  points  of  different  colors  or  different  shapes.  These 
plots  showed  at  a  glance  that  the  point  shapes  and  colors  would 
be  well  blended,  if  there  was  a  sufficient  number  of  points  on 
the  plot.  Accordingly,  all  brush  subtypes  were  considered  the 
same  with  regard  to  rate  of  spread. 

Similar  plots  were  also  made  to  show  successively  that 
brush  and  conifer  timber,  then  brush,  conifer  timber,  and  mixtures 
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of  both  do  not  differ  in  rate  of  spread  to  the  extent  that  is 
detectable  in  these  measurements. 

For  grass  and  hardwoods,  the  points  are  much  too  scarce  to 
yield  any  definite  result.  However,  there  is  some  evidence  that 
extrapolation  of  rate  of  spread  values  for  other  types  to  grass 
and  hardwood  might  not  cause  serious  error.  We  shall  assume  that 
this  is  the  case  until  more  data  are  available  for  grass  and 
hardwood. 

DISCUSSION 

There  resulted  three  series  of  plots  in  which  a,  r,  and  w 
were  taken  successively  as  the  running  variable  for  various 
combinations  of  ranges  of  the  other  two  variables.  One  w-plot 
is  shown  in  Fig.  9.  Each  plot  was  examined  in  an  effort  to 
derive  relations  between  R  and  the  running  variable.  It  was 
immediately  noticed  that  the  amount  of  scattering  is  consider¬ 
able.  Although  a  high  percentage  or  R- values  were  between 
zero  and  0.10  mph,  the  number  of  values  above  0.10  was  not 
negligible,  and  some  values  were  as  high  as  1.5  mph.  Thus, 
before  attempting  to  derive  relations  between  R  and  the  varia¬ 
bles,  some  consideration  of  the  causes  of  this  scattering 
seemed  necessary.. 

Low  precision  in  the  measurement  of  R,  a,  w,  r  could 
obviously  produce  some  scattering  but  did  not  seem  sufficient 
to  account  for  abnormally  high  values  such  as  the  ones  above 
a  few  tenths  of  a  mile  per  hour.  The  latter  were  probably 
due  to  long-range  spotting  (across  several  hundred  feet  or 
more).  Even  under  a  definite  set  of  conditions,  the  long-range 
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spotting  process  is  highly  random  with  respect  to  both  jump 
distance  and  jump  frequency.  Since  spread  by  spotting  is,  in 
general,  much  faster  than  ground  spread,  the  average  rate 
measured  over  a  certain  distance  can  be  expected  to  depend  on 
the  spotted  portion  of  that  distance.  Lack  of  information  on 
the  extent  of  spotting  made  it  impossible  to  correct  each-value 
for  contribution  from  this  mode  of  spread.  However,  the  fol¬ 
lowing  line  of  reasoning  was  used  in  an  attempt  to  derive 
reasonable  and  useful  results. 

Long-range  spotting  must  occur  predominantly  in  the  wind 
direction.  Thus,  if  c. -measurements  were  very  accurate,  and  the 
natural  fluctuations  were  not  too  great,  rates  of  spread  at  wind 
angles  greater  than  10— 20  deg  would  be  due  to  normal  spread 
alone.  However,  due  to  natural  fluctuations  and  measuring 
errors  in  wind  angles,  pure  normal  spread  could  be  found  only 
at  relatively  large  c.-values.  When  rates  of  spread  were  plotted 
against  a  for  all  r-  and  w-values  and  for  all  fuel  types,  it  was 
found  that 

(a)  The  number  of  points  higher  than  0.10  mph  in  the 
a- range  from  0  to  45  deg  was  greater  by  far  than 
at  higher  a-ranges. 

(b)  In  each  a-range  (0—10  deg,  10—20  deg,  etc.)  the 
density  of  points  (number  of  points  per  unit  area 
of  graph  paper)  decreased  in  going  from  low  to  high 
R-values,  the  rate  of  decrease  being  fastest  at 

R  ~  0.10. 

(c)  Some  abnormally  high  R  values  (>0.30  mph)  were 
scattered  throughcat  the  7,-scale,  even  at  its 
upper  end. 


URS  641-6 


9-9 


Observations  in  (a)  and  (b)  seemed  to  support  the  hypothesis 
that  rate  of  spread  values  above  0.10  mph  contained  significant 
contributions  from  spotting,  and  since  spotting  occurred  only  at 
small  wind  angles,  the  result  in  (c)  is  due  to  occasionally  large 
errors  in  a.  This  hypothesis  is  supported  by  results  shown  later 
in  this  report,  indicating  that  large  variations  in  wind  veloc¬ 
ity  and  relative  humidity  are  insufficient  to  raise  the  average 
rate  of  spread  above  0.10  mph,  even  when  the  average  rate  of 
spread  is  evaluated  by  including  all  values  above  0.10  mph. 

Thus,  in  determining  relations  between  rate  of  spread 
and  basic  variables,  all  points  above  0.10  mph  were  discarded. 
Analysis  of  these  points  would  help  our  understanding  of  the 
effect  of  basic  variables  on  long-range  spotting  were  it  not 
for  the  insufficient  number  of  such  points  and  the  lack  of 
other  important  information. 

Even  after  removal  of  all  points  above  0.10  mph,  each  plot 
of  R  against  some  variable  still  showed  considerable  scattering. 
However,  the  general  trend  qualitatively  agreed  with  the  usual 
views  on  fire  spread.  For  example,  in  a  plot  of  R  against  r 

(relative  humidity),  for  a  certain  range  of  a  and  w,  most  spread 
velocities  at  low  relative  humidity  are  higher  than  those  at  high 
relative  humidity,  and  an  upward  trend  is  noticeable  in  a  plot 
of  R  against  w  for  a  given  range  of  r  and  a  low  range  of  a 
(See  Fig.  9). 

Relations  between  R  and  the  variables  were  obtained  by  the 
usual  curve  fitting  process,  and  since  the  excessive  scattering 
made  it  unfeasible  to  determine  any  characteristic  other  than  a 
general  trend,  linear  relations  were  assumed  and  most  of  the 
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fitting  was  done  visually.  In  several  cases,  however,  least- 
square  fitting  was  also  carried  out,  which  showed  that  visual 
fitting  is  adequate  for  this  purpose.  For  the  case  of  Fig.  9 
the  resulting  least-square  equation  is 

R  -  0.01805  +  0.002232w 

The  most  probable  error  for  the  constant  0.01805  is  0.0027 
and  that  for  the  constant  0.002232  is  0.00025.  Regarding  the 
correlation  between  the  observed  rate  of  spread  Rq  and  w,  the 
standard  error  of  estimate  SR  ,  i.e.,  the  square  root  of  the 
mean  of  squared  deviations  of  R0  from  values  predicted  by  the 
least-squares  line  is  0.021  mph.  The  correlation  coefficient 
r  (defined  by  r  =  ^1  -  (Sro/cRq)2,  where  gRq  is  the  standard 
deviation  from  mean,  not  from  least-square  lines)  is  equal  to 
0.6.  Since  the  limits  of  r  are  1  for  complete  correlation  and 
0  for  complete  lack  of  correlation,  an  r-value  of  0.6  indicates 
that  the  assumption  of  linearity  is  not  too  unreasonable  in  view 
of  the  data  scattering. 

Each  of  the  three  series  of  plots  with  fitted  straight  lines 
gave  a  fitted  average  value  of  R  for  each  combination  of  a-ranges 
0-45,  45-135,  135-180  deg,  w-ranges  0-5,  5 -10,  10-15,  15-20 
20—25,  25  —30  mph,  and  r-ranges  0—15,  15“ 30,  30  — 45,  45  — 60  per 
cent.  Average  values  thus  obtained  are  shown  in  Table  3.  It  is 
seen  that  the  three  ways  of  plotting  lead  to  three  averages  which 
are  reasonably  close  to  each  other.  These  averages  are  averaged 
once  more  to  give  the  numbers  with  bars. 
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Notes:  In  each  box,  the  left  column  corresponds  to  a  =  0—45,  the  middle  to  a  =  45  —  135 

and  the  right  to  a  =  135—180  deg. 

Numbers  in  the  first  row  are  from  plots  of  R  vs  w,  those  in  the  second  row  from 
plots  of  R  vs  a,  those  in  the  third  row  from  plots  of  R  vs  w,  and  those  in  the 
last  row  are  averages  of  the  three  numbers  above  in  the  same  column. 
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Conclusions 


Examination  of  Table  3  leads  to  the  following  conclusions: 

(1)  In  each  range  of  relative  humidity  (r) ,  the  rate  of 
spread  increases  with  w  rather  markedly  for  a  between 
0  and  45  deg  and  less  markedly  for  a  between  45  and 
135  deg.  It  generally  decreases  at  increasing  w  for 
a  between  135  to  180  deg. 

(2)  At  low  wind  velocity  (w  from  0  to  10  mph) .  the  rate  of 
spread  decreases  with  increasing  relative  humidity,  the 
decrease  being  more  pronounced  in  the  0— 5  mph  range 
than  in  the  5  — 10  mph  range. 

(3)  At  higher  wind  velocity,  the  inhibiting  effect  of 
relative  humidity  becomes  very  small  and  should 
probably  be  ignored. 

(4)  For  w  between  0  and  5  mph,  the  rate  of  spread  does 
not  depend  on  wind  direction. 

Conclusion  (1)  is  in  accord  with  our  intuition.  As  the  wind 
becomes  stronger,  fire  spreads  faster  in  the  direction  of  the  wind, 
and  the  spread  pattern  becomes  more  and  more  asymmetric  (provided, 
of  course,  that  the  fuel  area  is  homogeneous).  Conclusion  2  is 
also  expected.  Conclusion  (3)  is  at  first  rather  surprising. 

Tt  seems  to  indicate  that,  at  moderate  and  high  wind  velocities, 
the  convection  heat  is  more  than  sufficient  to  drive  all  moisture 
content  off  the  fine  fuel  adjacent  to  the  fire  front  before  igni¬ 
tion  takes  place.  Conclusion  (4)  may  indicate  that  the  resolution 
of  w-measurements  is  no  better  than  5  mph. 
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Section  10 

ANALYSIS  OF  RATE -OF -SPREAD  DATA  IN  URBAN  FUEL 

A  sample  of  data  provided  by  the  U.S.  Forest  Service  is 
shown  in  Table  4  (taken  from  Ref.  1).  The  duration  of  meas¬ 
urement  can  be  found  from  the  times  of  start  and  finish  in 
column  3.  Column  4  gives  the  rate  of  spread  in  miles  per  hour. 
Wind  speeds  given  in  column  6  are  also  in  miles  per  hour. 

Column  10  indicates  the  fuel  type  for  light  wooden  struc¬ 
tures,  for  heavy  wooden  structures,  for  light  stone  or  concrete 
structures,  and  for  heavy  stone  or  concrete  structures,  indi¬ 
cated  by  the  numbers  1,  2,  3,  4,  respectively. 

The  percentage  of  total  area  occupied  by  fuel,  termed 
"built-upness,"  is  given  in  column  11.  We  prefer,  however,  to 
call  this  variable  "building  density"  or  "fuel  density,"  re¬ 
serving  the  term  "built-upness"  (degree  of  buildup)  for  the 
thickness  of  the  fuel  layer.  The  four  ranges  of  fuel  density 
(in  percent),  below  20,  20  —  29,  30  —  39,  40  and  above,  will  be 
denoted  by  I,  II,  III  and  IV,  respectively.  The  number  of  stories 
in  each  building  (a  measure  of  "built-upness"  is  given  in 
column  12. 

"V"  in  column  13  is  a  value  obtained  by  combining  fuel 
type  (T) ,  fuel  density  (B) ,  and  number  of  stories  (S)  in  a  cer¬ 
tain  manner.  This  value  has  not  been  used  in  our  analysis. 

The  direction  of  spread  with  reference  to  north  is  given 
in  most  cases  in  column  12  as  angle  to  wind.  The  angle  between 
the  spread  direction  and  wind  direction  is  found  from  the  spread 
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direction  in  column  15  and  the  wind  direction  in  column  7.  It 
was  understood  that  for  cases  in  which  the  spread  direction  is 
not  Indicated  in  column  15,  the  spread  direction  is  forward, 
i.e.,  with  the  wind. 

There  is  a  total  of  73  values  of  rate  of  spread,  of  which 
two  have  been  rejected  because  the  values  appear  to  refer  to 
spread  within  single  structures  and  another  because  of  unspeci¬ 
fied  type  of  spread  (ground  spread  or  spotting). 

METHODS  AND  RESULTS 

The  data  were  first  grouped  according  to  fuel  types  and 
fuel  densities.  The  numbers  of  data  points  in  each  group  are 
shown  in  Fig.  10,  with  a  breakdown  into  ground  spread  and  spot¬ 
ting.  There  is  only  one  data  point  for  each  type  of  spread  in 
group  IV-4.  Consequently,  this  group  was  later  incorporated 
into  group  IV- 3. 

Two  facts  were  immediately  noted.  First,  a  certain  relation 
exists  between  fuel  type  and  fuel  density.  Light  fuel  occurs 
predominantly  in  low-density  areas  and  heavy  fuel  in  high-density 
areas.  Thus,  care  must  be  taken  when  attributing  some  effect  to 
one  of  these  two  variables  to  the  exclusion  of  the  other.  Sec¬ 
ond,  ground-spread  fires  tend  to  occur  in  high-density  or  heavy 
fuel  areas  more  often  than  in  low-density  or  light  fuel  areas, 
whereas  the  reverse  is  true  for  spotting  fires.  The  number  of 
stories  naturally  tends  to  follow  the  fuel  type  and  may  well  be 
combined  with  the  latter  to  form  one  single  variable. 


Fuel  Type 
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Fig.  10.  Number  of  Data  Points  in  Each  Group 
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The  next  obvious  step  was  to  separate  ground-spread  data 
from  spotting  data  and  treat  them  separately,  because  these  two 
modes  of  spread  involve  different  mechanisms.  By  ground  spread, 
the  fire  reaches  an  adjacent  structure  through  a  relatively  short 
distance.  The  mechanism  is  spontaneous  ignition  or  pilot  ignition 
(ignition  of  preheated  fuel  by  embers).  Spotting  occurs  when  a 
structure  located  beyond  the  range  of  heat  and  embers  is  ignited 
by  flying  firebrands.  It  is  important  to  note  that  sufficiently 
long  spread  runs  are  likely  to  involve  both  types  :  spread. 

Thus,  for  most  runs  reported  as  ground  spread,  a  small  portion 
of  the  jumps  are  probably  made  by  spotting.  Similarly,  some 
ground  jumps  may  be  expected  to  occur  in  relatively  long-lasting 
fires . 

Ground  Spread 

We  shall  examine  ground  spread  data  first.  In  each  group, 
variables  which  can  have  some  influence  on  rate  of  spread  are 
wind  velocity,  wind  direction,  measured  time  of  spread,  and  fuel 
dryness.  There  are  only  9  cases  of  wet  conditions,  of  which  7 
are  in  group  IV-3,  and  2  in  group  III-2  (Fig.  11).  In  addition, 
there  exist  11  cases  with  unknown  dryness  (San  Francisco  fire, 
April  1906).  Two  of  these  are  in  group  III-3,  5  in  group  III-2, 

3  in  group  III-l,  and  1  in  group  II-2.  We  shall  later  find  a 
plausible  effect  of  fuel  dryness  in  group  IV-3  which  contains  8 
cases  of  dry  condition,  7  cases  of  wet  conditions,  and  no  case 
of  unknown  dryness.  However,  the  magnitude  of  this  effect  in 
group  IV-3  does  not  seem  to  interfere  with  an  independent  study 
of  the  other  three  variables.  From  this,  we  conclude  that  the 
effect  of  fuel  dryness  in  other  groups  also  must  be  of  little 
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Fig.  11.  Average  Values  for  Ground  Spread  (mph) 
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consequence  in  the  study  of  other  variables.  Furthermore,  most 
data  points  in  each  group  belong  to  the  same  type  of  dryness ; 
for  example,  3  out  of  3  data  points  in  group  III-l,  5  out  of  7 
in  group  III-2,  and  2  out  of  3  in  group  III-3  are  of  unknown 
dryness o  Group  II-2  has  only  2  points,  one  of  unknown  condition, 
the  other  of  dry  condition  «  However,  the  change  in  rate  of 
spread  between  these  2  data  points  can  be  easily  accounted  for 
by  other  effects. 

Consequently,  it  seems  justified  to  neglect  fuel  dryness 
at  this  stage  and  proceed  with  the  examination  of  wind  direction, 
wind  velocity,  and  measured  time. 

Effect  of  wind  direction  is  a  simple  problem  since  of  the 
41  cases  only  5  are  backward  spread  and  all  others  are  forward 
spread.  If  we  assume  that  for  a  given  wind  velocity,  the  di¬ 
rectional  effect  of  wind  is  the  same  for  both  rural  and  urban 
fuels,  the  5  backward-spread  values  can  be  converted  to  forward 
spread  values  by  using  ratios  previously  found  in  data  for  rural 
fuel. 


The  effect  of  wind  velocity  was  next  examined  by  plotting 
rate  of  spread  within  each  group  against  wind  velocity.  In  the 
most  populated  group,  IV-3,  the  8  dry  cases  show  an  upward  trend 
with  increasing  wind  velocity.  An  approximately  fitted  straight 
line  can  be  represented  by  the  equation 

R  =  0.027  +  1.25  x  10"3w 

where 

R  =  rate  of  spread 
w  wind  velocity 
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both  in  miles  per  hour.  This  equation  can  be  compared  with  its 
equivalent  in  rural  fuel: 

R  =  0.020  +  1.54  x  10"3 

for  w  >  10  mph  and  wind  direction  between  0  and  45  deg  (derived 
from  Table  3) . 

The  7  wet  cases  also  seem  to  indicate  an  upward  trend.  Un¬ 
fortunately,  all  the  7  are  grouped  in  a  narrow  range  of  wind 
velocity  and  the  scattering  is  such  that  the  rate  of  increase 
in  rata  of  spread  cannot  be  reliably  estimated.  In  this  case, 
as  in  the  dry  case,  the  scattering  of  data  points  can  be  reduced 
somewhat  by  rejecting  points  with  measured  time  less  than  1  hour 
(1  such  point  in  the  dry  case) .  Some  discussion  on  this  effect 
will  appear  soon. 

A15  the  3  groups  III-l,  III-2,  and  III- 3  also  show  a  defi¬ 
nite  upward  trend  although  with  more  scattering.  For  groups 
with  lower  fuel  densities,  the  dependence  on  wind  velocity  is 
completely  masked  by  some  other  effect.  We  shall  attribute  this 
fact  and  several  others  to  the  relative  extent  of  ground  spread 
and  spotting  which  varies  with  fuel  density. 

To  investigate  the  remaining  variable,  i.e.,  the  measured 
time,  it  is  necessary  to  keep  the  effect  of  wind  velocity  constant. 
We  have  seen  that  the  effect  of  wind  velocity  is  approximately  the 
same  in  urban  and  rural  fuel  at  least  in  the  cases  for  which  this 
effect  can  be  evaluated.  Therefore,  this  effect  can  be  kept  con¬ 
stant  in  analysis  of  urban  data  by  converting  rates  of  spread  at 
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any  wind  velocity  to  their  values  at  a  wind  velocity  of  10  mph, 
using  appropriate  ratios  previously  found  for  rural  fuel.  A 
velocity  of  10  mph  was  selected  as  basis  for  reference  because 
it  was  approximately  the  average  of  wind  velocities  recorded  in 
the  data.  Most  corrections  applied  were  less  than  50  percent. 

We  are  now  ready  to  examine  the  effect  of  measured  time. 
Values  in  each  group  were  plotted  against  the  time  during  which 
the  spread  had  been  observed.  For  the  largest  group,  i.e., 

IV-3,  most  points  are  concentrated  in  the  time  range  from  0~5  hr 
and  are  almost  uniformly  scattered.  Outside  this  time  range  there 
is  only  one  point,  at  10.2  hr,  which  is  slightly  lower  than  the 
average  of  points  in  the  0—5  hr  range.  However,  this  fact  was 
not  considered  sufficient  to  indicate  a  downward  trend  with  in¬ 
creasing  time  of  observation.  For  other  groups,  the  scattering 
also  tends  to  be  more  pronounced  at  short  times  of  observation, 
especially  when  all  groups  of  the  same  fuel  density  are  combined 
to  obtain  enough  data  points.  The  scattering  at  short  times  will 
be  accounted  for  in  a  subsequent  section. 

The  average  values  for  rate  of  spread  in  each  group  are  given 
in  Fig.  11.  Numbers  of  values  giving  the  averages  are  shown  in 
parentheses.  Average  values  shown  at  the  bottom  of  Fig.  11  are 
for  combined  fuel  types  in  each  fuel  density  range.  Averaging 
was  done  by  dividing  the  sum  of  distances  by  the  sum  of  corrected 
times  which  were  determined  from  the  corrected  rate  and  distance. 

We  have  said  that  7  of  the  15  values  in  group  IV-3  are  for 
wet  conditions.  (Roofs  and  exterior  walls  were  wet  due  to  rain 
or  melting  snow.)  The  average  rate  of  spread  for  the  wet 
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category  is  0.029  mph  and  that  for  the  dry  category  is  0.041. 

The  difference  seems  to  indicate  a  real  effect,  although  no 
further  evidence  is  available. 

Spotting 

Data  for  rate  of  spread  by  spotting  were  next  investigated. 
The  29  data  points  (including  those  with  measured  time  less  than 
1  hr  and  later  discarded)  were  grouped  according  to  fuel  type  and 
building  density  as  before.  There  was  no  case  of  backward  spread. 
The  wind  velocities  range  from  4  to  38  mph,  with  an  average  of 
17.6  mph,  whereas  in  ground  spread  data,  wind  velocities  vary 
from  3  to  28  mph  and  average  to  14.9  mph.  Thus,  the  occurrence 
of  spotting  does  not  seem  to  depend  strongly  on  high  wind  veloc¬ 
ity. 


Correction  for  wind  velocity,  as  has  been  done  with  ground 
spread  data,  is  not  possible  in  this  case  because  data  for  spot¬ 
ting  in  rural  fuel  have  not  been  thoroughly  studied,  and  also, 
the  effect  of  wind  on  spotting  may  not  be  the  same  for  both  types 
of  fuel.  Thus,  wind  velocity  and  measured  time  could  not  be 
separated . 

Rate  of  spread  values  were  plotted  against  measured  time 
for  each  fuel  density  range  and  for  all  fuel  types  combined. 
Again,  wild  fluctuations  occur  at  short  measured  times  (less 
than  1  hr).  However,  the  average  of  values  in  this  range  are 
consistently  higher  than  values  at  longer  times.  Thus,  there 
seems  to  be  some  tendency  to  report  cases  of  fast  jump  and  to 
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overlook  cases  in  which  the  fire  spent  some  time  before  jumping. 
To  minimize  this  bias,  values  with  times  of  measurement  less 
than  1  hr  have  been  discarded  (12  such  values) . 

Average  rates  of  spread  for  each  group  are  given  in  Fig.  12. 
Values  at  the  bottom  of  the  figure  are  averages  for  all  fuel 
types  in  each  fuel  density  range. 

DISCUSSION 

Many  results  of  this  analysis  can  be  interpreted  in  terms 
of  the  characteristics  of  urban  fuel  and  the  spread  processes 
in  this  medium. 

A  significant  characteristic  of  urban  fuel  is  discontinuity. 
Each  structure  is  surrounded  by  an  open  space  of  varying  dimen¬ 
sion,  in  contrast  with  the  essentially  continuous  forest  fuel. 

A  consequence  of  this  discontinuity  is  that  even  ground  spread 
proceeds  by  jumps.  Any  jump  may  or  may  not  take  place,  and  if  it 
does,  the  time  spent  before  the  jump  is  highly  variable,  depend¬ 
ing  on  the  distances  to  adjacent  structures,  the  orientation  of 
these  structures,  and  on  wind  velocity  and  direction.  Spotting 
is  still  more  stochastic.  For  an  ignition  to  occur,  a  firebrand 
must  be  still  burning  on  landing  and  must  land  on  a  surface 
which  happens  to  be  ignitable. 

Fuel  discontinuity  also  hinders  development  of  full  and 
uniform  fire  fronts  in  urban  areas.  Probability  considerations 
coupled  with  some  data  on  rate  of  spread,  burning  time,  and 
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Fig.  12.  Average  Values  for  Spotting  Spread  (mph) 
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fuel  distribution  seems  to  indicate  that,  except  for  extreme 
fuel  and  weather  conditions,  a  mass  fire,  at  any  given  time, 
is  made  up  of  a  number  of  small  islands  (a  few  structures) 
burning  simultaneously  and  propagating  independently.  Most 
fires  recorded  in  the  data  were  representative  of  these  islands. 
Since,  for  a  given  set  of  fuel  and  weather  conditions,  the  rate 
of  spread  in  nonextreme  cases  depends  on  the  size  of  individual 
islands  but  not  on  the  size  of  the  overall  conflagration,  rate 
of  spread  values  obtained  from  these  moderately  "big"  fires  can 
be  used  for  very  large  overall  fires  as  expected  during  nuclear 
attacks . 

When  a  spot  jump  takes  place,  the  fire  is  considered  to 
have  reached  the  newly  ignited  structure,  although  intermediate 
structures  will  burn  only  at  a  later  moment. 

Since  jumps  are  essentially  random  both  in  occurrence  and 
time,  it  is  easy  to  understand  that  measurements  on  short  sprea^ 
runs  involving  a  small  number  of  jumps  are  likely  to  yield 
widely  fluctuating  values  for  rate  of  spread  whereas  for  long 
runs,  the  overall  rate  of  spread  is  more  reproducible.  Con¬ 
sequently,  the  average  rate  in  each  group  was  obtained  by  divid¬ 
ing  total  distance  by  total  time  rather  than  by  adding  all  rate 
values  and  dividing  by  the  number  of  values.  In  this  manner, 
short  runs  are  combined  into  long  runs  and  fluctuations  are 
expected  to  be  smoothed  out. 

From  the  tables  of  average  rates  of  spread  for  both  ground 
spread  and  spotting  we  note  that,  within  each  fuel  type,  the 


rates  increase  with  decreasing  fuel  density.  Within  each  fuel 
density  range ,  however,  the  rates  do  not  show  a  definite  trend 
as  fuel  type  varies.  Consequently,  fuel  density  was  considered 
the  more  important  factor. 

The  increasing  rates  of  spread  when  fuel  density  decreases 
was  interpreted  as  due  to  the  increasing  number  of  spot  jumps 
compared  with  the  number  of  ground  jumps.  To  visualize  this 
process,  let  us  consider  a  case  of  very  high  fuel  density  and 
a  case  of  very  low  density.  If  the  density  a  was  near  unity, 
spacing  between  structures  would  be  practically  zero.  There¬ 
fore,  every  time  a  structure  burns,  at  least  one  adjacent  struc¬ 
ture  would  be  ignited,  i.e.,  the  probability  P  (a)  of  a  ground 

O 

jump  during  the  burning  time  would  be  one.  The  probability  of 
firebrands  being  emitted,  carried  a  given  distance,  and  remain¬ 
ing  active  on  landing  depends  on  the  type  and  volume  of  the  burn¬ 
ing  structure,  wind  velocity,  atmospheric  moisture,  and  so  on. 

For  the  moment,  let  us  assume  that  this  probability  is  not  zero 
and  denote  it  by  Pg .  Ignition  by  spotting  also  requires  that 
the  firebrand  lands  on  the  roofs,  so  that  the  probability  of  a 
spot  jump  is  a  function  of  a-  We  assume  that  this  probability 

is  proportional  to  a  and  we  represent  it  by  c.P  .  Thus,  in  a 

s 

fuel  area  of  very  high  density,  the  fire  never  misses  a  ground 
jump.  It  proceeds  from  A  to  B,  then  to  C,  and  so  on.  From  time 
to  time  it  also  makes  a  spot  jump  from  say  F  to  V.  If  the  fire 
is  observed  for  a  sufficiently  long  time,  the  average  rate  of 
spread  can  be  evaluated  by  counting  the  number  of  ground  jumps 
and  the  number  of  spot  jumps,  multiplying  these  numbers  by  their 
respective  distances  and  dividing  by  the  total  time. 
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As  the  fuel  density  a  decreases,  the  probability  of  spotting 

is  assumed  to  decrease  linearly  with  a.  As  to  the  probability 

of  ground  spread  P  (a),  evidence  suggests  that  the  decrease  is 

§ 

faster.  On  this  basis,  for  very  low  fuel  density,  a  large  per¬ 
centage  of  fires  would  die  out  soon  after  start.  But  those 
which  did  spread  to  some  distance  would  do  so  predominantly  by 
spotting.  Therefore,  high  rates  of  spread  would  be  observed. 

To  describe  this  process  with  a  mathematical  expression, 

let  us  consider  the  following  probabilities  for  each  burning 

structure:  the  probabilities  P  (a)  of  making  a  ground  jump  and 

§ 

[l  -  P  (a)]  of  not  making  it;  the  probabilities  aP  of  making  a 
g  s 

spot  jump;  and  [l  -  aP  ]  of  not  making  it.  All  possible  com- 

s 

binations  are: 


1  -  aP 


(a)  Pg(a) 

(b)  Pg(a) 

(c)  1  -  P  (a) 


(d)  1  -  P  fa)  1  -  aP 


(57) 


To  evaluate  overall  rate  of  spread,  spot  jumps  are  counted 

even  when  there  are  simultaneous  ground  jumps.  Therefore,  the 

probability  for  a  count  of  spotting  is  (a)  +  (c)  or  aP  ,  and  the 

s 

probability  of  a  count  of  ground  spread  is  (b)  or  P  (a) (1  -  aP  ) , 

g  s 

Given  that  the  fire  has  made  a  total  of  n  jumps,  the  number 
of  ground  jumps  is 

nPg(a)(l  -  aPs) 

P  (a)(l  -  aPs)  +  aPg 


(58) 
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and  the  number  of  spot  jumps  is 


naP 


Pg(a)(l  -  aPs)  +  aPs 


(59) 


The  average  distance  per  ground  jump  is  Va /a ,  where  a  is 
the  average  area  of  the  structures.  This  distance  includes  the 
spacing  and  the  average  dimension  of  the  structure  itself.  The 
distance  for  spot  jumps  is  highly  variable,  but  we  consider  a 
mean  value  d  . 


The  total  distance  of  spread  is 


m/a /a  P  (a)  (1  -  aP  )  +  nd  aP 

P  (a)  (1  -  aP  )  +  aP 
g '  s  s 


(60) 


The  time  per  jump  is  a  stochastic  variable  and  depends  on 
fuel  density.  We  assume,  however,  an  average  time  t^  for  both 
type  of  jumps.  Then  the  total  time  is  nt.,  and  the  average  rate 
of  spread  is  given  by 


_  Va/cT  P  (a)(l  -  aP  )  +  d  aP 

r  =  - ta _ _ 2 _ 2 _ 2. 

[P  (a) (1  -  aP  )  +  aP  ]t. 
g  s  sJj 


(61) 


For  very  high  fuel  density,  the  rate  of  spread  is  obtained 

from  this  expression  by  making  P  (a)  =  1.  Unless  P  is  zero, 

S  s 

the  fire  proceeds  by  both  ground  jumps  and  spot  jumps.  For  suf¬ 
ficiently  low  fuel  density,  P  (a)  =  0.  The  fire  then  proceeds 

S 

exclusively  by  spotting  and  the  rate  of  spread  is  simply  d  /t., 

®  J 

By  comparing  the  formulas  for  rate  of  spread  in  the  two 
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cases  it  is  easy  to  see  that  rate  of  spread  can  be  much  greater 
at  low  density  than  at  high  density. 

In  addition  to  the  increase  in  rate  of  spread  with  de¬ 
creasing  fuel  density,  several  other  results  can  also  be  ex¬ 
plained  by  this  simple  model. 

We  have  seen  that  more  spotting  fires  were  observed  in  low 
density  areas  than  in  high  density  areas,  whereas  for  ground- 
spread  fires,  the  reverse  was  true. 

We  also  noted  that  when  rates  of  spread  were  plotted  against 
wind  velocity,  the  degree  of  scattering  increased  rapidly  when 
going  from  one  density  range  to  the  next  lower  range,  and  for 
the  two  lowest  density  ranges,  the  scattering  completely  masked 
any  dependence  on  wind  velocity.  This  is  due  to  the  fact  that 
as  density  decreases,  the  relative  probability  for  spot  jumps 
increases.  Since  the  distance  per  spot  jump  is,  in  general, 
much  larger  than  the  distance  per  ground  jump,  the  average  rate 
of  spread  for  insufficiently  long  runs  is  greatly  affected  by 
the  occurrence  or  nonoccurrence  of  spot  jumps. 

The  effect  of  dryness  could  not  be  quantitatively  investi¬ 
gated  except  for  group  IV-3  in  which  wet  condition  seemed  to  de¬ 
crease  the  rate  of  spread  to  some  extent.  If  enough  data  points 
were  available,  it  would  be  interesting  to  see  the  effect  of  wet 
conditions  at  low  fuel  densities.  Exterior  wetness  is  expected 
to  decrease  the  rate  of  ground  spread  because  more  heat  is  re¬ 
quired  to  raise  the  temperature  of  exterior  walls  and  because  o^ 
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a  decrease  in  the  amount  of  pilot  ignition  due  to  embers  being 
extinguished.  Spotting  also  becomes  less  likely  if  roofs  are 
wet.  The  effect  of  wetness  on  rate  of  spread  is  probably  more 
drastic  with  spotting  than  with  ground  spread. 

Tables  2  and  3  refer  separately  to  fires  which  have  been 
labeled  as  ground-spread  or  spotting,  respectively.  Distances 
covered  during  each  observation  are  rather  short.  On  the  basis  of 
the  model  just  described,  both  types  of  spread  should  occur  in 
any  sufficiently  long  spread  run  according  to  their  own  proba- 
'ilities.  Averages  in  each  table  were  obtained  by  adding  the 
distances  and  dividing  by  the  total  time,  thus  combining  many 
short  fires  to  a  long  one.  If  we  now  combine  Tables  2  and  3  in 
the  same  way,  the  resulting  data  would  characterize  long,  spread¬ 
ing  fires  and  would  give  better  average  velocities  for  fires 
which  spread  indefinitely  if  fuel  conditions  allow  them  to 
spread.  The  combined  averages  are  0.28,  0.29,  0.07,  and  0.03 
for  fuel  densities  I,  II,  III,  and  IV  (in  that  order)  for  all 
fuel  types  in  each  density  range. 

Some  further  comments  on  these  final  results  seem  to  be 
relevant.  The  apparent  constancy  of  the  first  two  figures  can 
be  explained  by  various  arguments:  it  may  be  a  real  effect, 
i.e.,  the  mean  rate  of  spread  may  level  off  as  fuel  density 
decreases  below  II;  data  for  density  range  l  may  be  insufficient, 
hence  the  first  value  is  less  reliable  than  others;  actual,  fuel 
density  for  range  L  may  not  be  much  less  than  20  percent,  a  value 
which  would  put  these  cases  in  range  LI,  etc. 
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As  mentioned  previously,  only  the  most  important  factor, 
namely,  fuel  density,  could  be  studied  quantitatively.  The 
effect  of  exterior  humidity  is  real  but  could  not  be  evaluated. 

The  final  values  for  mean  rate  of  spread  characterize  peacetime 
fires  spreading  through  significant  distance.  These  fires  may 
have  occurred  under  the  more  severe  fire  weather  conditions  since 
it  is  easier  to  put  out  fires  under  poor  fire  weather  conditions. 
Fires  from  nuclear  attacks  will  spread  under  weather  conditions 
that  are  characteristic  of  the  locality,  the  year,  and  the 
season.  If  the  effect  of  weather  on  rate  of  spread  is  sufficiently 
important,  the  sample  of  peacetime  fires  may  contain  more  fast 
fires  than  does  the  population  of  fires  from  nuclear  attacks, 
therefore  the  mean  rates  of  spread  given  above  may  be  on  the 
high  side  when  used  for  nuclear  attack  fires.  It  is,  of  course, 
impossible  to  make  any  correction  for  this  fact.  However,  for 
conservative  prediction  of  rate  of  spread  for  fires  from  enemy 
attacks,  the  uncertainty  may  not  be  of  great  consequence. 
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Section  11 

VARIOUS  ESTIMATED  DATA 

In  this  section  we  group  various  types  of  information  derived 
from  very  limited  observed  data,  supplemented  in  some  cases  by 
subjective  judgment. 

SPREAD  PROBABILITY  IN  URBAN  FUEL 

The  U.S.  Strategic  Bombing  Survey  studied  the  efficiency  of 
37  linear  miles  of  fire  breaks  of  various  widths  in  Nagoya,  strong 
winds  or  fire  storms.  Fire  was  stopped  over  34.8  percent  of  the 
total  length  of  breaks  in  the  width  range  from  65  to  150  ft  and 
75  percent  of  the  total  length  for  widths  ranging  from  150  ft 
upwards  (Ref,  2).  This  seems  to  indicate  that  for  a  spacing  of 
107  ft  (midpoint  between  65  and  150  ft),  the  probability  of  spread 
is  approximately  0.65,  and  for  some  spacing  beyond  150  ft,  it  is 
0.25.  This  spacing  is  probably  not  far  from  the  point  at  200  ft 
i.e.,  halfway  between  150  and  250  ft.  The  latter  spacing  is  con¬ 
sidered  as  the  limit  for  fire  spread  and  is  justified  only  by  a 
feeling  acquired  through  personal  experience  and  training. 

These  two  points  and  the  fact  that  fire  spread  is  certain 
when  spacing  is  zero  allow  us  to  draw  curve  A  in  Fig.  13.  It  must 
be  emphasized  that  this  curve  should  be  considered  as  an  upper 
limit  of  the  probability  of  spread  for  two  reasons:  first,  fires 
were  not  considered  stopped  in  places  where  burned  structures 
were  found  on  both  sides  of  the  fire  breaks,  although  both 
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Fig.  13.  Probability  of  Spread  Versus  Distance 
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sides  of  an  adequate  fire  break  might  have  been  hit  directly  by 
incendiary  bombs;  second,  fire  might  have  jumped  across  a  fire 
break  ai_  one  point  and  spread  along  some  distance  on  the  other 
side  of  the  break.  The  stopping  efficiency  of  the  break  would 
be  underestimated  if  fire  burned  on  both  sides  —  at  points  B,  C, 

D,  E  on  one  side  and  at  corresponding  points  B' ,  C’  ,  D' ,  E'  on 
the  other  side  —  but  the  fires  at  these  latter  points  originated 
from  A'  instead  of  B,  C,  D,  E. 

Examination  of  fire  spread  data  in  incendiary-attacked 
German  cities  has  led  to  the  conclusion  that  under  normal  fire 
conditions,  a  10-ft  space  between  two  brick  buildings  had  about 
50  percent  change  of  preventing  fire  spread  (Ref.  2).  The  prob¬ 
ability  curve  for  this  special  case  (spread  between  two  brick 
buildings)  is  perhaps  not  too  far  off  from  curve  B  in  Fig.  13, 
which  is  obviously  a  lower  limit  of  the  curve  being  sought, 
since  brick  buildings  are  generally  less  vulnerable  to  fire  spread 
than  other  types  of  buildings  commonly  existing  in  urban  areas. 

The  solid  curve  in  this  figure,  drawn  between  these  two 
limits,  represents  the  best  guess  and  is  believed  to  yield  ac¬ 
ceptable  results  when  applied  to  urban  fuel  in  general. 

WEATHER  CONDITIONS  FOR  NO-SPREAD  IN  RURAL  FUEL 

The  list  of  weather  conditions  for  no-spread  in  rural  fuel 
and  the  sets  of  data  following  have  been  provided  by  the  U.S. 
Forest  Service  (Ref.  1). 
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All  Fuels.  Over  1  in.  of  snow  on  the  ground  at  the  nearest 
weather  reporting  stations. 

Grass .  Relative  humidity  above  80  percent. 

Brush  or  Hardwoods.  0.1  in.  of  precipitation  or  more  with¬ 
in  the  past  7  days  and  wind 

•  0—3  mph;  relative  humidity  60  percent  or 
higher,  or 

•  4—10  mph;  relative  humidity  75  percent  or 
higher ,  or 

•  11—25  mph;  relative  humidity  85  percent  or 
higher. 

Co;  •ftfM :  Timber.  (a)  1  day  or  less  since  at  least  0.25  in. 
of  precipitation  and  wind 


•  0—3  mph;  relative  humidity  50  percent  or 
higher,  or 

•  4~10  mph;  relative  humidity  75  percent  or 
higher,  or 

•  11—25  mph;  relative  humidity  85  percent  or 
higher. 

(b)  Or,  2—3  days  since  at  least  0.25  in.  of  precipita¬ 
tion  and  wind 


•  0—3  mph;  relative  humidity  60  percent  or 
higher,  or 

•  4—10  mph;  relative  humidity  80  percent  or 
higher,  or 

•  11—25  mph;  relative  humidity  90  percent  of 
higher. 


UR S  641-6 


11-5 


(c)  Or,  4—5  days  since  at  least  0.25  in.  of  precipita¬ 
tion  and  wind 

•  0—3  mph;  relative  humidity  80  percent  or  higher. 

(d)  Or,  6—7  days  since  at  least  0.25  in.  of  precipita¬ 
tion  and 

•  0—3  mph;  relative  humidity  90  percent  or  higher. 
WEATHER  CONDITIONS  FOR  EXTINCTION 

Rural  Fires 


Extinction  after  total  burning  time  is  not  considered. 

Grass .  Measurable  precipitation  at  the  three  nearest 
weather  stations. 

Brush  or  Hardwoods.  0.1  in.  of  precipitation  or  more  at 
the  three  nearest  weather  stations. 

Conifer  Timber.  (a)  0,5  in.  of  precipitation  or  more  at 

the  three  nearest  weather  stations. 

(b)  Or.  0.25  to  0.5  in.  of  precipitation  at  the  three 
nearest  weather  stations  and  no-spread  conditions  for  the  fol¬ 
lowing  two  12-hr  periods. 

(c)  Or,  no-spread  conditions  for  eight  consecutive  12-hr 
periods  and  measurable  precipitation  at  the  three  nearest  weather 
stations  during  any  two  12-hr  periods. 
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Urban  Fires 


It  is  noticed  that  extinction  of  urban  fires  requires  rather 
severe  weather  conditions. 

Light  Residential.  1.0  in.  of  precipitation  at  the  Weather 
Burean  Station  and  no-spread  conditions  for  36  consecutive  hr. 

Heavy  Residential.  1.5  in.  of  precipitation  at  the  city 
Weather  Bureau  Station  and  no-spread  conditions  for  72  consecu¬ 
tive  hr. 

Commercial.  2.0  in.  of  precipitation  at  the  city  Weather 
Bureau  Station  and  no-spread  conditions  for  seven  consecutive 
days . 

City  Center  or  Massive  Manufacturing.  2.0  in.  of  precipita¬ 
tion  at  the  city  Weather  Bureau  Station  and  no-spread  conditions 
for  two  consecutive  months. 

BURNING  TIMES 


Total  burning  time  is  defined  as  the  period  during  which  a 
large  fire  might  remain  stationery  yet  be  capable  of  resuming 
active  spread  if  burning  conditions  changed  for  the  worse.  Burn¬ 
ing  times  for  various  classes  of  fuel  are  listed  below. 

Rural  Fuels 


Fuel  Type 


Total  Burning  Time 


Grass  30  min 

Light  brush  16  hr 

Medium  brush  36  hr 

Heavy  brush  72  hr 

Timber  7  days 
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Urban  Fuels 


Fuel  Type 

Total 

Burning 

Light  residential 

36  hr 

Heavy  residential 

72  hr 

Commercial 

7  days 

City  center  and  massive  mfg. 

2  mos . 

Time 


1 
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Section  12 

FIRE  SPREAD  VULNERABILITY  OF  RESOURCES 

The  general  approach  to  the  problem  of  fire  damage  assessment 
involves,  first,  estimation  of  the  general  burnout  of  fuel  areas 
and,  second,  determination  of  the  fire  damage  to  resources  by  means 
of  fire  spread  vulnerability  factors.  This  section  deals  with  the 
second  aspect  of  the  problem. 

Fire  spread  vulnerability  can  be  defined  in  simple  terms  as 
follows:  Given  that  the  general  area  is  swept  by  fire,  is  the 
resource  destroyed  or  what  is  its  chance  of  being  destroyed? 

The  solution,  however,  is  made  difficult  by  the  limited  amount 
of  appropriate  data.  In  this  study,  methods  have  been  selected 
which  can  be  used  for  evaluating  vulnerabilities  on  the  basis  of 
some  existing  data  and  which  can  be  readily  improved  as  more 
extensive  and  more  accurate  data  become  available. 

The  vulnerability  of  resources  to  fire  spread  clearly  depends 
on  three  groups  of  factors:  susceptibility  of  the  resource  to 
fire  (normally  involving  the  structural  characteristics  and  con¬ 
tents),  the  characteristics  of  the-  adjacent  fuel  (as  an  ignition 
source),  and  the  spacing  between  the  resource  and  the  adjacent 
fuel. 


It  seems  advantageous  to  study  first  the  vulnerability  for 
one  type  of  resource  in  order  to  establish  the  effects  of  spacing 
and  surrounding  fuel  on  the  ignition  event.  The  result  is  a 
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measure  of  that  part  of  vulnerability  due  to  the  environment  of 
the  resource.  Attention  is  next  given  to  the  question  of  how 
the  vulnerability  changes  with  the  structural  characteristics 
of  the  resource. 

The  type  of  resource  chosen  for  studying  the  environment 
dependence  of  vulnerability  is  the  common  type  of  structure  in 
urban  areas.  These  structures  are  characterized  by  substantial 
combustible  surface  exposed  to  heat  and  firebrands.  The  reason 
for  this  choise  is  that  most  data  on  ignition  events  refer  to 
this  type  of  exposed  structures. 

With  respect  to  environment  or  type  of  fuel  in  adjacent 
areas,  three  situations  may  be  distinguished:  (1)  the  fuel  is 
entirely  or  predominantly  urban;  (2)  the  fuel  is  entirely  or 
predominantly  rural;  and  (3)  the  fuel  is  an  approximately  equal 
mixture  of  rural  and  urban  types.  In  the  following  subsections, 
we  shall  successively  discuss  the  urban,  rural,  and  mixed  en¬ 
vironments  . 

URBAN  ENVIRONMENT 

The  characteristics  of  surrounding  fuel  can  be  specified 
in  three  ways : 

A0  The  size,  structural  materials  of  each  structure, 
and  its  distance  to  the  resource  are  given. 

B.  Only  the  number  of  adjacent  structures  and  their 
distances  to  the  resource  are  specified. 

C.  The  only  information  is  that  the  resource  is  within 
an  urban  area  of  a  given  class  and  given  building 
density. 
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Case  A 

When  all  characteristics  of  surrounding  fuel  are  completely 
specified,  a  deterministic  approach  is  obviously  indicated.  Two 
alternative  methods  are  proposed.  In  one  method,  the  ignition 
event  is  predicted  on  the  basis  of  the  intensity  of  radiated 
heat,  adjusted  upward  to  include  the  estimated  effect  of  con- 
vected  heat  and  embers.  The  other  method  uses  a  formula  which 
gives  a  plausible  measure  of  total  heat  and  embers  evolved  from 
the  burning  structure,  the  minimum  measure  required  for  ignition 
being  derived  empirically  or  subjectively. 

It  will  appear  that  either  method  can  be  improved  by  fur¬ 
ther  observed  data  as  they  become  available..  Since  it  is  not 
certain  which  one  will  be  first  to  receive  backup  data  and  be¬ 
come  more  reliable,  this  report  will  describe  both. 

Method  of  Adjusted  Radiation  Intensity 

The  radiant  energy  as  a  function  of  distance  from  a  radi¬ 
ating  surface  can  be  measured  or  calculated  if  the  source  char¬ 
acteristics  and  geometry  are  well  specified.  The  variation  of 
conveoted  heat  with  distance  presents  a  more  complicated  prob¬ 
lem  and  has  not  been  thoroughly  investigated.  Flame  contact 
and  embers  (falling  at  short  distance)  are  predominantly  random. 
Their  treatment  as  deterministic  is  justified  by  the  following 
consideration.  Since  at  distances  within  which  these  flames 
and  embers  are  active  the  fuel  is  preheated  by  radiated  and 
convected  heat,  ignition  requires  only  a  little  tongue  of  flame 
or  one  or  two  embers  to  fall  on  a  piece  of  fuel.  Thus,  the 
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action  of  flames  and  embers  is  primarily  to  extend  the  effect 
of  radiated  and  convected  heat  to  distances  at  which  this  heat 
could  not  by  itself  produce  ignition.  Spotting  firebands  fall¬ 
ing  beyond  the  heat  range  is  a  purely  stochastic  phenomenon. 

Since  nothing  is  known  about  spotting  except  that  its  probabil¬ 
ity  is  usually  low  compared  with  that  of  continuous  spread,  we 
leave  spotting  outside  the  scope  of  this  study. 

Several  investigators  (Refs.  3,  4,  5)  have  reported  their 
results  on  radiated  heat.  At  any  distance  c  on  the  normal  pass¬ 
ing  through  the  center  of  a  rectangular  radiating  surface,  the 
ratio  of  the  intensity  I  to  the  intensity  IQ  (at  zero  distance 
from  the  source)  is  given  by  the  so-called  configuration  factor 
§,  which  depends  not  only  on  the  distance  c  but  also  on  the 
shape  of  the  radiating  surface,  i.e.,  on  the  ratio  N  of  the 
longer  dimension  to  the  shorter  dimension.  In  Fig.  14  (from 
Ref.  5),  the  factor  cp,  which  is  one-fourth  ?,  is  plotted  against 
NA/c,  where  A  is  half  the  shorter  dimension.  Thus  cp  is  the  con¬ 
tribution  from  each  of  the  four  sections  of  the  radiating  sur¬ 
face  obtained  by  drawing  two  lines  through  the  center  of  the 
surface,  one  lengthwise  and  the  other  crosswise. 

The  factors  §  and  cp  apply  to  fully  radiating  surfaces. 

When  a  surface  is  made  up  of  radiating  sections  (windows,  com¬ 
bustible  walls)  alternating  with  nonradiating  sections,  an  exact 
treatment  would  require  each  radiating  section  to  be  considered 
separately  and  the  total  effect  obtained  by  summing  ail  contribu¬ 
tions.  In  practice,  however,  it  is  much  simpler  and  sufficiently 
accurate  to  multiply  I  by  the  fraction  of  total  area  represent¬ 
ing  window  and  combustible  area. 


figuration  Factor  cp  for  Differently  Shaped  Radiators 
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One  more  correction  is  needed  to  account  for  the  fact  that 
the  area  actually  radiating  at  any  time  is  not  necessarily  the 
entire  face  of  the  structure.  If  one  horizontal  dimension  is 
much  larger  than  the  other,  the  fire  may  well  burn  in  only  one 
section  at  a  time.  Obviously,  there  is  no  fixed  value  for  the 
area  of  the  radiating  section,  but  it  is  believed  that  a  reason¬ 
able  maximum  value  for  typical  cases  is  the  product  of  the 
height  and  twice  the  depth  (assumed  to  be  the  minimum  horizontal 
dimension)  of  the  structure.  If  the  width  of  the  radiating  sur¬ 
face  is  less  than  twice  the  depth  of  the  structure  or  if  the 
radiating  surface  is  along  the  depth  of  the  structure,  no  cor¬ 
rection  is  needed. 

The  intensity  I  at  the  radiating  surface  depends  on  the 

o 

fire  load  or  total  amount  of  -combustible  in  the  enclosure.  For 

-2 

fire  loads  greater  than  about  5  ib.-ft  ,  a  common  case,  the  in- 

-2  -1 

tensity  is  about  4  cal-cm  ~-see  .  Smaller  fire  loads  produce 
intensities  of  about  2  cal-cm  ^-sec  K 

The  minimum  heat  flux  required  to  ignite  wood  has  been 

-2  -1 

measured  (Ref.  5).  At  0,8  cal-cm  -sec  or  more,  spontaneous 

ignition  occurs  within  about  2  min.  When  the  heat  flux  is  be- 

-2-1 

tween  0.8  and  0.3  cal-cm  -sec  ,  pilot  ignition  can  take  place, 
if  after  5  —  10  min  of  heating,  a  little  flame  or  ember  comes 
near  the  heated  surface. 

There  are  no  data  on  the  amount  of  convected  heat  and  its 
variation  with  distance  or  on  the  number  of  embers  given  off 
by  a  structure  of  a  given  size.  Consequently,  their  effects 
must  be  taken  into  account  in  somewhat  arbitrary  ways:  first, 
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the  intensity  at  zero  distance  is  increased  from  4  to  6  cal-cm-  - 

-1  *  -2 
sec  for  fire  loads  equal  to  or  greater  than  5  lb-ft  and  from 

-2  -1 

2  to  3  cal-cm  -sec  for  smaller  fire  loads,  Second,  flames  and 

embers  are  assumed  to  be  always  available  for  pilot  ignition  at 

-2  -1 

distances  for  which  heat  intensity  exceeds  0.30  cal-cm  -sec  , 
provided  the  exposed  structure  has  some  significant  ignitable 
area.  With  the  further  assumption  that  -the  flux  of  convected 
heat  varies  with  distance  in  the  same  manner  as  the  radiated  com¬ 
ponent,  it  is  possible  to  evaluate  the  distance  cmax  for  each 

ignition  source  at  which  the  heat  intensity  is  equal  to  the  cr.it- 

-2  -1 

ical  value  0.30  cal-cm  -sec 

Method  of  Vulnerability  Number 

First,  a  function  Vg  is  selected  to  represent  the  igniting 
capability  of  the  ignition  source  at  various  distances ^  Then, 
the  minimum  value  this  function  must  have  for  ignition  is  empiri¬ 
cally  determined.  For  a  very  small  source,  an  appropriate  func¬ 
tion  is  probably 


(62) 


where  V  is  the  volume  of  the  ignition  source  (in  cubic  feet)  and 
d  the  separation  (in  feet). 


While  the  inverse  square  relation  is  probably  correct  for 
very  small  V,  there  are  two  errors  when  V  does  not  satisfy  this 
requirement.  The  first  error  is  of  a  geometrical  nature.  Since 
the  distance  d  is  measured  from  the  center  of  the  face  of  the 
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ignition  source  to  some  point  P  on  the  normal  to  the  center,  the 
simple  formula  is  correct  only  for  the  volume  element  AVo  exactly 
at  the  center  of  the  face.  For  all  other  volume  elements,  the 
true  distance  to  point  P  is  greater  than  d  by  an  amount  which 
increases  as  the  volume  elements  get  further  away  from  the  center 
of  the  face.  The  correct  value  of  Vg  for  point  P  is  given  by 


V 

e 


=  4 


z  D  ,-W/2  .  H/2 

;  i 

•  ! 

I  !  , 

0  0  O 


_ dx  dy  dz _ 

,  ,  ,  s  2  ,  2  ,  2 

(d  +  x)  +  y  +  z 


(63) 


where  D,  W,  and  H  are  the  depth,  the  width,  and  the  height  of  the 
ignition  source,  respectively,  and  x,  y,  z  are  distances  in  the 
direction  of  D,  W,  H.  Since  exact  evaluation  of  this  integral  is 
complicated,  the  following  approximate  formula  is  preferred, 


D  ...  -1  W  .  -1  H, 

V  =  -.77- — 77  tan  777  tan  77 
e  d(d  +  D)  2d  2d 


(64) 


V 


d  (1+7  WH 

d 


.  -1  -W  „  -1  H 

tan  2d  tan  2d 


The  second  error  is  due  to  the  fact  that  heat,  embers,  and 
firebrands  from  one  layer  of  the  source  may  be  partially  blocked 
by  all  layers  lying  in  front.  The  extent  of  this  effect  may  vary 
with  each  process,  i.e.,  radiation,  convection,  embers ,  -firebrands , 
with  the  fuel  arrangement  in  the  source,  and  with  the  phase  of 
combustion.  A  quantitative  correction  is  not  possible  at  the 
present  state  of  the  art.  Formula  (64),  ‘however,  will  be  on  the 
safe  side  (overestimate)  when  predicting  fire  spread  vulnerability. 
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Serious  error  may  result,  however,  from  using  Formula  (64) 
if  the  outside  walls  (or  a  significant  fraction  thereof)  are 
incombustible  and  do  not  collapse  until  the  end  of  the  combustion, 
because  in  this  case  all  or  part  of  the  igniting  flux  is  blocked 
off.  It  is  reasonable  to  assume  that  the  fraction  of  the  "igniting 
flux"  that  gets  through  is  equal  to  the  fraction  of  window  and 
combustible  area.  Denoting  this  fraction  by  a ,  we  write  the 
formula  as 


V 

e 


(65) 


Finally,  we  must  account  for  the  fact  that  a  structure  ignited 

in  the  normal  way  frequently  burns  in  only  one  section  at  a  time, 

especially  when  one  horizontal  dimension  is  much  larger  than  the 

other.  As  in  the  method  of  adjusted  ’adiation  intensity,  we  may 

limit  the  width  W  of  the  radiating  face  to  twice  the  depth  D  of 

the  structure,  so  that  the  maximum  volume  of  the  active  ignition 
2 

source  is  2HD  . 

The  minimum  V value  required  for  ignition  cannot  be  estab¬ 
lished  objectively,  because  observed  data  are  still  very  scanty. 
Hence,  only  a  subjective  estimate  is  available  at  present.  For 
an  ignition  source,  let  us  consider  a  one-story,  50  by  50  ft  by 
12-ft  high  structure  with  50  percent  window  area.  On  the  basis 
of  intuition  and  limited  data  on  urban  fires  (Ref.  2),  one  would 

be  inclined  to  believe  that  ignition  takes  place  at  a  distance  of 

-3 

30  ft  or  less.  The  corresponding  value  of  Vg  is  1.4  x  TO 

With  either  of  these  two  methods,  the  destruction  or  sur¬ 
vival  of  the  resource  is  predicted  by  adding  the  heat  intensities 
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or  the  V  values  from  all  adjacent  structures  acting  on  each  face 
of  the  resource  structure  and  comparing  the  highest  total  heat 
intensity  of  Vg  value  on  any  single  face  with  minimum  heat  or 
required  for  ignition. 

Addition  of  heat  intensities  or  Vg  values  from  adjacent  struc¬ 
tures  is  based  on  the  assumption  that  these  structures  burn  simul¬ 
taneously.  This  leads  to  the  question:  which  ones  will  burn 
simultaneously?  This  depends  on  the  course  of  the  fire  front. 
Suppose  one  face  of  the  resource  is  exposed  to  a  parallel  row  of 
structures  A,  B,  C,  .  .  ..  All  these  structures  will  burn  simul¬ 
taneously  if  a  fire  front  approaches  these  structures  in  a  direction 
perpendicular  to  the  row.  However,  this  is  not  necessarily  true 
if  the  fire  spreads  along  the  row  of  structures.  In  this  case  the 
number  of  structures  burning  simultaneously  must  be  determined 
from  the  rate  of  spread  and  the  burning  time  of  fuel.  Unless  a 
separate  vulnerability  determination  is  made  for  each  of  the  two 
cases,  it  is  recommended  that  the  average  number  of  simultaneous 
burning  structures  be  used. 

Case  B 


In  this  case  the  number  of  adjacent  structures  and  their  dis¬ 
tances  to  the  resource  are  specified.  Since  the  structural  char¬ 
acteristics  are  unknown,  it  is  necessary  to  assume  that  the 
distribution  of  these  characteristics ,  i.e.,  size,  construction 
materials,  number  of  stories,  etc.,  is  typical  of  most  urban 
areas.  Hence  the  vulnerability  problem  becomes  a  stochastic  one. 
Given  that  all  surrounding  structures  burn  (simultaneously  or  not), 
we  must  find  the  probability  that  the  resource  is  ignited.  Needed 
is  the  probability  that  in  a  typical  urban  area  a  burning  structure 
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ignites  another  structure  at  a  specified  distance.  Limited  infor¬ 
mation  of  this  type  has  been  given  in  the  previous  section. 

The  stochastic  vulnerability  of  the  resource  is  the  proba¬ 
bility  that  the  resource  is  ignited  by  one  or  more  surrounding 
structures.  If  the  spacings  are  d^,  d^,  d^,  .  .  .  and  the  corre¬ 
sponding  probabilities  of  ignition  are  P(d^),  P(d2),  P(d^),  .  .  ., 
as  given  in  Fig.  13,  the  vulnerability  of  the  resource  is  given  by 

1  -  [1  -  P(d1)][l  -  P(d2)][l  -  P(d3)]  .  .  .  (66) 

where  the  subscripts  of  d  extend  to  all  surrounding  structures 
(except  those  which  will  certainly  not  burn  in  case  of  general 
fire) . 


Case  C 

Not  only  structural  characteristics  but  spacings  also  are 
left  unspecified.  In  order  to  apply  the  data  of  Fig.  13,  we  must 
independently  determine  the  distribution  of  spacings  between 
adjacent  structures  by  studying  many  maps  of  urban  areas.  This 
work  is  best  done  by  dividing  urban  areas  into  a  number  of  cate¬ 
gories,  for  example:  general  residential,  residential  tract, 
commercial,  industrial,  etc.  Suppose  the  spacing  distribution 
is  expressed  as  number  (integral)  of  structures  at  each  spacing 
range  from  a  given  structure  and  denoted  by  n^,  n^,  n^,  .  .  . 
for  the  ranges  0-5,  5  —  10,  10—15,  .  .  .  ft,  respectively.  The 

probability  that  the  resource  is  ignited  by  one  or  more  structures 

."th¬ 
in  the  i  range  is 


P.  =  i  -  [i  -  p(i)i 


ni 


(67) 
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where  p(i)  is  the  average  ignition  probability  for  the  i  range, 
a  quantity  easily  determined  from  Fig.  13.  The  vulnerability  is 
then  given  by 


or 


1  -  n  (1 

i=l 


1 


co  n  • 

n  [i  -  p(i)] 


i=l 


(68) 


(69) 


RURAL  ENVIRONMENT 

We  are  here  concerned  with  an  area  in  which  urban  fuel  is 
insignificant  in  comparison  with  rural  fuel.  The  resource  is 
probably  separated  from  the  rural  fuel  by  a  zone  of  open  space, 
parking  areas,  roads,  lawns,  etc.  The  rural  fuel  surrounding 
this  zone  may  be  continuous  or  may  consist  of  patches  of  various 
sizes  separated  by  open  space.  Since  many  types  of  data  for  a 
stochastic  treatment  in  this  case  are  not  available,  fuel  data 
must  be  completely  specified,  for  example,  by  use  of  a  map  showing 
the  resource  and  the  surrounding  fuel. 

As  for  urban  environment  above,  the  heat  intensity  on  each 
face  of  the  resource  structure  is  calculated  by  means  of  either 
of  the  two  methods  described  in  those  paragraphs,  and  the  deter¬ 
ministic  vulnerability  of  the  structure  is  decided  on  by  consider¬ 
ing  the  face  which  receives  the  highest  intensity. 

In  applying  these  methods,  one  must  know  the  height  of  the 
radiating  face,  i.e. ,  the  probable  height  of  flames  from  burning 
rural  fuel,  the  dimensions  of  the  source  (size  of  fire  burning  at 
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the  same  time),  and  the  heat  intensity  at  the  radiating  surface 
Iq.  Flame  height  can  be  estimated  by  experienced  personnel  in 
terms  of  the  height  of  the  fuel.  The  width  and  depth  of  the 
source  must  be  determined  for  each  of  two  extreme  cases  (fire 
spread  perpendicular  or  parallel  to  a  given  face  of  the  resource), 
as  already  discussed  at  the  end  of  Case  A. 

The  intensity  of  radiated  heat  I  can  be  reasonably  estimated 
-2-1  0 

at  2  cal-cm  -sec  on  the  basis  of  experiments  on  structures 

(Ref.  5),  since  for  rural  fuel  the  ventilation  is  complete  and 

2 

the  fire  load  is  normally  less  than  5  lb/ft  .  For  medium  brush 

2 

the  fire  load  is  24  tons  per  acre  or  approximately  1  lb/ft 
(Ref.  1). 

MIXED  ENVIRONMENT 

When  comparable  proportions  of  urban  and  rural  fuels  exist 
in  the  environment,  a  deterministic  approach  seems  to  be  the  only 
feasible  one,  and  data  specification  must  be  such  that  this 
approach  can  be  used. 

Either  of  the  two  methods  of  Case  A  is  applied  to  each  struc¬ 
ture  and  each  patch  of  rural  fuel,  and  the  overall  effect  on  each 
face  of  the  resource  structure  is  computed.  The  computation  is 
likely  to  be  lengthy  in  some  cases.  It  might  be  possible  to  treat 
all  structures  irradiating  a  given  face  of  the  resource  as  a 
single  source  and  all  patches  of  rural  fuel  as  another  single 
source.  Appropriate  factors  (similar  to  window  factors)  must  be 
chosen  to  account  for  the  fact  that  only  certain  sections  of  the 
overall  source  constitute  the  actual  radiating  surface. 
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Section  13 

WEATHER  DATA  FOR  FIRE  SPREAD  PREDICTION 

SPECIFICATION  OF  THE  DATE  OF  ATTACK 

Since  fire  spread  is  affected  by  weather  and  since  for  a 
given  area  the  weather  varies  with  the  day,  the  month,  and  the 
year,  the  extent  of  fire  spread  predicted  clearly  depends  on  the 
time  of  attack  specified.  There  are  three  ways  of  specifying  the 
time  of  attack:  (1)  a  specific  date  in  the  past  is  selected; 

(2)  the  current  date  is  selected;  (3)  any  date  in  the  future  is 
considered  an  equally  likely  date  of  attack. 

ACQUISITION  OF  WEATHER  DATA 

Case  I 

A  specific  date  in  the  past  may  be  selected  for  one  of  three 
reasons:  (a)  at  that  date  and  during  the  subsequent  interval, 

weather  conditions  were  typical  for  the  area  of  interest;  (b)  the 
weather  had  some  characteristics  desirable  for  the  enemy  when 
planning  an  attack;  (c)  information  on  the  fire  damage  due  to  a 
certain  type  of  weather  is  wanted. 

In  this  case,  weather  data  do  not  present  a  major  problem. 

The  U.S.  Department  of  Commerce  Weather  Bureau  provides  detailed 
information  on  all  weather  factors  for  most  cities  and  surrounding 
areas.  Weather  data  for  the  specified  period  and  specified  area 
can  be  used  directly.  Some  handling  is  needed,  however ,  if  the 
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records  give  more  details  than  can  be  used  in  the  calculation. 

For  example,  in  many  cases  hourly  wind  velocity  and  direction  are 
given.  The  calculation  for  a  large  area  or  a  long-lasting  fire 
would  be  to  lengthy  if  new  data  for  wind  velocity  and  direction 
were  introduced  every  hour.  In  such  cases  the  wind  data  must  be 
reduced  to  average  values  for  day  and  night  so  that  calculation 
need  not  be  interrupted  for  intervals  of  approximately  12  hr. 

Case  II 

Selection  of  the  current  date  as  date  of  attack  may  be 
prompted  by  reasons  similar  to-  tl  ise  mentioned  in  the  previous 
case.  Data  on  current  weather  can  be  obtained  from  various 
weather  stations.  Weather  for  the  immediate  future  must  be  pre¬ 
dicted  on  the  basis  of  current  weather.  The  reliability  of  such 
prediction  depends  on  the  meteorological  characteristics  of  the 
area  and  especially  on  how  far  into  the  future  the  prediction  is 
carried  out.  At  present,  long-range  forecasts  can  cover  a  period 
of  5  days  with  reasonable  accuracy. 

Case  III 

The  date  of  attack  is  most  likely  to  be  left  indeterminate 
in  the  future.  The  task  in  this  case  is  to  estimate  the  mean 
fire  effect  on  the  assumption  that  the  attack  will  be  completely 
random  in  time. 

Whereas  current  weather  can  serve  as  a  basis  for  making  daily 
forecasts  or  forecasts  for  a  future  period  of  a  few  days,  long- 
range  prediction  requires  a  statistical  approach  based  on  weather 
that  has  been  observed  over  a  long  period  in  the  past. 
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One  obvious  way  to  evaluate  the  mean  fire  effect  would  be  to 
carry  out  a  fire  spread  calculation  for  each  day  in  the  past  taken 
as  attack  day  and  for  a  period  long  enough  to  cover  all  types  of 
weather  that  can  occur  in  the  area.  The  fire  effect  would  then 
be  averaged  over  all  attack  days.  However,  the  level  of  effort 
involved  in  this  calculation  is  clearly  prohibitive. 

A  simple  approach  is  to  select  2  months  with  extreme  weather 
by  inspection  of  a  period  of  n  years  and  calculate  fire  spread  for 
each  month  according  to  proceuures  to  be  given  later  in  connection 
W" uh  the  use  of  weather  data  for  fire  spread  prediction.  The 
result  would  give  an  upper  and  a  lower  limit  for  fire  damage. 

How  satisfactory  this  prediction  is,  depends  on  the  requirements 
of  the  model  user  and  on  how  widely  these  limits  are  separated. 

The  third  approach  is  essentially  to  predict  the  mean  weather 
for  any  future  year  and  calculating  the  fire  effect  accordingly. 
This  is  hot  rigorously  valid  because  the  use  of  mean  weather  will 
yield  an  exact  mean  of  fire  effect  only  if  the  latter  is  linearly 
related  to  all  weather  factors  (wind,  relative  humidity,  precipi¬ 
tation,  etc.).  Nevertheless,  the  method  seems  to  be  satisfactory 
on  two  counts:  First,  analysis  of  observed  data  on  past  fires 
indicates  that  relations  between  weather  and  certain  character¬ 
istics  of  fire  spread  are  not  too  far  from  linearity.  These 
relations  are  between  rate  of  spread  on  the  one  hand  and  wind 
velocity  and  relative  humidity  on  the  other.  Second,  the  non¬ 
linear  relation  between  weather  and  the  rural  fires ’  capability 
of  spread  will  receive  a  special  treatment,  as  will  be  discussed 
below* 
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In  attempting  to  determine  a  mean  weather  pattern  for  fire 
spread  prediction,  several  points  must  be  kept  in  mind: 

•  Weather  cannot  be  represented  by  a  single  number.  As 
many  numbers  are  required  as  there  are  meteorological 
phenomena  which  characterize  weather.  Therefore  it  is 
not  a  simple  matter  to  define  a  weather  pattern  that  is 
typical  of  the  area. 

•  The  statistical  study  of  past  weather  for  the  purpose  of 
fire  spread  prediction  must  center  on  those  phenomena 
that  can  significantly  affect  fire  spread.  Records  of 
past  fires  have  shown  that  the  general  fire  behavior  is 
most  sensitive  to  wind  velocity  and  direction,  relative 
humidity,  and  precipitation. 

•  For  a  given  area,  some  correlation  may  exist  between  var¬ 
ious  weather  factors.  For  example  low  relative  humidity 
may  be  associated  with  wind  in  certain  directions.  This 
correlation  may  be  taken  into  account  by  studying  the 
seasons  separately  and  possibly  by  distinguishing  day 
weather  and  night  weather  in  cases  where  the  weather 
seems  to  follow  characteristic  patterns  for  days  and 
nights . 

o  The  variability  of  weather  can  play  an  important  role. 

It  has  been  mentioned  that  for  rural  fires  weather  condi¬ 
tions  can  be  divided  into  three  groups.  One  group  makes 
the  fire  spread,  another  group  prevents  it  from  spread¬ 
ing,  and  the  third  group  extinguishes  it.  Variations 
within  the  same  group  would  produce  no  significant  effect 
in  fire  spread.  However  the  more  often  weather  conditions 
pass  into  the  no-spread  or  the  extinction  group,  the  less 
fire  damage  should  be  expected  for  a  random  attack. 


All  these  factors  make  it  difficult  to  define  and  calculate 
a  mean  weather  pattern,  i.e.,  a  mean  combination  of  all  weather 
factors  and  their  changes.  An  approximate  mean  weather  pattern 
can  be  obtained,  however ,  by  selecting  from  the  weather  record  a 
single  year  in  which  wind  velocity,  wind  direction,  relative 
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humidity,  precipitation  and  the  variations  of  each  appear  to  be 
most  typical.  The  actual  weather  in  that  year  is  then  used  to 
calculate  the  mean  fire  effect.  Selection  of  the  most  typical 
year  can  be  done  in  two  ways.  One  is  qualitative  and  simple, 
the  other  somewhat  more  quantitative  and  more  complicated. 

The  first  method  of  selection  requires  plotting  each  day's 
wind  velocity  (for  each  predominant  direction)  relative  humidity, 
and  precipitation  on  separate  graphs  and,  by  visual  inspection, 
selecting  a  year  which  seems  most  typical  with  respect  to  all 
four  weather  factors  and  their  variations.  To  include  all  types 
of  weather  patterns  that  may  occur  in  the  area,  the  years  examined 
must  cover  one  average  long  trend  cycle.  In  many  cases  this  cycle 
varies  from  about  10  to  20  years.  This  simple  method  is  probably 
sufficient  in  an  area  where  the  weather  pattern  does  not  vary 
appreciably  from  year  to  year. 


The  second  method  aims  at  establishing  a  measure  of  "weather 
normality"  for  each  season  of  each  year.  Normal  weather  is  made 
up  of  average  wind  velocity,  average  relative  humidity,  average 
precipitation,  and  average  cycle  of  each  for  a  period  of  10  to 
20  years.  A  measure  of  weather  normality  for  each  season  of  each 
year  is  obtained  by  combining  the  deviations  of  season  averages 
from  normal  weather.  The  most  normal  season,  i.e.,  the  season 
with  minimum  combined  deviations,  is  selected. 


The  main  steps  in  the  calculation  are  as  follows: 

(1)  A  season  average  of  wind  velocity  (for  each  predominant 

H  h  Tr\n  ^  n  q  aK f  q  a  a^U.  «  n  « -C  ^  «  -t. 

**■*•*■  v*  *-'-'**/  -«-v  wuuai.auu  lui  cauu  ocaouu  UJ4  caCll  C  d  X  U GdX 

year.  This  average  is  denoted  by  WF(i),  Wg (i) ,  W  (i), 
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and  Wyj(i)  for  the  year  i  and  for  spring,  summer,  autumn, 
and  winter  respectively.  It  is  calculated  from  the 
formula 

WF(i)  =  I.  Wp(j )  for  spring  (70) 

i 

and  from  similar  formulas  for  other  seasons.  Here  WF(j) 
is  the  average  wind  velocity  for  the  day  j.. 

Similar  season  averages  are  also  obtained  for  relative 
humidity  and  precipitation.  They  are  denoted  by  Hp(i), 
Hg(i),  HH(i),  Hj^Ci)  for  relative  humidity  and  by  Pp-(i), 
Ps (i),  Py(i),  Pw(i)  for  precipitation.  Wind  velocity 
is  expressed  in  mph,  relative  humidity  in  percent,  rain¬ 
fall  in  inches,  and  snowfall  in  inches  of  equivalent 
rainfall. 

(2)  The  season  averages  are  again  averaged  over  the  entire 
period  of  n  (10  to  20)  years  according  to  the  formulas 

_  1  — 

W_  =  —  )  W„(i)  and  so  on  for  other  seasons 
r  n  .  r 

i 

H„  =  —  T  H„(i)  and  so  on  for  other  seasons  (71) 
r  n  V  r 

L 

P_  =  —  )  P„(i)  and  so  on  for  other  seasons 
F  n  'r  Fs 

L 

(3)  Deviations  of  seasonal  averages  from  the  period  average 
are  then  computed  from  the  formulas 


%<*)  =  v1) 

-WF 

and 

so 

on 

for 

other 

seasons 

DHF(i)  =  5Fa) 

-«F 

and 

so 

on 

for 

other 

seasons 

(72) 

DPF(i)  =  Vi} 

-?F 

and 

so 

on 

for 

other 

seasons 
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In  these  equations,  the  first  subscript  on  D  indicates 
wind,  relative  humidity,  or  precipitation  and  the  second 
subscript  indicates  the  seasons. 


(4)  The  variability  of  each  weather  factor  is  measured  by 

comparing  the  average  cycle  (short  range)  for  each  season 
of  each  year  to  the  average  cycle  of  the  period  of  n  years. 


The  average  cycle  for  each  season  is  obtained  by  dividing 
90  by  the  number  of  cycles  counted.  One  cycle  is  the 
number  of  days  during  which  wind,  relative  humidity,  or 
precipitation  makes  a  complete  uptrend  and  a  complete 
downtrend,  with  minor  fluctuations  lasting  one  or  two 
days  ignored. 


Seasonal  average  cycles  are  denoted  by  Wj,(i),  hp(i), 
Pp(i)  for  wind,  relative  humidity,  precipitation  and 
are  given  by 


WF(1) 


90 

°W 


yt> 


90 

CH 


90 


(73) 


where  cy,  cp,  cp  are  the  numbers  of  cycles  in  wind, 
relative  humidity,  and  precipitation  actually  counted. 


Seasonal  average  cycles  are  then  averaged  over  the 
period  of  n  year: 

w_  =  —  Y  w_(i) 

F  n  Lr  P-  ' 

l 

Kf  =  |  £  hp(i)  (74) 

i 

P F  “  n  ? 
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Deviations  of  each  seasonal  average  from  the  grand 
average  for  the  period  are  computed  from 

DwF(i)  =  «f<1)  '  "F 

DhF(1)  =  SF(i)  -  ¥F  (75) 

DpF(1)  =  PF(i)  "  PF 

(5)  The  extent  to  which  the  overall  weather  of  each  season  in 
the  year  i  deviates  from  normal  conditions  is  expressed 
by  the  numbers  dp(i),  dg(i),  d^(i),  and  d^i)  defined  by 

dF(i)  “  DWFCi)  +  aDHF(i)  +  BDPF(i)  +  YDwF(i)  +  6DhF(i) 

+  F(i)  (76) 

where  a,  3,  y>  6,  e  are  normalizing  factors  such  that 

Wp  =  aHp  =  8Pp  =  ywf  =  6hp  =  spp  (77) 

These  normalizing  factors  imply  that  wind  velocity, 
relative  humidity,,  precipitation,  and  their  variabilities 
are  given  equal  importance  in  making  up  the  overall 
weather. 

The  four  seasons  having  minimum  d(i)  but  not  necessarily 
belonging  to  the  same  calendar  year,  will  be  selected 
and  the  actual  weather  patterns  in  these  seasons  will  be 
used  to  calculate  fire  spread. 


USE  OF  WEATHER  DATA 


When  a  specific  date  in  the  past  or  the  current  date  is  speci¬ 
fied  as  attack  date,  it  is  simply  required  to  evaluate  the  fire 
effect  due  to  the  particular  weather  pattern  following  the 
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specified  date.  With  an  unspecified  attack  date,  all  days  of  a 
future  period  of  many  years  are  equally  Ikkely  to  be  attack  days, 
and  the  problem  is  not  to  find  the  fire  effect  for  one  particular 
attack  day,  but  to  find  the  mean  fire  effect  for  all  possible 
attack  days.  By  finding  a  past  year  that  best  represents  the 
normal  weather  of  the  area,  we  have  avoided  the  immense  problem 
of  calculating  the  fire  effect  for  each  and  every  day  of  the  entire 
period  and  have  equivalently  changed  the  unspecified  attack  date 
problem  into  one  which  specifies  the  year  and  season  of  attack. 

It  is  now  required  to  evaluate  a  mean  fire  effect,  given  that  the 
attack  occurs  on  each  day  of  that  year  with  equal  probability. 

This  can  be  accomplished  by  two  simple  methods,  one  for  urban  fuel 
and  the  other  for  rural  fuel. 

1.  Since  the  spread/no- spread  capability  of  urban  fires  is 
essentially  weather- independent ,  the  day  of  attack  and 
its  weather  affect  the  fire  spread  only  to  the  extent 
that  the  rate  of  spread  depends  on  wind  velocity  and 
possibly  on  precipitation.  Therefore  the  attack  can 
occur  on  any  day,  and  the  mean  fire  effect  is  obtained 

by  a  single  fire  spread  calculation,  provided  average  wind 
velocity  and  average  precipitation  are  used  to  determine 
the  rate  of  spread. 

2.  For  rural  ^uel,  the  weather  of  che  attack  day  determines 
whether  the  fire  will  spread  or  not  and  how  long  it  will 
spread.  Thus  the  ultimate  result  depends  strongly  on 
the  attack  day.  However  the  fire  spread  calculation  can 
also  be  independent  of  the  attack  days  if  a  mean  spread 
duration  is  calculated.  The  mean  spread  duration  is  the 
mean  number  of  consecutive  days  of  spread  condition.  For 
the  same  weather  pattern,  this  mean  duration  varies  with 
each  type  of  rural  fuel  (brush,  grass,  timber).  The 
general  procedure  for  rural  fuel  is  to  allow  the  fire  to 
spread  for  a  number  of  days  equal  to  the  mean  spread  dur¬ 
ation  at  a  rate  determined  by  the  average  wind  velocity 
and  average  relative  humidity. 
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For  both  urban  and  rural  fires,  a  calculation  is  made  for 
each  of  the  four  seasons.  The  results  may  or  may  not  be  combined, 
depending  on  each  particular  case. 

PROBLEM  OF  WIND  DIRECTION 

In  a  given  area  and  during  a  given  season,  the  wind  is  likely 
to  blow  in  several  different  directions.  Averaging  wind  direc¬ 
tions  may  in  some  cases  lead  to  meaningless  or  erroneous  results. 
Suppose  the  wind  is  either  eastward  or  southward,  and  each  direc¬ 
tion  is  generally  stable  for  several  days.  Then  the  area  extending 
eastward  or  the  one  extending  southward  from  the  initial  fire  will 
burn,  depending  on  the  attack  day.  Averaging  the  two  wind  direc¬ 
tions  would  lead  to  the  result  that  the  area  extending  southeastward 
burns,  which  it  never  does  (see  Fig.  15A), 

Suppose  however  that  the  wind  in  this  case  changes  direction 
frequently,  say,  every  few  hours.  The  spread  pattern  can  be 
pictured  as  in  Fig.  15B.  It  is  obvious  that  the  use  of  an  average 
wind  direction  would  lead  to  approximately  the  same  result. 

These  considerations  suggest  some  practical  solutions  to  the 
problem  of  wind  direction.  In  studying  weather  data,  some  atten¬ 
tion  must  be  given  to  wind  directions,  their  distribution,  and  the 
frequency  of  their  changes.  If  the  changes  in  direction  are  small 
(less  than  90  deg)  or  if  the  wind  shifts  frequently  between  two 
widely  divergent  directions,  the  average  direction  is  recorded 
and  used  in  the  calculation. 

If  Che  wind  changes  periodically  (every  week  or  more)  between 
two  or  more  predominant  directions ,  the  total  time  spent  in  each 


13-12 


URS  641-6 


direction  must  be  recorded.  This  time  is  used  to  evaluate  the 
probability  that  fcr  random  attack  the  fire  will  spread  in  each 
direction.  For  the  case  of  Fig.  15A,  if  the  wind  blows  eastward 
for  a  total  of  60  days  and  southward  for  a  total  of  30  days,  the 
final  result  is  that  the  eastward  area  will  burn  with  a  probability 
0.66,  and  the  southward  area  will  burn  with  a  probability  0.33. 


The  same  approach  can  be  extended  to  the  case  in  which  the 
wind  changes  between  two  opposite  directions.  However  if  these 
changes  are  frequent,  two  different  situations  must  be  considered. 

(a)  The  wind  blows  in  each  opposite  direction  for  an  average 
time  T^  less  than  the  burning  time  of  the  fuel.  The 
fire  will  spread  in  both  directions  with  an  overall  rate 
approximately  equal  to  half  the  rate  which  corresponds 
to  the  actual  wind  velocity. 

(b)  The  wind  blows  in  each  direction  for  an  average  time  T2 
greater  than  the  burning  time  of  the  fuel.  In  these 
cases  there  will  be  no  significant  spread  in  either 
direction. 


In  both  (a)  and  (b) ,  side  spread,  i.e. ,  at  right  angle  to 
the  wind  directions,  must  be  examined.  The  extent  of  this  spread 
depends  on  wind  velocity  and  on  urban  and  rural  fuel  types. 
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Section  14 

LOCAL  APPLICATION  OF  FIRE  SPREAD  MODELS 

Theoretical  models  have  been  formulated  in  a  general  manner, 
largely  independent  of  dimension  scale.  The  size  of  the  area 
which  must  be  surveyed  for  fuel  characteristics  and  the  effort 
required  for  data  processing  and  model  computation  did  not  have 
to  be  considered. 

Methods  for  large-scale  prediction  have  been  obtained  by 
implementing  these  theoretical  models  with  input  data  at  a 
level  of  detail  determined  by  the  available  resources  for  the 
three  tasks,  data  collection,  data  processing,  and  model  compu¬ 
tation.  The  level  of  detail  in  the  output  information  is  deter¬ 
mined  by  the  level  of  detail  in  the  input  data.  The  character¬ 
istics  of  the  output  information  for  large-scale  application  have 
been  discussed  at  some  length  under  output  requirements  (Section  3) . 

Application  of  fire  spread  models  at  the  local  scale  involves 
a  higher  level  of  detail  in  the  input  data  rather  than  a  change  in 
the  basic  models.  Following  are  general  procedures  for  data  col¬ 
lection  and  model  application  for  cases  in  which  the  fuel  area  is 
a  large  or  moderate-size  city. 

DATA  COLLECTION 

Cities  consist  mainly  of  urban  fuel.  We  do  not  however 
exclude  the  possibility  of  isolated  patches  of  rural  fuel  large 
enough  to  be  treated  as  such.  The  only  required  data  for  these 
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patches  Is  fuel  type  —  brush,  grass,  timber,  and  hardwood.  Major 
fuel  gaps  are  unlikely  in  an  area  of  this  size  and  the  effect  of 
fuel  density  could  not  be  determined  from  observed  data  currently 
available. 


For  urban  fuel,  the  following  data  will  be  needed:  type, 
fuel  density,  mean  roof  area,  or  spacing  distribution. 

Urban  fuel  types  of  importance  at  this  scale  are  light 
residential,  heavy  residential,  commercial,  city  center,  and 
industrial.  Quantitative  definitions  of  the  five  types  must  be 
given  before  an  actual  survey  can  be  carried  out  successfully. 

What  is  city  center  fuel  in  one  city  may  be  the  same  as  commer¬ 
cial  fuel  in  another  city.  Each  fuel  type  is  characterized  by 
a  number  of  factors  important  in  fire  spread  but  which  cannot  be 
investigated  separately.  These  factors  are:  structural  materials, 
fuel  loading  or  amount  of  fuel  in  each  unit,  distribution  of 
number  of  stories,  types  of  outside  walls ,  and  percentage  of 
window  area. 


Other  factors  such  as  mean  roof  area,  fuel  density,  spacing 
distribution  may  also  be  related  to  fuel  type.  They  may  however 
vary  considerably  within  each  type.  This  possibility  must  be 
examined  by  selecting  a  sample  of  each  type  and  measuring  these 
factors.  If  variations  within  each  type  are  sufficiently  small, 
the  results  of  these  measurements  are  permanent  data  and  can  be 
used  for  any  area  with  known  fuel  type. 


Fuel  density  is  the  ratio  of  total  plan  area  to  the  total 
area,  which  includes  normally  occurring  empty  spaces,  such  as 
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streets,  alleys,  lawns,  yards,  and  parking  areas  but  not  partic¬ 
ularly  large  open  spaces  that  are  specified  separately. 

The  mean  rocf  area  can  be  determined  from  fuel  density  if 
the  total  number  of  structures  is  also  counted.  It  is  needed 
for  calculating  spread/no- spread  requirements  in  urban  fuel. 

(See  Section  4.) 

Spacing  distribution  refers  to  adjacent  structures  only.  It 
is  expressed  as  the  probability  that  an  adjacent  structure  is  in 
each  range  of  distance  from  any  structure  in  the  area.  A  struc¬ 
ture  is  considered  as  adjacent  to  another  if  at  least  half  of  it 
can  be  seen  from  the  second  structure. 

The  first  step  in  the  survey  work  is  to  identify  sections 
that  can  be  labeled  as  one  of  the  five  fuel  types  or  as  empty 
space.  Large  empty  spaces  that  might  occur  in  each  section  must 
also  be  specified. 

In  the  second  step,  fuel  density,  mean  roof  area,  and  spac¬ 
ing  distribution  are  determined  for  each  section.  Convenient 
methods  for  these  measurements  are  still  to  be  worked  out. 

Sanborn  maps  or  aerial  photos  are  likely  to  be  needed. 

Finally  the  entire  ai*ea  is  divided  into  grid  squares  100  to 
300  ft  on  each  side,  depending  on  the  particular  case.  These 
squares  are  called  cells  and  will  be  used  to  map  the  course  of 
the  fire  in  time  and  space. 
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MODEL  APPLICATION 

For  rural  sections,  the  capability  of  spread  depends  on 
previous  and  current  weather  and  must  be  predicted  by  referring 
to  tables  of  spread/no- spread  weather  conditions  for  each  type 
of  fuel.  Procedures  for  following  the  fire  from  cell  to  cell 
are  the  same  as  those  given  previously  for  large-scale  applica¬ 
tions.  However,  the  cell  size  and  time  increment  can  be  reduced 
by  factors  appropriate  to  the  extent  of  fuel  area  and  to  the 
level  of  detail  in  the  input  data. 

For  urban  sections,  spread  capability  is  predicted  by  using 
calculation  procedures  described  previously  in  connection  with 
model  IIC.  The  course  of  the  fire  can  be  established  in  the 
same  manner  as  for  rural  sections. 
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Section  15 

CONCLUSIONS  AND  SUGGESTIONS 

A  semiempirical  method  for  fire  spread  prediction  is  now 
available.  With  the  currently  existing  data  the  method  gives  an 
approximate  prediction  of  fire  spread  in  rural  and  urban  fuels. 

The  accuracy  of  prediction,  however,  can  be  significantly  improved 
by  applying  the  method  according  to  a  special  scheme  and  by  acquir¬ 
ing  more  empirical  data,  as  we  shall  now  suggest. 

SPECIAL  SCHEME  FOR  MODEL  APPLICATION 

The  process  of  collecting  fuel,  topography,  and  weather  data 
generally  constitutes  a  major  task.  It  has  been  suggested  that 
fuel chta  be  determined  permanently  for  each  specific  type  of 
fuel,  for  example,  light  residential  or  heavy  residential,  by 
surveying  some  area  samples  containing  only  one  type.  Data  spec¬ 
ification  for  each  fire  is  done  simply  by  giving  the  proportion 
of  each  type  in  the  total  area  and  certain  fuel  characteristics 
that  cannot  be  regarded  as  inherent  in  each  type. 

In  general  the  level  of  detail  of  the  output  information 
follows  that  of  the  input  data,  which  —  in  turn  —  depends  on  the 
size  of  the  total  area.  Therefore  the  area  surveyed  should  not 

t 

exceed  greatly  the  extent  of  expected  fire  spread.  It  seems 
good  economy  to  carry  out  data  collection  and  fire  spread  calcu¬ 
lation  simultaneously. 
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For  rural  areas,  a  detailed  survey  is  first  made  for  a  dis¬ 
tance  (extending  from  the  initial  fire)  that  can  be  expected  to 
be  covered  by  fire  before  a  weather  change.  Fire  spread  calcula¬ 
tion  is  made  for  the  section  thus  surveyed.  A  second  section  is 
surveyed  only  if  fire  is  expected  to  reach  it  during  the  next 
weather  period.  Similarly,  for  urban  areas,  the  section  surveyed 
each  time  is  determined  by  the  answers  to  the  questions:  Can  fire 
spread  through  it?  What  distances  can  the  fire  travel  at  the 
speed  predicted  by  the  current  weather?  By  means  of  this  scheme, 
it  is  possible  to  keep  the  surveyed  area  at  a  minimum. 

ACQUISITION  OF  MDRE  EMPIRICAL  DATA 


A  semiempirical  model  is  as  good  as  its  parametric  data.  In 
Part  II  of  this  report,  existing  data  on  rate  of  spread  in  rural 
and  urban  fuels  have  been  analyzed,  and  some  additional  data  on 
weather  conditions  for  spread  or  no-spread,  for  extinction  in 
rural  fuel,  and  on  spread  probability  in  urban  fuel  have  been 
estimated.  A  quick  look  at  the  contents  of  Part  II  would  reveal 
the  need  for  more  empirical  data  not  only  for  better  implementa¬ 
tion  of  this  method  but  for  general  understanding  of  fire  spread 
also.  The  more  important  ones  are  briefly  indicated  below: 

•  Additional  Investigation  on  Rate  of  Spread  Data  in  Rural 
Fuel .  A  fair  amount  of  data  on  rate  of  spread  has  been 
obtained.  For  ground  spread  the  dependence  of  rate  on 
wind  velocity  and  direction  and  on  relative  humidity  has 
been  established  in  some  detail.  The  effect  of  long-range 
spotting,  however,  has  not  been  evaluated.  This  would 
require  further  investigation  of  the  present  data.  This 
investigation  must  probably  be  backed  up  by  a  small  amount 
of  new  data. 
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•  Additional  Data  on  Rate  of  Spread  in  Urban  Fuel,  Effects 
of  wind  and  humidity  (exterior)  on  rate  of  spread  cannot 
be  reliably  evaluated  from  present  data. 

•  Burning  Time  of  Rural  and  Urban  Fuels.  Some  estimates  for 
rural  fuel  have  been  given  in  Part  II.  More  detailed 
information  from  direct  measurements  on  both  rural  and 
urban  fires  is  needed. 

•  Spread  Probability  in  Urban  Fuel.  This  is  the  probability 
that  for  a  given  weather  condition  and  a  given  type  of 
urban  area,  fire  will  spread  from  one  structure  to  another 
at  a  known  distance.  It  is  required  in  Models  II  and  IIA 
for  derivation  of  the  parameter  A/|i  ,  in  Model  I1C  for 
derivation  of  spread/no-spread  conditions,  and  in  methods 
for  calculating  fire  spread  vulnerability.  This  proba¬ 
bility  can  be  obtained  in  two  ways: 

1.  The  probability  of  no-spread  can  be  studied  by  examining 
maps  of  burnout  areas  in  World  War  II  incendiary-attacked 
cities  or  in  any  large-scale  urban  fire.  For  a  given 
type  of  fuel  area,  the  pairs  of  burned  and  unburned  struc¬ 
tures  (usually  near  the  burnout  contour)  with  a  given 
spacing  are  counted  and  tabulated.  It  can  be  shown  that 
the  distribution  of  spacxngs  (a  characteristic  of  the 
area,  not  of  the  fire)  must  be  combined  with  these  data 

to  give  the  probability  of  no-spread. 

2.  Alternatively  and  perhaps  more  basically,  the  maximum 

spacing  required  for  ignition  to  occur  is  determined  for 
each  narrow  range  of  responsible  factors.  These  factors 
are:  type  of  structural  materials,  fire  load,  and  size. 

Data  on  maximum  spacings  can  be  used  as  such  in  certain 
applications,  for  example,  in  determining  the  fire  spread 
vulnerability  of  resources.  They  can  also  be  combined  with 
the  distributions  of  type,  fire  load,  and  size  in  a 
stochastic  fuel  area  to  obtain  the  probability  of  spread 

as  a  function  of  distance.  Data  on  maximum  spacings  can 
be  obtained  from  maps  of  burnout  areas  if  documentation 
is  sufficiently  detailed  or  from  careful  observation  of 
real  or  experimental  fires. 
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Appendix  A 
FIRE  FRONT  MODELS 

BASIC  FIRE  FRONT  MODEL  (MODEL  I) 

The  mechanisms  of  the  model  and  the  definition  of  terms  has 
been  given  in  Section  6,  Part  I. 

A  simple  argument  on  the  relation  between  adjacent  proba¬ 
bility  states  leads  to  the  equation 

Pn  (t  +  dt)  =  Pn_j  (t)  Xdt  +  Pn  (t)  [l  -  (X  +Mf)  dt]  (78) 

which  reduced  to 

-  *P-l<t>  -  «P.(t)  «  =  *■  *  M, 

For  n  =  0,  this  equation  becomes 

=  -  ep.(t)  (80) 

with  solution 

Po(t)  =  e“Ut  (81) 

For  n  0,  we  use  the  Laplace  transform  Pn(s)  for  Pn(t)  and  write 

SPr<  (S)  -  Pn  (0)  -  XPn_,  (S)  -  01^  (S) 

Since  for  n  £  0,  P^(0)  -  0,  we  get  the  solution 


(82) 
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Ms)  = 


(s  +  0) 


n  +  i 


(83) 


or 


r„  (t)  ^  ^ 

The  expression  for  Qn ( t )  can  be  easily  obtained  either  by 
direct  integration  of  P  (t)  or  through  Laplace  transformation. 

The  probability  that  at  time  t  the  fire  front  is  still  alive 
in  some  cell  is  given  by 


P0  (t)  +  (t)  +  P2  (t)  +  .  .  .  ad  inf. 


or 


o ,-ot 


1  +  xt  +  am  +  imL  + 

2!  3! 


(84) 


-Pft 

Since  this  probability  is  equal  to  e  ,  the  meun  lifetime 
t  of  the  fire  is  given  by 


-u.  t 


d  e 


-Pi 


dt 


dt  =  — 


1_ 

Mr 


(85) 


Thus  the  mean  ultimate  burnout  distance  n  (°°) ,  found  to  be 
is  the  product  of  the  mean  velocity  and  mean  lifetime. 


APPLICATION  TO  INHOMOGENEOUS  CELLS  (MODEL  IA) 


By  Laplace  transform,  the  equation 


0£  ^n-1,11  In— 1  (Q  “  (^n,n+l  +  O') 


(86) 
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becomes 


S  J'Jn  (s)  ^-n-l,n  Pn-1  (®)  ~  (^-n,n+l  "**  Mn)  Pn  (®) 


(87) 


P„  <S)  «  — ,  Pn -i  (S) 

S  +  7-n.ntl  +  Pn 


(88) 


whose  solution  is 


- 1,  i 


n'M  - 1, 

rrx“" 


i  +  ."i 


(89) 


and  since 


P0(3> 


s  +  xoi  +  MD 


(90) 


we  have 


Pn(s)  = 


i — n 

-01  +  Mol  1 


^  i-l,i 


s  +  X 01  +  Mo  1  I  s  +  X  i ,  i  +.1  +  Mi 


n  __ 

5  (S  *  Xi,  i+1  +  Mi) 


(91) 


The  Laplace  transform  of  Q^(t)  i-s 


9n  (s) 


_  A«„  P„  (s) 


(92) 


Therefore 


nx,.,  i 

Qn  (s)  =  ^  - 1 - - 


n(s  +  Ai>i+1  +  Mi) 


(93) 


Terms  in  Q  (t)  corresponding  to  factors  0  (s  +  \  .  +  U,. )  involve 

n  r ,  j.  ■  j. 

e_ct:,  where  c  is  a  constant,  and  therefore  vanish  at  t  =  ^ 
Retaining  only  nonvanishing  terms,  we  have 
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Qn{°°)  = 


n  +  Mi) 


^  n-1 

\i,n+i  +  Mn  i=0  ^i,i+l  + 


(94) 


APPLICATION  TO  SPOTTING  (MODEL  IB) 


Starting  from  the  basic  equation 


=  XiP„-i(t>  +  x2p„.2(t)  +  . . .  +  x„P0(t) 


-  (n,  +  x,  +  x2  +  . .  ,)P„(t) 


we  obtain  by  Laplace  transformation 


SPn(s)  -  pn(0)  -  Xj  Pn_!(s)  +  A2Pn_2(s)  +  .  .  .  +  X„  po(s) 


where 


-  I  Mr  +  ^  Xjj  Pn(s) 


(95) 


(96) 


P„(0)  = 


n  =  0 
n  ^  0 


For  n  =  0 


-(nf  +  z\)t  i 

oW  =  e  \  y  and  p0(s)  =  ;  — 


s  +  nt  +  SXj 


(97) 


For  n  0 


„  _  Xl  R.-l(s)  +  x2  P n— 2  (s)  +  .  .  .  +  Xn  pn(s) 

A(S)  8  +  Ht  +  2X, - 


(98) 
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Pn+i(s)  = 


_  xi  Pn  (S)  +  x2  pn_i  (s)  +  .  .  ,  +  Xn  +  t  p0  (s) 


s  +  ptf  + 


Defining  the  generating  function  p(z)  and  X(z)  as 


p(z)  =  p0(s)  +  Pi(s)z  +  p2(s)z2  +  ...  +pn(s)zn 


(99) 


X(z)  =  0  +  Xjz  +  X2z2  +  .  .  .  +  Xnzn 


(100) 


Eq ,  (98)  can  be  written 


p(z)  = 


s  +  +  2X{  -  X(z) 


(101) 


Since 


Qn  (t)  -  /»,  f  PnWM 


(102) 


we  can  write 


9n(s)  = 


Mi  Pn  (s) 


(103) 


9n(s)  = 


Al  9n_l(s)  +  ^qn-2(s)  +  .  .  .  +  Xnq0(S) 
s  +  +  SXj 


(104) 


q0(s)  = 


~  where 
s  o 


S  =  s  +  +  2x5 


Mt 

q,<s> 


Mi  M  .  Mi^'i 

q2(s)  7§r  ~75<~ 


(105) 


A-6 


URS  641-6 


Factors  in  S  will  lead  to  terms  in  Q^(t)  involving  negative 
exponentials  in  t,  hence  vanishing  at  t  =  », 


Qo(°°) 


_ a±_ 

Mf  + 


(106) 


Qi  H 


Mf*i 
(Mf  + 


^•1  Qo(°o) 

Mf  +  2^i 


(107) 


Q2  (°°) 


Mf  A2 

(Mf  +  2Ai)2  + 


Mj  Xi 

(Mf  +  S^i)3 


+  A2Q0(qq) 

Mi  + 


(108) 


Qn  H 


^lQn— 1(°°)  +  ^Qn-2(°°)  +  ♦  •  •  jMjjoH 

Mf  + 


(109) 


(°°) 


A 

It  can  be  easily  seen  that  the  generating  function  Q(z)  of 
is  given  by 


Q  (z) 


Mr 

Mr  +  -  X  (z) 


(HO) 


although  this  relation  is  probably  not  as  useful  as  the  relation 

for  Q  («>)  . 
n 

SPECIAL  CALCULATION  OF  MEAN  LIFETIME  OF  FIRES 


To  obtain  the  parameter  we  must  know  the  mean  lifetime 
of  the  fire,  i.e.,  we  must  observe  the  lifetimes  or  the  spread 
distances  and  velocities  of  many  actual  fires  burning  under  the 
same  or  similar  fuel  and  weather  conditions.  Such  information 
does  not  exist  in  records  of  past  fires  and  is  not  likely  to 
become  available  in  the  future.  In  one  special  case,  however, 
the  mean  lifetime  can  be  calculated  directly,  as  we  shall  now 
discuss . 
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Suppose  the  weather  is  steady  and  corresponds  to  spread  condi¬ 
tion  for  the  type  of  fuel  under  consideration.  The  fire  spreads 
until  it  encounters  a  fuel  break  that  it  cannot  jump.  A  fuel 
break  here  is  a  cell  void  of  fuel  or  containing  protected  fuel. 

We  call  such  a  cell  empty  and  any  other  cell  full.  The  empty 
cells  are  randomly  distributed.  Let  e  be  the  fraction  of  empty 
cells  and  1  -  e  the  fraction  of  full  cells  or  "fuel  density"  of 
the  area. 

There  are  two  ways  of  specifying  the  fuel  density:  (1)  the 
fraction  e  of  empty  cells  is  given  for  an  unlimited  distance 
(or  area).  A  consequence  of  this  specification  is  that  the 
probability  of  having  a  continuous  series  of  full  cells  is  finite 
even  for  very  long  series;  (2)  the  fraction  e  has  been  measured 
over  a  given  finite  distance.  In  this  case  the  largest  possible 
spread  distance  is  the  total  number  of  cells  minus  the  number  of 
empty  cells.  The  probability  of  spreading  to  a  larger  distance 
is  zero. 


We  shall  now  calculate  the  mean  spread  distance  in  each  case 
and  obtain  the  mean  lifetime  by  dividing  the  mean  distance  by 
the  mean  velocity. 


(1)  With  the  exception  of  cell  0,  which  has  been  ignited 
initially,  all  other  cells  have  the  probability  e  of 
being  empty  and  1  -  e  of  being  full.  Thus  the  proba¬ 
bility  of  spreading  to  zero  cell  is  e  and  the  probabil¬ 
ity  of  spreading  to  one  or  more  cells  is  1  -  e  (here 
denoted  by  f ) .  The  probability  of  spreading  to  one  cell 
only  is  fe  and  so  on.  The  mean  spread  distance  D  in 
terms  of  cell  width  is 


»m  = 


fe  ±  2  fe  +  3  fe  +  ... 
e  +  fe  +  fe  +  ... 


1  -  e 


e 


(111) 


(2)  In  a  sequence  of  N  cells,  v  cells  are  empty  and  are 
randomly  distributed,  with  v/N  =  e, 
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The  number  of  different  ways  a  group  of  v  cells  can  be 
selected  from  N  cells  is 


N\  _  N! 
v  )  v  !  (N  -  v) ! 


(112) 


For  the  fire  to  spread  to  cell  n  and  then  stop,  cells 
1,  2,  .  .  n  must  be  full  and  cell  n  +  1  must  be  the 

first  cell  to  be  empty.  The  remaining  v  -  1  empty  cells 
are  distributed  randomly  among  a  total  of  N  -  n  -  1 
cells.  Out  of  N  -  n  -  1  cells,  a  group  of  v  ■  1  cells 
can  be  selected  in 


N  -  n  -  l\  (N  -  n  -  1) !,  m3) 

v  -  1  /  (u  -  1) !  (N  -  n  -  v) .' 


different  ways  with  the  obvious  condition  n  £  N  -  v. 

Therefore  the  probability  of  spreading  to  n  and  only 
n  cells  is 


(N  -  n  -  1) !  v  !  (N  -_v)l 
(v  -  1) !  (N  -  n  -  v)  l  N ! 


(N  -  n  -  -1)  !(N  -  u)  lv 
(N  -  n  -  v) !  N ! 


(114) 


This  probability  can  also  be  expressed  as 


F(F  -  1)(F  -  2)  .  .  .  (F  -  n  +  1) 
(N  -  1)(N  -  2)  ...  (N  -  n) 


(115) 


where  F  represents  N  -  v,  the  total  number  of  full 
cells . 


fi 
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It  is  easily  seen  that  for  a  given  v/N,  if  N  and  F 
tend  to  infinity,  we  obtain  the  same  expression  as 
before,  i.e.,  fne. 

The  probability  distribution  for  n  =  0,  1,  2,  .  .  . 
can  be  found  from  this  formula  and  can  be  used  to 
calculate  the  mean  spread  distance  Dm. 
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Appendix  B 
FUEL  STATE  MODELS 


BASIC  FUEL  STATE  MODEL  (MODEL  II) 

The  mechanism  of  the  model  and  the  definitions  of  U..(t), 

ij 

F^j(t),  B^.(t)  have  been  given  in  Section  6. 

The  probability  that  an  unignited  cell  ij  is  ignited  by  at 
least  one  of  its  eight  neighbors  during  time  dt  is 

Aij  dt  “  (A,  +  A2  +  .  .  .  +  Ag)  dt  (116) 


The  three  equations  of  the  model  are 


=  -  U„(t)  z  r.(t)A. 


-  U.j(t)  Z  r.(t)  A„  F,j 


(117) 

(118) 


Ms  Fu(t) 


(119) 


where  ij  is  any  cell  outside  the  initial  fire  and  the  u’s  indicate 
the  eight  cells  surrounding  i j . 

An  expression  for  B^ (®)  can  be  obtained  as  follows: 

Define  the  quantities  F^  (t)  and  B^  (t)  as 


F,*  (t)  =  2F„  (t)  A„ 


(120) 
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B,*  (t)  =  2Bu  (t)  Au 


(121) 


Equation  (119)  can  be  written 


9BU  (fc) 
8 1 


=  Ms  Fu  (t) 


(122) 


Multiplying  both  sides  by  and  summing  we  get 


9B*  (t) 

at 


=  Ms  F|*j  (t) 


(123) 


and  division  of  Eq.  (117)  by  Eq.  (123)  gives 


9Ujj(t) 
9Bfi  (t) 


-  *  u- (t) 


(124) 


which  integrates  to 


\?| j  (t) 


1  * 

-  Bfjd) 


(125) 


V* 


B^.  (0)  is  obviously  zero  and  since  cell  ij  is  not  part  of  the 

initial  fire,  U..(0)  is  1.  Therefore 

ij 


-IT  Bi'jW 

U,j(t)  =  e  s 


(126) 


By  use  of  Eq.  (119)  and  the  relation 


Uij(t)  +  Fu(t)  +  Bij  (t)  =  1 


(127) 


U..(t)  can  be  expressed  in  terms  of  B..(t),  and  Eq.  (126)  becomes 

J 
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M. s  0t 


=  1  -  Bjj(t)  -  e 


(128) 


As  t  tends  to  » ,  B„(t)  and  B^_.  (t)  approach  their  limits  which 
are  equal  to  or  less  than  1  and  b B„(t)/c>t  approaches  zero; 
therefore 


B,jH  =  1  -  e 


—  b  * 
-*1.  % 


or 


_7 

B|jH  =  1  -  e  s 


(129) 


LINEAR  STOCHASTIC  VERSION  (MODEL  IIA) 


Solution  for  Bn(°°) 


The  assumed  mechanism  can  be  described  by  the  equations 


(t) 

at 


=  -  Un  (t)  F..J  A 


(130) 


3F„  (t) 

■"TT"  «  Un  (t)  Fn_|  A  -  Fn  (t) 


(131) 


9Bn  (t) 

at 


=  MsFn(t) 


(132) 


with 


Fn{0)  = 


for  n  =  0 
for  n  ^  0 


By  changing  n  to  n  -  1  in  Eq.  (132)  and  multiplying  by  A  we 


have 


4 
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9Bn_t  (t)  A 


dt 


-  Ms  A  Fn_!  (t) 


(133) 


and  dividing  Eq.  (130)  by  Eq.  (133), 


aun  (t)  =  _1 

^WA  (i. 


Un  (t) 


(134) 


which  gives 


-aB  (0  a 

U„  (t)  =  e  a  -  ^ 


Upon  substitution  for  U^(t)  we  have  from  127  and  132 


(135) 


_1_  9Bn(t) 
Ms  at 


—  a  B  (t) 
n— 1 


=  1  "  Bn  (t)  -  e 


(136) 


and  for  t  = 


Bn  (°°)  =  1  -  e 


-a  b  («) 
n~l 


(137) 


Solution  for  F-|  (t)  and  B1  (t) 

Explicit  expressions  for  F^(t)  and  B^(t)  can  be  obtained  as 
follows : 


According  to  Eq.  (130), 


9Ui  (t) 
Ui(t) 


-n  t 

Ae  s  dt 


(138) 


t.tH  ioh 

«*  L  t  A.  ^  k  K 


nr  t  un  o 

wo 


vj  j  j.uuu.^iauiUil 


U,  (t)  =  e- 


->ae 


-“j 


(139) 


k 
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Substituting  for  ll^t)  in  Eq.  (131)  we  get  successively 

8F.  ft\  t  _ji  t 

— J-ii  =  Ae-a  eae  s  e  s  -  ms  Fi  (t) 


Fj  (t)  =  A  e-a  e 


-li  t  -P  t  _(!  t 

fae  s  /  eae  s  dt  +  ke  s 


(140) 


where  k  is  a  constant  of  integration.  Letting 


x  =  e 


and  dx  =  -  m  x  dt 


(141) 


and  integrating,  we  obtain 


Ft(t)  =  ae  a  e-^  Mst  +  a(l  -  e"^)  +  ~~  (l  -  e“  ^ 

*  rfe  •  •  ■] 


(142) 


which  can  be  integrated  term  by  term  to  give  B^(t)  according  to 
Eq.  (132). 


Mean  Lifetime  of  Individual  Cell 


Consider  the  probability  F^(t)  that  cell  n  is  burning  at 
time  t.  At  time  t  +  dt,  this  probability  is  given  by 


Fn  (t  +  dt)  =  Fn  (t)  (1  ~  dt) 


(143) 


d  Kn  (t)  •  -  Ms  Fn  (t)  d  t 


If  n  is  0,  i.e.,  the  initially  ignited  cell,  Fq(0)  =  1  and 


-/»  1 

F0(t)  *  e  •" 


(144) 
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For  any  other  cell  n,  let  9  denote  the  time  after  ignition,  then 

-  ii  o  ,  ,  . 

Fn(0)  =  e  s  (145) 

The  probability  that  cell  n  goes  out  between  0  and  9  +  d9  is 

F  (0)u  d9 .  The  mean  lifetime  r  is  therefore  given  by 
n  s  J 

/»  _u  0  1 

0  e  s  nsdQ  =  —  (146) 


Mean  Rate  of  Spread 

Once  ignited,  a  cell  has,  for  each  time  increment  dt,  a 
probability  Adt  of  igniting  the  next  cell  and  a  probability  p  dt 

b 

of  dying.  Therefore  the  time  for  jumping  from  one  cell  to  the 
next  is  a  stochastic  variable.  It  is  required  to  calculate  the 
mean  jump  time  given  that  the  fire  has  made  a  large  number  of 
jumps . 

Let  P(t)  be  the  probability  that  the  jump  takes  place  before 
time  t,  then 


P  (t  +  dt)  =  P(t)  +  [  1  -  P(t)j  Adt  (147) 

The  term  a  dt  does  not  appear  in  this  equation  because  it  is  known 
s 

that  none  of  the  cells  involved  died  before  the  jump.  The  solu¬ 
tion  is 


P(t)  =  1  -  e“Al  (148) 

The  probability  that  the  jump  occurs  between  t  and  t  +  dt  is 


dP(t)  =  Ae~At  dt 


(149) 
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Hence  the  mean  jump  time  t  j  is 

OO  4 

tj  =  J  t  Ae-At  dt  =  ^ 
0 

and  the  mean  rate  of  spread  R  is  given  by 


(150) 


(151) 


Final  Mean  Spread  Distance 

The  recursive  formula  for  burnout  probabilities  is 

-a  b  a 

Bn+i  =  1  -  e  n  a  =  -  (152) 

For  sufficiently  large  a,  approaches  a  limit  at  which 

B  =  B  , -  =  B  , 0  =  .  .  . .  The  limit  is  the  solution  of  the 
n  n+1  n+2 

equation 


x  -  1  -  e"ax 


or 


where 


a 


1  -  o' 


(153) 


x  -  ax 

The  solution  z  =  0  or  x  =  0  always  exists.  In  addition,  a  finite 
solution  exists  if 

1  „ 

~  ~  1  or  a  -  1 

cl 
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It  is  evident  that  when  the  limit  of  B  is  finite  the  sum  TIB  is 

n  -j.  n 

infinite.  We  now  show  that  when  the  limit  of  B  is  0  or 

n 

when  a  <  1,  the  sum  is  finite.  Equation  (152)  above  can  be 
written  as 


Bn+l 


aB„ 


azB;  + 

2 !  3 ! 


The  ratio  B  ,-,/B  is 
n+1  n 

B.+1  _  ^  a3Bz 

Bn  2:  '  3! 


(154) 


(155) 


Since  B  tends  to  zero  as  n  increases,  the  ratio  tends  to  a. 
Therefore  if  a  is  less  than  1,  the  ratio  is  also  less  than  1  and 
by  the  ratio  test  the  series  converges. 


Burning  Time  and  Ignition  Probability 


We  have  been  dealing  with  the  probability  Adt  that  ignition 
occurs  in  stochastic  fuel  during  the  time  increment  dt  given 
that  the  ignition  source  is  still  burning.  The  parameter  A  is 
also  called  probability  of  ignition  per  unit  time.  Evidently 
the  longer  the  fuel  burns,  the  greater  the  ultimate  probability 
of  ignition.  The  relation  between  the  stochastic  burning  time 
and  ignition  probability  is  derived  as  follows: 


Let  P(t)  be  the  probability  that  ignition  has  occurred  before 
time  t.  The  probability  that  ignition  occurs  between  t  and 
t  i  dt  is 


-  /»  1 

fiP(t)  =  o  s 


-  pm 

\_  / 


Ad  t 


1 


(156) 
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-nst 

where  the  factor  e  "  is  the  probability  that  the  ignition  source 
is  still  burning  at  time  t. 


Integration  from  0  to  t  gives 


(157) 


Since  stochastic  burning  time  may  extend  to  some  high  value,  the 
probability  that  ignition  ultimately  occurs  is 


P(co) 


1  -  c"a 


(158) 


It  is  this  probability  that  we  usually  derive  from  data  on  past 
fires . 


;  LINEAR  DETERMINISTIC  VERSION  (MODEL  IIB) 

I 

i 

The  basic  equation  representing  the  simplified  mechanism 

i 

i  of  the  model  is 


dX 

dt 


=  < 


k  W  (t) 

for 

W(t)  <  w0 

kW0 

for 

o 

Al 

J 

(159) 


where  X,  W(t),  W  have  been  defined  in  Section  6. 

The  width  of  the  burning  fire  is  the  distance  traveled  from 
the  time  a  fuel  element  is  ignited  to  the  time  that  same  element 
goes  out.  Therefore 


,  r  I  42?  Vw,  „  „  f‘ 


w  m  _  . 

J,.T  \  d  or 


w  tn\  /-l  n 

j  ^  ~ 

t  -  7 


n  £n\ 


m  * 
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or 

^2  =  k  [w(t)  -  W(t  -  T)]  (161) 

It  is  convenient  to  assume  that  W(0)  has  been  chosen  to  be  equal 

to  W  . 
o 

The  initial  fire  (assumed  to  be  ignited  simultaneously 

throughout  the  region  W(0))  will  keep  burning  until  t  =  t,  the 

burning  time  of  fuel.  At  that  time  it  suddenly  goes  out.  During 

the  period  0  -  t,  the  width  grows  at  the  constant  rate  WQk,  so 

that  at  time  t  =  t,  after  the  initial  fire  goes  out,  the  fire 

width  becomes  W  kT . 

o 

If  WQkr  is  equal  to  or  greater  than  Wq  (i.e„,  kT  £  1),  the 

velocity  remains  constant  at  W^k  and  from  Eq.  (160)  the  front 

width  remains  constant  at  W  kT . 

o 

For  the  case  W  kT  <  W  ,  the  value  of  W(t)  can  be  best  calcu- 
o  o 

lated  numerically  for  each  successive  time  increment  At. 

The  final  burnout  distance  is  readily  derived  in  the  follow¬ 
ing  manner. 

Divide  t  into  v  equal  intervals  of  length  d  =  t/v.  Let 
be  the  distance  covered  during  the  i^  time  interval.  From  the 
above  discussion 

Xj  =  X2  =  .  .  .  =  X  v  =  X0 

where  is  WQkT/v.  X  X^+£>  •  •  •  can  be  expressed  according 
to  Eq.  (161).  Thus 
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Xi  =  x0 

x2  =  x0 

•  •  • 

=  x0 

Xy+t  =  die  (Xy  +  Xy-1  +  .  .  „  +  Xt) 

Xy+2  ~  ^  ^  (^"y+l  +  Xy  +  .  .  .  +  X2) 

X  j  =  d  k  (X  i_t  +  X  j_2  +  .  .  .  +  Xt_y) 


(162) 


Summing  each  column  we  get 


2  Xj  =  vX0  +  dk 


2  Xj  +  2  Xj  + 

V  1  V- 1  1 


.  .  +  £  X 


1 

tA‘J 


(163) 


To  the  right  hand  side  we  add  and  subtract  the  quantity 


dk 


v-i  v-2 

2  Xj  +  2  Xt  + 


i 

.  +  2  Xi 
1 


(164) 


which  is  equal  to 


dkX0 


Equation  (163)  now  becomes 


2  X ,  -  vXr  -  d  k  X 
l 


i  -  i'  yv0  -  u -  +  dk  v  2  Xj 

^  1 


(165) 


Denoting  the  final  distance  EX^  by  D,  we  find  upon  substi¬ 


tution  of  r/v  for  d  and  W  kT /v  for  X 

o  o 


D  -  W°kr 
1  -  kr 


1  - 


kr  v-1 
2  v 


(166) 
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As  v  tends  to  infinity,  i.e.,  the  time  interval 
tesimal,  the  last  factor  in  the  bracket  becomes 
expression  for  D  is 


W0  lcr 
1  -k  t 
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becomes  infini- 
,  and  the  final 


D 


1 


kr 

2 


(167) 
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FOREWORD 


This  document  summarizes  the  information  contained  in 
URS  641-6,  Prediction  Models  for  Fire  Spread  Following  Nuclear 
Attack,  January  1965,  which  is  bound  separately. 
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Summary  Report 
of 

PREDICTION  MODELS  FOR  FIRE  SPREAD 
FOLLOWING  NUCLEAR  ATTACKS 


OBJECTIVE 

The  main  objective  of  the  study  is  to  develop  an  improved 
method  for  predicting  the  extent  of  fire  spread  beyond  the  area 
directly  ignited  by  nuclear  explosion  and  to  specify  the  data 
required  for  use  with  the  method  and  ways  of  acquiring  these 
data.  Since  the  ultimate  purpose  is  assessment  of  fire  damage 
to  resources,  the  method  also  includes  general  procedures  for 
deriving  damage  information  from  fire  spread  information. 

The  study  has  been  aimed  at  a  method  that  is  both  applicable 
with  the  currently  existing  data  and  capable  of  subsequent:  improve¬ 
ment  as  more  and  better  data  become  available. 

OUTPUT  REQUIREMENTS 

Estimates  of  fire  damage  to  resources  are  needed  primarily 
for  preattack  planning  and  possibly  for  early  postattack  appli¬ 
cations  and  can  be  performed  at  three  geographical  levels: 

i 

I  national,  regional,  and  local. 

Each  type  of  application  calls  for  a  different  level  of 
detail  in  the  output  information.  National  application  seems 
to  require  only  gross  information,  for  example,  the  average  loss 
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from  a  large  number  of  weapons.  For  regional  use,  a  somewhat 
greater  level  of  detail  is  needed.  Local  use  would  require  the 
highest  level  of  detail. 

DESCRIPTION  OF  THE  PROBLEM 

Mass  fires  are  characterized  by  rate  of  spread,  heat  and 
flame  intensity,  size  of  fire  burning  at  a  given  time,  final  burn¬ 
out  area,  and  percentage  of  burned  fuel.  Due  t:o  the  intimate 
relations  between  these  characteristics,  the  final  burnout  area 
(required  for  estimating  fire  damage)  cannot  be  easily  singled 
out  as  a  subject  of  investigation.  Thus  the  prediction  method 
is  likely  to  involve  a  representation  of  the  fire  in  time  and  in 
space  and  require  data  on  rate  of  spread,  burning  time  of  fuel, 
igniting  capability,  and  ignition  susceptibility. 

Factors  that  can  influence  fires  are  called  fire  spread 
variables  and  are  grouped  into  three  general  classes:  Fuel  vari¬ 
ables,  fire  spread  variables,  and  topography  variables. 

Fuel  variables  include  fuel  types ,  fuel  build-up  or  thick¬ 
ness  of  the  fuel  layer,  fuel  density  or  fraction  of  total  area 
occupied  by  the  fuel,  spatial  distribution  of  fuel,  fuel  age, 
and  fuel  fineness. 

Weather  variables  are  moisture  content  of  fuel,  wind  veloc¬ 
ity  and  direction,  and  temperature. 

The  topography  group  of  variables  consists  of  slope,  relief, 
and  possibly  altitude. 
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The  effect  of  a  variable  may  be  real  but  indiscernible  among 
other  more  important  ones,  or  it  may  remain  sensibly  constant 
over  the  usual  range  of  fire  spread  conditions. 

METHODS  OF  APPROACH 

A  purely  theoretical  approach,  relying  extensively  on  physi¬ 
cal  laws,  is  considered  to  be  impractical  in  the  present  stage  of 
knowledge.  A  purely  empirical  method  could  be  used  to  predict 
fire  spread  by  locating,  in  the  body  of  data  on  pa^r  fires,  a  set 
of  conditions  which  matches  the  future  situation  and  by  using 
the  past  to  predict  the  future.  However,  without  some  theoreti¬ 
cal  back-up,  it  would  be  difficult  to  untangle  the  numerous  fac¬ 
tors  and  efficiently  organize  the  empirical  data.  Thus  a  semi- 
empirical  approach  seems  to  be  most  appropriate.  In  this  approach 
a  mathematical  model  describes  various  characteristics  of  fire 
spread  in  terms  of  some  parameters  to  be  evaluated  from  empirical 
data . 


Whether  the  model  (and  its  output  information)  is  deter¬ 
ministic  or  stochastic  depends  largely  on  how  fuel,  weather,  and 
topography  conditions  are  specified.  Although  for  most  cases 
they  are  all  specified  in  a  stochastic  manner,  simplicity  is 
gained  by  treating  weather  conditions  as  deterministic. 

MATHEMATICAL  MODELS 


Two  families  of  mathematical  models  have  been  developed: 


€1  C* _ _  - - 
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have  been  derived  from  the  parent  model.  Specific  versions  serve 
a  variety  of  purposes:  illustrating  particular  features  of  the 
parent  model,  deriving  some  general  properties  which  could  not  be 
conveniently  studied  with  the  parent  model,  providing  more  con¬ 
venience  in  specific  applications.  Table  1  lists  various  models 
with  their  main  characteristics  and  applications. 

Table  1 

SYNOPSIS  OF  MODELS 

Model _ Characteristics  _  Appropriate  Fuel 

FIRE  FRONT  MODELS 

I  Stochastic.  One-dimensional.  Param-  Preferably  urban 

eters  are  \  and  p£;  the  latter  not 

available;  required  data  difficult  to 
obtain. 

IA  Same  as  I  but  with  more  detailed  input  Preferably  small 
information  on  fuel  characteristics.  or  moderate  urban 

IB  Same  as  I  but  with  provision  for  long-  As  above 
range  spotting  possibility. 

FUEL  STATE  MODELS 

II  Stochastic.  Two-dimensional.  Param-  General 

eters  are  A  and  pg.  Data  for  evalu¬ 
ating  pc  not  available. 

o 

IIA  Same  as  II  but  one-dimensional.  General,  prefer¬ 

ably  urban 

IIB  Same  as  II  but  deterministic  and  one-  General,  prefer - 
aimensional.  Existing  data  for  eval-  ably  rural 
uating  parameters  t,  k,  W0  are  incom¬ 
plete, 

IIC  Same  as  II  but  deterministic.  Appli-  Rural  and  urban 

cation  uses  existing  data  supplemented 
by  some  subjective  estimates. 
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Fire  Front  Models 


The  entire  fire  area  or  fire  front  is  pictured  as  a  random 
walker  moving  along  a  strip  of  fuel  area  divided  into  small 
square  sections  called  cells.  At  any  moment,  one  of  three  pos¬ 
sible  events  may  take  place:  either  the  fire  front  dies  (is 
burned  out)  or  it  moves  one  cell  forward  or  it  stays  where  it 
is.  The  occurrence  of  these  events  is  governed  by  their  rela¬ 
tive  probabilities,  which  are  determined  by  fire  spread  variables, 
fuel  density,  distribution  and  type,  wind  velocity,  relative 
humidity,  .  .  .  combined  into  two  main  parameters  \  and  p£. 

These  two  parameters  are  defined  as  the  probability  rates  (or 
probability  per  unit  time)  of  spreading  and  of  dying,  so  that 
within  the  time  increment  dt,  the  probability  of  spreading  to  the 
next  cell  is  \dt  and  that  of  stopping  permanently  is  p^dt.  It 
can  be  shown  that  \  is  numerically  equal  to  the  mean  rate  of 
spread  in  cells  per  unit  time  and  l/p£  is  equal  to  the  mean 
lifetime  of  the  fire  as  a  whole.  These  interpretations  provide 
convenient  methods  for  deriving  X  and  p^  from  observed  data,  pro¬ 
vided  such  data  exist.  Mathematical  relations  based  on  this 
mechanism  give  the  probability  of  finding  the  fire  still  burning 
and  that  of  finding  the  fire  burned  out  at  a  given  location  and 
a  given  time.  In  addition, the  mean  rate  of  spread  in  terms  of 
cells  per  hour  is  readily  calculated. 

With  slight  modifications,  the  mechanism  can  be  applied  to 
a  strip  of  fuel  area  in  which  the  cells  differ  significantly 
with  respect  to  fuel  characteristics  or  to  cases  in  which,  by 

t-Ko  -r\y~r\n  a  c  c  r\  f  Inno.rariop  cnn^finff.  t*hp  fi'rp  r>3n  -i  1  min  nvpr  nnP 
v  w  ~  -*o  ■ — "o~  wt - “oi - - - j - r  -  -  —  - 

or  more  cells. 
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Fuel  State  Models 


These  models  stochastically  or  deterministically  describe 
the  events  in.  each  cell  of  the  fuel  area.  At  any  time,  a  given 
cell  is  in  one  of  three  states:  the  unignited  state, including 
the  early  phase  in  which  the  cell,  though  ignited,  is  still 
incapable  of  igniting  other  cells;  the  flaming  state, in  which 
the  burning  intensity  is  sufficient  to  ignite  other  cells;  and 
the  burnout  state.  Whether  and  when  a  cell  passes  from  one  state 
to  the  next  depends  on  the  particular  set  of  fuel,  weather,  and 
topography  variables.  The  mathematical  formulation  however  does 
not  involve  these  variables  directly.  All  transitions  and 
derived  quantities  are  expressed  in  terms  of  two  parameters,  A 
and  ps, which  are  defined  in  a  manner  similar  to  X  and 
described  above.  They  can  be  analytically  related  to  the  burn¬ 
ing  time  of  the  cell  (not  of  the  entire  fire)  and  to  the  mean 
rate  of  spread  in  cell  widths  per  unit  time. 

The  same  mechanism  gives  rise  to  four  models,  each  requir¬ 
ing  a  slightly  different  method  of  implementation  and  serving 
a  specific  purpose. 

In  the  linear  version,  an  unignited  cell  may  pass  to  the 
burning  state  with  probability  Adt  per  time  increment  dt,  provided 
the  preceding  cell  is  in  F  state.  The  two-dimensional  version 
computes  this  transition  probability  by  taking  into  account  the 
possibility  of  ignition  from  more  than  one  neighboring  cell. 
Either  version  can  handle  the  situation  in  which  the  fire  can 
jump  over  one  or  more  cells, provided  A  values  can  be  determined 
according  to  jump  distances. 


URS  641-6 


7 


For  the  two-dimensional  version,  no  analytical  solution  can 
be  obtained  except  for  the  probability  of  any  cell  being  burned 
out  at  time  equal  to  infinity,  from  which  the  mean  number  of 
burnout  cells  or  the  mean  spread  area  is  readily  computed.  The 
linear  version,  however,  gives  several  interesting  analytical 
results.  The  rate  of  spread  is  A  cell  widths  per  unit  time  (as 
in  the  linear  model).  The  mean  burning  time  of  cells  is  1/u 

O 

time  units.  The  mean  final  spread  distance  is  infinite  if  A/p 

O 

is  equal  to  or  greater  than  1.  For  A/ps  less  than  1,  this  dis¬ 
tance  is  finite  and  can  be  calculated. 

The  fuel  state  model  includes  a  deterministic  version  which 
is  readily  applicable  to  normally  continuous  fuel  in  wildland 
areas.  The  rate  of  spread  is  assumed  to  be  proportional  to  the 
total  heat  flux  (with  embers  and  firebrands)  emerging  from  the 
burning  zone.  It  is  further  assumed  that  the  heat  supply  per 
unit  burning  area  is  independent  of  distance  from  the  front. 

A  cut-off  is  made,  however,  when  the  width  of  the  fire  exceeds  a 
certain  value  W  .  There  are  three  parameters:  WQ,  the  burning 
time  t  of  the  fuel,  and  a  spread  constant  k. 

As  in  the  stochastic  model,  the  final  spread  distance  is 
infinite  for  kT  >  1  and  finite  for  kT  <  1,  in  which  case  it  is 
easily  calculated.  The  model  also  gives  the  rate  of  spread  and 
fire  width  as  a  function  of  time. 

Finally,  a  practical  method  is  set  up  for  applying  the  fuel 
state  model  with  parametric  data  that  are  currently  available 
or  can  be  reasonably  estimated.  It  consists  of  two  steps. 

First,  it  predicts  whether  fire  can  spread  through  the  area 
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under  the  specified  conditions.  For  rural  fuels,  spread  capabil¬ 
ity  is  largely  determined  by  weather,  since  fuel  distribution  is 
mormally  continuous.  For  urban  fuels,  the  controlling  factor  is 
fuel  density.  Weather  conditions  for  spread  in  rural  fuels  and 
fuel  conditions  for  spread  in  urban  fuel  are  estimated  from  some 
limited  data. 

Second,  if  fire  spread  is  predicted,  its  progress  is  fol¬ 
lowed  from  cell  to  cell  and  from  one  time  increment  to  the  next 
by  means  of  rate-of-spread  data  until  weather  conditions  or  fuel 
conditions  change. 

PARAMETERS 

Each  parameter  is  a  complicated  quantity  representing  the 
combined  effects  of  several  basic  variables.  Its  evaluation, 
therefore,  requires  extensive  empirical  data.  The  parameter  Uf 
could  be  derived  from  data  on  the  mean  lifetime  of  the  fire. 
Unfortunately  no  such  data  can  be  found.  Evaluation  of  u 

b 

requires  data  on  burning  time  of  fuel.  Some  estimates  do  exist 
but  are  not  considered  sufficiently  accurate  for  the  purpose. 

A  and  X  are  directly  given  by  the  mean  rate  of  spread.  Informa¬ 
tion  on  rate  of  spread  exists  in  fair  amounts  for  rural  fuel 
but  is  rather  limited  for  urban  fuel. 

RATE  OF  SPREAD  DATA 

Statistical  analysis  of  rate  of  spread  data  provided  by 
the  Forest  Service  gives  relationships  between  rate  of  spread 
and  basic  variables. 
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For  .ural  fires,  rate  of  spread  increases  rather  markedly 
with  wind  velocity  for  forward  wind,  less  markedly  for  lateral 
wind,  and  generally  decreases  with  increasing  wind  velocity  for 
back  wind.  For  wind  velocity  from  0  to  10  mph,  rate  of  spread 
expectedly  decreases  with  increasing  relative  humidity.  How¬ 
ever,  it  becomes  independent  of  relative  humidity  as  wind  veloc¬ 
ity  increases  beyond  10  mph.  A  table  ii.  Part  II  gives  rate  of 
spread  values  for  various  combinations  of  wind  velocity,  wind 
direction,  and  relative  humidity. 

For  urban  fires,  effects  of  exterior  humidity  and  wind  can¬ 
not  be  evaluated  quantitatively.  Strongest  influence  seems  to 
come  from  fuel  density  and  fuel  type.  An  important  result,  at 
first  startling,  is  that  observed  rate  of  spread  increases 
rapidly  as  fuel  density  decreases.  The  explanation  is  that  for 
fires  which  do  not  go  out  in  low-density  fuel,  the  relative 
tendency  for  ground  spread  decreases  much  faster  than  the  tend¬ 
ency  for  spotting. 

VARIOUS  ESTIMATED  DATA 

The  probability  of  fire  spread  from  one  structure  to  another 
in  a  random  fuel  area  is  estimated  from  very  limited  data.  The 
approximate  curve  obtained  shows  that  at  distances  up  to  about 
20  ft,  fire  spread  is  certain  in  every  case,  and  at  300  ft  it 
hardly  ever  takes  place.  Until  more  detailed  information  is 
available  —  for  example,  one  curve  for  each  type  of  fuel  area  — 
this  curve  is  needed  for  estimating  fire  spread  vulnerabilities 
of  resources  and  fuel  requirements  for  spread  or  no-spread. 
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More  or  less  subjective  estimates  on  weather  conditions  for 
no-spread  in  rural  fuel,  weather  conditions  for  fire  extinction, 
and  burning  times  of  fuel  are  presented. 

FIRE  SPREAD  VULNERABILITY  OF  RESOURCES 

General  procedures  are  given  for  evaluating  the  chance  of 
fire  destruction  for  resources  located  in  areas  swept  by  fire. 
Distinction  is  made  of  three  types  of  environments:  urban,  rural, 
and  mixed.  The  procedures  also  depend  on  the  specified  details 
of  the  environment. 

WEATHER  DATA  FOR  FIRE  SPREAD  PREDICTION 

A  method  is  worked  out  for  obtaining  a  set  of  weather  con¬ 
ditions  that  is  consistent  with  the  meteorological  character¬ 
istics  of  .the  area  under  consideration  and  with  the  fact  that 
the  attack  may  take  place  at  any  time  in  the  future.  Depending 
on  the  available  effort,  the  weather  is  predicted  on  a  statisti¬ 
cal  basis,  or  an  average  weather  pattern  is  selected  by  visual 
inspection  of  weather  records,  or  a  worst  and  a  best  weather 
patterns  are  selected  to  obtain  an  upper  and  lower  limits  of 
fire  damage. 

LOCAL  APPLICATION  OF  FIRE  SPREAD  MODELS 

A  brief  discussion  is  devoted  to  the  use  of  the  fire  spread 
models  for  cities  of  large  or  moderate  size  and  to  the  problem 
of  data  collection  at  this  geographical  scale. 
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CONCLUSIONS  AND  SUGGESTIONS 


Major  accomplishments  of  this  program  are  summarized,  and 

|  suggestions  for  obtaining  certain  important  empirical  data  are 

?  given. 


